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The  principal  objectives  of  our  research  work  (Grant  Number  AFOSR-87-0342)  have  been 
two-fold.  First,  we  have  obtained  the  time  -  dependent  behavior  of  electron  swarms  in  various 
atomic  and  molecular  gases  using  a  novel  numerical  procedure  that  we  recently  developed.  The 
electrons  injected  in  the  gaseous  source  excite  and  ionize  the  atoms  as  well  as  molecules  and,  more 
importantly,  cause  dissociation  and  vibrational  excitation  of  the  molecular  species.  Second,  we 
have  calculated  the  cross  sections  and  the  rates  )f  negative  ion  production  via  the  process  of 
dissociative  electron  attachment  to  various  light  molecules.  We  have  also  investigated,  in  detail, 
the  role  played  by  initial  rovibrational  excitation  of  the  molecule  in  enhancing  the  rates  of 
production  of  negative  ion  beams.  Since  the  processes  of  dissociative  attachment  and  of 
vibrational  excitation  are  complementary  in  nature  (both  proceed  via  the  formation  of  an 
intermediate  resonant  anion  state),  we  have  also  investigated  the  isotope  effect  for  the  vibrational 
excitation  of  molecular  hydrogen  and  its  five  isotopes.  As  part  of  this  project  we  have  discovered  a 
very  useful  scaling  law  for  these  excitation  cross  sections.  Using  this  scaling  property  it  is 
possible  to  obtain  the  cross  sections  for  vibrational  excitation  or  deexcitation,  by  electron  impact, 
of  heavier  isotopes  of  from  the  corresponding  cross  sections  for  molecular  hydrogen. 

During  three  years  of  our  investigations  we  made  the  following  specific  advances: 

(a)  Time  Dependent  behavior  of  Electron  Swarms  in  Various  Gases. 

With  the  ultimate  aim  of  investigating  the  temporal  evolution  of  the  electron  velocity 
distribution  as  it  attains  its  equilibrium  in  a  hydrogen  source,  we  have  been  devising  stable 
numerical  algorithms  and  developing  corresponding  computer  codes  for  obtaining  the  electron 
velocity  distribution  in  a  general  gaseous  source.  It  should  be  remarked  that  the  procedure  which 
we  have  developed1,2  is  general  enough  in  that  it  can  be  easily  adapted  for  obtaining  velocity 
distribution  of  positron  swarms  or  ion  swarms  in  a  gas  as  long  as  the  relevant  collision  cross 
sections  are  available.  In  order  to  obtain  the  electron  velocity  distribution  majority  of  previous 
investigators  have  attempted  to  solve  the  Boltzmann  equation  either  analytically  or  numerically. 
Analytical  solutions  of  the  Boltzmann  equation  have  made  use  of  simple  model  collision  cross 
sections  which,  though  reasonable,  are  not  necessarily  veiy  accurate.  The  use  of  actual  cross 
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sections,  either  experimentally  observed  or  theoretically  calculated,  in  the  Boltzmann  equation 
forces  one  to  apply  numerical  techniques  for  its  solution.  The  procedure  used  most  frequently  in 
the  past  is  to  expand  the  angular  dependence  of  the  distribution  function  as  a  power  series  in  the 
Legendre  polynomials,  truncate  the  series  after  a  few  (typically  two)  terms  and  solve  the  resulting 
set  of  equation  numerically.  This  expansion  procedure  has  some  inherent  problems.  For  example, 
for  high  values  of  E/N  (the  ratio  of  the  applied  electric  field  to  the  gas  number  density)  a  large 
number  of  terms  in  the  Legendre  expansion  may  contribute  to  the  distribution  function.  A  large 
number  of  terms  in  the  Legendre  expansion  may  also  be  necessary  at  low  values  of  E/N  for  a  gas 
with  relatively  high  inelastic  cross  sections.  Inclusion  of  a  large  number  of  terms  in  the  expansion 
can  make  this  procedure  computationally  expensive  and  sometimes  numerically  unstable. 
Furthermore,  no  systematic  study  of  the  convergence  behavior  of  the  expansion  techniques  for  an 
arbitrary  gas  is  available.  We  felt  that  it  is  desirable  to  have  alternative  procedures,  especially  the 
ones  which  do  not  involve  any  expansion  of  the  distribution  function,  for  obtaining  the  electron 
velocity  distribution  function  in  a  gas  mixture. 

The  starting  point  of  our  procedure  is  to  realize  that  in  the  absence  of  collisions  with  the 
ambient  gas  particles  the  electron  swarm  is  accelerated  constantly  by  the  applied  electric  field  E  in 
such  a  manner  that  in  a  time  interval  At  the  velocity  of  each  electron  is  changed  by  an  amount  Av  = 
-eE/m.  Thus  the  difference  between  f(v,r),  the  distribution  function  at  time  f,  and  f(v  +  Av,t  + 

A f),  the  distribution  function  at  the  later  time  f + At,  must  be  due  to  collisions  of  electrons  with  gas 
particles.  (The  frequency  of  electron-electron  collisions  is  relatively  small.)  Thus  one  can  write 
f{v  +  Av,  t  +  At)  =  /(v,f)  +  R(\,t)  At,  (1) 

where  R(v,t)  is  the  collision  term  containing  all  the  relevant  cross  sections.  The  standard 
Boltzmann  equation  is  obtained  by  simply  expanding  the  left  hand-side  of  Eq.  (1)  to  first  order  in 
At  and  then  taking  the  limit  At  0.  Thus,  for  the  purposes  of  obtaining  the  electron  velocity 
distribution  function  Eq.  (1)  is  more  fundamental  than  the  standard  Boltzmann  equation  and  is 


certainly  computationally  more  convenient  than  the  Boltzmann  equation. 


Now  choosing  the  z  axis  to  be  along  the  applied  electric  field  E,  the  electrons  have  a 
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constant  acceleraaor  along  the  z  direction  and  zero  acceleration  in  the  transverse  direction.  Thus, 
in  Cartesian  coordinates  Eq.  (1)  is  rewritten  as 

flvx,  vy)  v7  -  (eE/m)At,f  +  At]  =  /(vx,  vy,  vz,  t)  +  R(y,t) At  (2) 

Because  of  die  axial  symmetry  along  vz,  f  is  stored  in  a  two-dimensional  array  only  as  a  function 
of  vz  and  vx  with  the  step  size  Av  given  by  Av  =  (eE/m)Ati  The  temporal  evolution  of  the 
distribution  function  then  involves,  in  pan,  a  shifting  of  the  array  containing  f  along  vz  at  each  time 
interval  At  Tins  shifting  of  the  array  essentially  achieves  the  same  effect  as  the  acceleration  of  the 
electrons  due  to  the  applied  electric  field  without  accumulating  any  round-off  error.  This  shifting 
procedure  in  itself  greatly  improves  the  numerical  stability  in  obtaining  the  electron  distribution 
function  as  compared  to  the  alternative  approach  of  evaluating  derivatives  numerically  in  the 
traditional  solution  of  the  Boltzmann  equation.  A  criterion  for  numerical  stability  is  that,  for  each 
value  of  Ivl,  die  time  step  At  should  be  chosen  small  enough  such  that  the  distribution  function  is 
larger  than  the  corresponding  collision  term  R(v,f)Af,  that  is, 

f(v,r)  >  R(v,f)Af  for  all  v.  (3) 

By  meeting  this  criterion  one  can  avoid  generating  physically-meaningless  negative  values  of  the 
electron  velocity  distribution  function.  Interestingly,  the  standard  Courant-Lcwy-Friedrichs 
condition,  namely,  that  the  step  sizes  Av  and  At  must  satisfy  the  inequality: 

(eE/m)Af  <  Av,  (4) 

which  is  a  necessary  condition  for  the  numerical  stability  of  a  partial  differential  equation  is  not 
necessary  in  our  case. 

A  step-by-step  procedure  for  attaining  the  equilibrium  velocity  distribution  on  the  computer 
is,  then,  as  follows: 

Step  1.  Start  with  a  given  distribution  function  at  time  t  which  is  stored  in  a  two-dimensional  array 
f(vx>vz)-  ft  could  either  be  an  analytical  function  (for  example,  a  Maxwellian  or  a 
Druyvesteyar.  or  a  delta  function)  representing  die  distribution  function  at  the  starting  time  I 
=  0  or  be  a  numerically  generated  function  at  some  earlier  time  f0. 
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Step  2.  Compute  the  collision  terms  R(vx,v7J  for  each  value  of  vx  and  vz.  Multiply  the  collision 
terms  by  Ar  and  add  to  the  distribution  function  from  which  they  were  obtained  [see  Eq. 
(2)]. 

Step  3.  Shift  the  resulting  array  along  the  vz  index  [f(vx,v7)  —  f(vx>  vz  +  Av)]  to  obtain  the  new 
distribution  function  at  the  later  time  t+  At 

Step  4.  Calculate  various  swarm  parameters  corresponding  to  the  later  time  r+ At  using  the  new 
distribution  .'unction.  Repeat  from  step  1  unless  the  swarm  parameters  stop  changing  in 
time  which  indicates  that  the  equilibrium  has  been  reached. 

We  have  used  this  procedure  for  obtaining  the  average  energy,  drift  velocity  and  ionizauun 
rate  for  electron-neon1  and  electron-argon2  systems  for  various  values  of  E/N.  The  equilibrium 
values  of  various  swarm  parameters  for  electrons  in  gaseous  neon  for  a  few  different  values  of  E/N 
are  given  in  Table  I  below. 

TABLE  I 


E/N 

Ionization  rate 

Drift  velocity 

Average  energy 

(in  Td) 

(in  107  s-l) 

(in  I07cms'1) 

(in  eV) 

141 

3.0074 

3.1110 

15.298 

190 

4.7470 

3.8049 

17.176 

237 

6.7463 

4.4547 

19.013 

283 

9.1608 

5.0139 

20.962 

353 

12.979 

5.8074 

23.787 

424 

17.113 

6.5471 

26.657 

495 

21.437 

7.2327 

29.573 

566 

26.013 

7.8715 

32.610 

This  unique  procedure,  that  we  have  described  above  for  obtaining  the  electron  energy 
distribution  in  a  gas,  is  computationally  stable  and  analyti  rally  exact.  In  particular,  this  procedure 
provides  a  detailed  time  evolution  of  various  swarm  parameters,  such  as  the  average  electron 
energy  or  the  rate  of  ionization  of  ambient  gas  particles  by  electron  impact  etc.,  whose  final 
equilibrium  values  can  be  experimentally  determined.  These  final  equilibrium  values  depend  only 
on  the  ratio  of  the  applied  electric  field  to  the  gas  number  density.  E/N.  and  are  independent  of  the 
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inidal  electron  energy  distribution  function.  Recently,  we  also  used  our  innovative  procedure  for 
obtaining  the  positron  velocity  distribution  function  in  various  rare  gases3. 

The  ultimate  aim  of  our  investigations  is  to  calculate,  in  the  future,  the  electron  energy 
distribution  in  a  realistic  mixture  of  atomic  and  molecular  hydrogen  using  all  the  theoretical  tools 
that  we  have  been  developing  and  numerically  testing  for  noble  gases.  For  this  purpose  we  will  be 
requiring  a  critical  compilation4  of  data  pertaining  to  cross  section  for  elementary  processes  in 
atomic  and  molecular  hydrogen.  In  the  case  of  noble  gases,  however,  we  have  discovered  that  the 
transient  behavior  of  any  swarm  parameter  as  it  evolves  to  its  final  equilibrium  value  depends  on 
the  initial  conditions  even  though  the  final  equilibrium  value  that  the  parameter  attains  is 
independent  of  the  initial  conditions.  -ave  further  discovered  that  this  transient  behavior  of  any 

swarm  parameter,  A(t),  can  be  accurately  fitted  by  a  sum  of  two  exponentially  decaying  factors  as 
follows: 

A(t)  =  c,exp(-a1t)  +  c2exp(-a2t).  (5) 

The  decay  constants  oq  and  are  determined  by  the  applied  electric  field.  We  have  noted  that  a 
single  exponential  function  does  not  provide  a  good  fit.  This  observation  indicates  that  perhaps  the 
initial  electron  swarm  could  be  regarded  as  consisting  of  two  components  with  each  component 
having  its  own  decay  time  to  reach  equilibrium  independently.  Figure  1  shows  the  temporal 
evolution  of  various  swarm  parameters  for  electrons  in  gaseous  neon,  for  E/N  =  566  Td.,  for  two 
different  initial  conditions.  In  this  figure,  'diamonds'  represent  the  actually  calculated  points  while 
the  solid  lines  are  the  two-exponential  fits  based  on  the  above  relation.  The  fitting  is  indeed 
superb.  Note  that  the  decay  constants  oq  and  are  the  same  for  all  swarm  parameters  for  a  given 
set  of  initial  conditions. 

At  the  present  moment  we  are  attempting  to  generalize  our  computer  code  to  include  (i) 
anisotropic  scattering  of  electrons  by  target  gas  atoms  via  differential  cross  sections,  (ii)  a  mixture 
of  gases  since  a  realistic  calculation  of  electron  swarm  parameters  in  a  nydrogen  source  must 
include  a  mixmre  of  atomic  and  molecular  hydrogen,  and  (iii)  electron  swarms  subjected  to  time  - 
dependent  electric  fields. 


c’/NEON 
E/N  «  566  Td 
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i.4 


4.6 


C  #  .  |  | 

E/N  ■  566  Td  1 


1  *r>  *^a  tr<*ncipnt  Jv*h»vinr  of  plf*cn*on  swnrm  nnrnrneters  in  neon  using  two 

i  i^uiv  nimij'  uv  1*1  iv  iuv  i  — 

exponentially  decaying  factors  (E/N  =  566  Td). 
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(b)  Nonlocal  effects  in  dissociative  electron  attachment  to  H2  and  its  isotopes. 


The  production  of  negative  ions  via  the  process  of  dissociative  electron  attachment  to 
molecular  hydrogen  (and  its  isotopes),  e*  +  Hj  —  H  +  H\  is  known  to  occur  when  the  incoming 
electron  and  the  neutral  target  molecule  form  an  intermediate  resonant  anion  state,  which  has  the 
ability  of  decaying  by  autodetaching  the  electron.  The  motion  of  the  nuclei  in  this  resonant  state  is 
governed  by  a  wave  function  q(R)  which  satisfies  an  integrodifferential  equation  involving  a 
complex,  nonlocal  potential.  In  solving  this  equation  one  can  make  use  of  a  local  approximation  in 
which  one  assumes  that  the  set  of  vibrational  levels  which  are  accessible  for  a  given  incident 
electron  energy  can  be  regarded  as  complete.  While  this  approximation  may  yield  acceptable 
values  for  the  cross  sections  well  above  threshold,  it  is  expected  to  be  less  accurate  near  threshold 
where  there  are  fewer  energetically  open  vibrational  channels.  During  the  present  investigations 
we  solved  the  full  nonlocal  integrodifferential  equation  for  the  nuclear  wave  function  q(R)  near 
threshold  using  a  newly  developed  technique  which  involves  matrix  inversion.  This  investigation 
has  allowed  us  to  compare  the  resulting  nonlocal  attachment  cross  sections  with  those  that  utilize 
the  nuclear  wave  function  ^(R)  which  is  obtained  using  the  local  approximation  to  the  full 
integrodifferential  equation.  Besides  providing  an  assessment  of  the  effect  of  the  local 
approximation  on  the  electron  attachment  cross  sections,  the  present  calculations  will  also  provide 
cross  sections  for  H"  production  which  are  more  accurate  than  previously  calculated5.  The  nuclear 
wave  function  ^(R)  satisfies  the  following  integrodifferential  equation: 

[  Tfj(R)  +  V'(R)  -  E]$(R)  =  -V(E-EvJ(,R)xV|J|(R)  -  f  dR'  K(R,  R')  i;  (R')  , 

(6) 


with  the  kernel  K(R,  R’)  given  by 

V(e,R) 
e  +  iq  ’ 

(7) 

where  TN  is  the  kinetic  energy  of  the  nuclei,  V"  is  the  effective  potential  energy  in  which  they  move 
in  the  resonant  state  and  E  is  the  total  energy  of  the  interacting  system.  The  quantity  V(e,  R)  is 
t  a  erred  to  us  die  interaction  matrix  element,  it  represents  the  euupling  between  the  discrete  and 


K(R.R')  =  Ix‘j(R')Xvj(R) 

v.J 


lim 
i-  +o 

0 


de^Ml 


E-E 


vj- 
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continuum  states.  This  coupling  leads  to  an  energy  shift  of  the  resonance  (the  principal  part 
integral  on  the  right  hand  side  of  (7))  and  provides  a  width  to  the  resonance  (the  imaginary  part  on 
the  right  hand  side  of  (7))  which  determines  the  lifetime  of  the  resonance. 

In  Figure  2  the  nonlocal  cross  section  for  attachment  to  H2  in  the  J  =  0,  v  =  0  level  is 
compared  to  its  local  counterpan.  It  can  be  seen  that  the  difference  between  the  two  cross  sections 
at  a  given  energy  is  small,  on  the  order  of  10  percent  near  threshold  where  the  difference  is  largest. 
That  the  nonlocal  cross  section  is  larger  than  the  local  one  reflects  the  fact  that  in  the  local 
calculation,  all  vibrational  channels  are  considered  to  be  open,  whether  or  not  they  are  actually 
energetically  accessible.  The  two  cross  section  curves  evidently  do  not  merge  into  each  other  even 
at  higher  energies  because,  in  the  present  calculations,  excitation  of  the  continuum  levels 
(corresponding  to  H  +  H  +  e)  was  not  taken  into  account. 

The  most  significant  difference  between  the  local  and  nonlocal  cross  sections  is  the 
conspicuous  step  structure  in  the  nonlocal  cross  section  in  the  energy  range  below  about  4.5  eV. 
This  step  structure  could  not  be  clearly  seen  in  the  available  experimental  data6  since  the  energy 
resolution  of  the  apparatus  used  was  of  the  order  of  0. 1  eV.  These  steps  occur  at  energies  for 
which  a  new  vibrational  channel  opens  up.  Thus,  for  example,  the  first  step  in  the  cross  section 
for  attachment  to  H2  in  v  =  0,  J  =  0  level  at  about  3.83  eV  corresponds  to  the  opening  of  the  v  =  10 
channel,  while  the  last  is  at  4.39  eV,  where  the  v  =  13  level  becomes  energetically  accessible.  It 
might  be  tempting  to  attribute  the  loss  of  attachment  flux  at  these  steps  merely  to  diversion  into  the 
newly  opened  vibrational  channel,  but  Figure  3  shows  that  this  is  an  oversimplification.  The 
opening  of  a  new  vibrational  level  is  accompanied  by  upward  jumps  in  the  cross  sections  for  the 
vibrational  excitation  of  previously  open  channels.  The  magnitude  of  the  upward  jumps  becomes 
smaller  as  we  consider  levels  farther  removed  in  energy  from  the  newly  opened  level.  The  jumps 
are  such  that  ihe  sum  of  ail  discontinuities,  including  both  in  the  dissociative  attachment  and  in  the 
vibrational  excitation  cross  sections,  is  indeed  zero.  The  disappearance  of  the  step  stru  .ire  in  the 
local  treatment  is  not  very  surprising  since  the  local  approximation  entails  replacing  tiie  quantity  E  - 
Evj  in  the  denominator  on  the  right  hand  side  of  Eq.  (7)  by  the  incident  electron  energy  and 
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Figure  2.  Local  and  nonlocal  cross  sections  for  attachment  to  H2  in  the  lowest  rovibrational  level. 
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Figure  3.  Various  nonlocal  cross  sections  for  scattering  from  the  lowest  rovibrational  level  of  H2, 
including  that  for  dissociative  attachment  and  those  for  excitation  of  the  v  =  10  to  v  =  1 3  levels. 
The  cross  sections  for  vibrational  excitation  of  v  =  10, 1 1,  !2  and  13  levels  arc  multiplied  by 
factors  of  15,  25,  50  and  100,  respectively. 
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assuming  that  the  set  of  vibrational  levels  over  which  one  sums  can  be  regarded  as  complete; 
naturally,  in  such  an  approximation  all  sensitivity  to  the  number  of  open  channels  is  washed  out. 

Finally,  Figure  4  shows  the  peak  total  attachment  cross  sections  for  various  rotationless 
levels  of  H2  and  HD  against  electron  energy  on  a  logarithmic  plot.  We  have  discovered  a  tight 

upper  bound  on  the  attachment  cross  sections  in  terms  of  the  incident  electron  energy  e  = 

The  dot-dashed  curve  on  the  top  represents  1/k2,  which  is  seen  to  provide  an  upper 
bound  to  the  attachment  cross  sections  for  all  six  isotopes  of  H2  for  all  values  of  the  incident 

n 

electron  energies  considered.  This  is  consistent  with  an  observation  of  Gauyacq'  that  the 
dissociative  attachment  cross  section  is  of  the  form 

°DA  00  =  ~2  (k)  - 

k 

where  Pdc[  (k),  representing  the  probability  that  the  electron  does  not  autodetach  in  the  resonant 
state,  is  less  than  unity.  In  fact,  1/k2  provides  a  tighter  upper  bound  than  rc/k2.  Details  of  these 
investigations  will  soon  be  appearing8  in  Physical  Review  A. 

(c)  Isotope  effect  in  vibrational  excitation  of  molecular  hydrogen. 

The  process  of  vibrational  excitation  is  closely  related  to  the  process  of  dissociative 
attachment  in  that  it,  too,  proceeds  via  an  intermediate  resonant  anion  state.  Using  the  nonlocal 
equation  described  above,  we  have  also  obtained  cross  sections  for  vibrational  excitation  of 
molecular  hydrogen  via  the  2ZU+  resonance  for  electron  energies  up  to  5  eV.  Our  objective  in 
studying  vibrational  excitation  has  been  twofold:  to  compare,  as  with  dissociative  attachment,  the 
nonlocal  cross  sections  with  those  obtained  using  the  local  approximation  and  also  to  examine  the 
dependence  of  the  cross  sections  on  the  nuclear  mass  (that  is,  the  isotope  effect). 

Figure  5  shows,  as  a  typical  example,  cross  sections  for  excitation  from  the  ground  level 
(vj  =•  0)  to  the  vf  =  4  level  using  both  the  local  and  the  nonlocal  formulations  of  the  theory.  The 

nonlocal  effects  are  seen  to  be  small  here,  although  differences  of  up  to  a  factor  of  two  are  seen  for 
excitation  of  the  highest  vibrational  levels.  The  prominent  feature  here,  however,  is  the  well- 
defined  structure  in  the  cross  sections  which,  we  emphasize,  is  by  no  means  a  nonlocal  effect, 
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Figure  4.  Total  cross  sections  for  attachment  to  H2  (solid  curve)  and  HD  (dashed  curves)  in 
various  rotationless  vibrational  levels.  The  dot-dashed  curve  on  the  top  represents  an  upper  bound 
provided  by  1/k^  where  /i\"/2mc  is  the  energy  of  the  incident  electron. 


Figure  5.  Cross  sections  for  vibrational  excitation  (v  =  0  —  4)  of  I  I2  using  both  nonlocal  and  local 
formalisms  of  the  theory. 
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being  present  in  the  local  cross  section  as  well.  Tills  structure  is  seen  only  for  electron  energies 
below  the  threshold  for  dissociative  attachment  and,  for  excitation  from  the  ground  level,  is  entirely 
absent  for  Vf<  3,  becoming  more  pronounced  with  increasing  Vf.  This  behavior  is  in  qualitative 
agreement  with  experimental  observations9,  as  well  as  in  quantitative  agreement  with  theoretical 
results  reported  by  other  authors10.  We  attribute  the  peaks  in  the  cross  sections  to  the  existence  of 
singularities  in  the  scattering  amplitude  in  the  complex  energy  plane.  These  singularities  are 
associated  with  the  bound  vibrational  levels  of  the  resonant  anion  potential  energy,  and  coincide 
with  them  (in  energy)  in  the  limit  as  the  level  width  tends  to  zero.  We  have  confirmed  that  as  the 
width  is  made  smaller,  the  peaks  in  the  cross  section  become  sharper  and  occur  at  values  of  the 
electron  energy  which  approach  those  at  which  an  anion  bound  state  could  be  excited.  The  same 
structure  is  also  seen  in  excitation  from  higher  vibrational  levels,  and  it  occurs  there  for  lower 
values  of  vp,  appearing  in  the  superelastic  cross  sections  by  Vj  =  2.  The  exact  reason  for  the 
behavior  of  the  structure  as  a  function  of  Vf  is,  at  present,  under  investigation. 

We  have  also  obtained  cross  sections  for  excitation  from  various  initial  levels  to  higher 
vibrational  levels  for  all  heavier  isotopes  of  hydrogen,  namely,  for  HD,  D2,  HT,  DT,  and  T2.  As 
part  of  our  research  endeavors  we  have  discovered  a  very  useful  scaling  law  for  these  excitation 
cross  sections.  Using  this  scaling  propeity  it  is  possible  to  obtain  the  cross  sections  for  vibrational 
excitation  or  deexcitation,  by  electron  impact,  of  heavier  isotopes  of  IL,  from  the  corresponding 

cross  sections  for  molecular  hydrogen.  Now,  using  separately  the  potential  curves  of  linear 
harmonic  oscillators  as  well  as  of  Morse  oscillators  to  mimic  the  actual  potential  curves  of  the 
neutral  molecule  H2  and  its  resonant  state  H2",  we  have  derived  analytically  that,  in  the  impulse 
limit  of  a  resonance,  an  arbitrary  vibrational  excitation  cross  section  o(v,  —  vf)  is  proportional  to 
M  - 1  vi  -  vf  1/2  w|iere  M  js  (he  reduced  nuclear  mass.  [The  relative  masses  of  the  isotopes  H2,  HD, 
D2,  HT,  DT,  orT2are  1.00,  1.33,  1.50,  2.00,  2.40  and  3.00,  respectively.)  That  this  scaling  law 
is  obeyed  quite  well  by  vibrational  excitation  uuss  sections  (calculated  using  the  actual  potential 
curves  of  H2  and  H2‘)  can  be  seen  in  Figure  6  where  the  quantity 
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<*x(V|-  Vf) 

[Mxl 

°Hj(vi  —  0  ) 

KJ 

(8) 


(with  X  =  H2,  HD,  D2,  HT,  DT,  and  T2)  is  shown,  for  v;  =  2  and  incident  electron  energy  e  =  5 
eV,  as  a  function  of  Vf.  Our  calculated  excitation  cross  sections  are  in  quite  good  agreement  with 
this  scaling  behavior  for  all  values  of  Vj  and  vf  up  to  and  including  5.  This  is  because  the  resonant 
state  under  present  consideration,  namely  the  2IU+  state  of  H2’,  has  a  relatively  large  width  (or 
small  lifetime)  and,  therefore,  can  be  construed  as  the  impulse  limit  of  the  resonance.  It  is  rather 
remarkable,  as  seen  in  Figure  6,  that  the  above  isotope  scaling  law  for  vibrational  excitation  cross 
sections  of  molecular  hydrogen  is  aoplicable  for  cross  section  values  varying  over  several  orders  of 
magnitude!  We  intend  to  write  the  details  of  these  investigations  soon  in  the  form  of  a  paper  for 
publication  in  a  refereed  journal. 

Finally,  an  updated  list  of  all  the  publications  and  presentations  earned  out  under  the  tenure 
of  the  present  Grant  is  provided  in  the  appendices  A  and  B. 
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Appcndix  B:  Presentations 

During  the  tenure  of  the  Grant,  AFOSR  -  87  -  0342,  the  following  presentations  were  made  at 
various  national  and  international  meetings.  Abstracts  of  son, 2  of  these  presentations  are 
attached  to  the  present  Report. 

1.  "Formation  of  ground  and  excited  states  of  antihydrogen",  Bull.  Am.  Phys.  Soc.  22,  991 

(1988)  (with  Sultana  N.  Nahar);  presented  at  the  1988  annual  meeting  of  the  Division  of 
Atomic,  Molecular  and  Optical  Physics,  Baltimore,  Maryland,  April  18-20, 1988. 

2.  "Relativistic  approach  for  e  scattering  from  argon",  Bull.  Am.  Phys.  Soc.  32,  935  (1988) 

(with  Sultana  N.  Nahar);  presented  at  the  1988  annual  meeting  of  the  Division  of  Atomic, 
Molecular  and  Optical  Physics,  Baltimore,  Maryland,  April  18-20, 1988. 

3.  "An  exact  numerical  solution  of  the  time  dependent  Boltzmann  equation",  Bull.  Am.  Phys. 

Soc.  24, 295  (1989)  (with  P.J.  Drallos);  presented  at  the  41st  Annual  Gaseous  Electronics 
Conference,  Minneapolis,  Minnesota,  October  18-21, 1988. 

4.  "Particle  beam  production  via  dissociative  election  attachment  to  molecules",  presented  at  the 

S PIE's  OE/LASE'89,  Los  Angeles,  California,  January  15-20,  1989. 

5.  "High-field  time-dependent  positron  velocity  distribution  functions",  Bull.  Ain.  Phys.  Soc. 

24, 1409  (1989)  (with  P.J.  Drallos);  presented  at  the  1989  annual  meeting  of  the  Division 
of  Atomic,  Molecular  and  Optical  Physics,  Windsor,  Ontario,  May  17-19, 1989. 

6.  "Energetics  of  negative  ion  formation  via  dissociative  attachment  of  e  +  LiH".  Bull.  Am. 

Phys.  Soc.  24, 1401  (1989)  (with  H.H.  Michels);  presented  at  the  1989  annual  meeting  of 
the  Division  of  Atomic,  Molecular  and  Optical  Physics,  Windsor,  Ontario.  May  17-19, 
1989. 

7.  "Nonlocal  effects  in  dissociative  electron  attachment  to  H2",  (with  D.E.  Atems);  presented  at 

the  Sixteenth  International  Conference  on  the  Physios  of  Electronic  and  Atomic  Collisions, 
New  York,  July  26-  August  1,  1989. 

8.  "Temporal  evolution  of  electron  swarms  in  argon  and  neon"  (with  P.J.  Drallos);  presented  at 

the  Sixteenth  International  Conference  on  the  Physics  of  Electronic  and  Atomic  Collisions. 
New  York,  July  26-  August  1,  1989. 

9.  "An  exact  theory  for  the  transient  behavior  of  electron  swarm  parameters"  (with  P.J.  Drallos); 
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A  novel  algorithm  for  calculating  the  time  evolution  of  the  electron  energy 
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We  are  presenting  a  novel  numerical  technique  for  obtaining  the  time  evolution  of  the  electron 
velocity  and  electron  energy  distribution  functions  in  the  presence  of  a  uniform  electric  field. 

Using  this  technique,  the  various  swarm  parameters  can  be  evolved  for  sufficiently  long  times 
so  that  equilibrium  can  be  reached  without  incurring  any  numerical  instabilities.  Results  are 
presented  for  electron  swarms  in  gaseous  neon  for  various  values  of  E /N. 


In  nearly  all  aspects  of  gaseous  electronics,  the  electron 
energy  distribution  is  of  fundamental  importance.  Knowl¬ 
edge  of  the  electron  energy  distribution  function  ( EEDF)  is 
usually  gained  through  its  relation  with  the  electron  velocity 
distribution  function  (EVDF)  which  is  a  solution  of  the 
Roltzmann  equation  for  a  given  set  of  collision  cross  sec¬ 
tions.  The  Boltzmann  equation  can  be  solved  analytically  for 
an  EEDF  for  only  a  few  simple  cases.1  In  practice,  however, 
numerical  methods  have  to  be  used  for  its  solution.  In  tradi¬ 
tional  techniques  for  solving  the  Boltzmann  equation,  the 
distribution  function  is  expanded  in  the  Legendre  polynomi¬ 
als — either  a  two-term  expansion  or  a  multiterm  expan¬ 
sion— and  the  resulting  set  of  equations  arc  solved  numeri¬ 
cally  for  the  equilibrium  EEDF  after  making  some 
simplifying  approximations  for  the  collision  cross  sections. 

However,  these  expansion  techniques  do  have  some  in¬ 
herent  drawbacks.2  For  example,  the  two-term  expansion 
method  breaks  down  for  high  values  of  E/A' but  works  well 
for  low  E /N  values.  The  EEDF  under  high  E  /N  situations 
can,  in  principle,  be  obtained  by  taking  more  terms  in  the 
Legendre  expansion  of  the  Boltzmann  equation.  It  is  not 
always  clear  when  the  multiterm  expansion  is  to  be  preferred 
over  the  two  term,  nor  just  how  many  terms  in  the  expansion 
of  the  EEDF  are  to  be  included.  The  convergence  behavior 
of  a  multiterm  expansion,  especially  for  high  values  of  E  /N, 
is  also  not  completely  understood.  Furthermore,  the  cross 
sections,  which  arc  adjusted  to  reproduce  the  experimental 
swarm  parameters  in  a  two-term  expansion,  do  not  yield  the 
same  EEDF  or  the  corresponding  swarm  parameters  when 
used  with  a  multiterm  expansion  and  vice  versa.  It  would  be 
desirable,  then,  to  have  a  procedure  for  obtaining  the  equilib¬ 
rium  EEDF  and  equilibrium  EVDF  that  does  not  involve 
Legendre  expansion  of  the  distribution  function.  Here  we 
have  utilized  a  finite-difference  scheme,  previously  used  by 
Tagashira  and  co-workers,3,4  in  which  no  expansion  of  the 
EVDF  in  spherical  harmonics  is  needed  and  the  time  evolu¬ 
tion  of  die  EVDF  to  its  equilibrium  value  is  obtained  expli¬ 
citly.  Thus,  the  question  of  how  many  terms  should  be  taken 
in  the  expansion  of  the  EVDF,  which  is  the  solution  of  the 
Boltzmann  equation,  is  completely  avoided. 

In  this  communication  vve  present  a  novel  algorithm  for 
obtaining  the  time  evolution  of  the  EVDF  nuincricaiiy.  This 
is  accomplished  by  beginning  with  the  Boltzmann  equation 
which  can  be  written  as 

^-—1.  +  a-VJ"(v,/)  =  R(v,/),  (1) 


where  a  =  —  (eE/m )  is  the  acceleration  of  the  electrons  of 
mass  m  due  to  the  electric  field  E,  and  R(\,t)  is  the  collision 
term  containing  all  of  the  relevant  cross  sections.  Now  we 
multiply  both  sides  of  Eq.  ( 1 )  by  a  finite  time  interval  A /  and 
then  add J\v,t)  to  each  side  which  yields, 

/(v,r)  +  ^At  ~  +  =/(v,/)  +E(v,/)A/, 

(2) 

with  Av  =  aAr  The  terms  on  the  left-hand  side  of  Eq.  (2) 
can  be  combined  to  yield  the  final  result, 

/(v  +  Av,r  +  At)  =/(v,r)  +  R(v.t)At.  (3) 

Equation  ( 3 )  describes  the  time  evolution  of  the  EVDF 
and  can  thus  be  used  to  obtain  the  time  evolution  of  the 
EEDF  and  the  various  electron  swarm  parameters.  Note 
that  Eq.  (3)  is  entirely  equivalent  to  Eq.  ( 1 )  and  is  perhaps 
even  more  fundamental  than  Eq.  ( 1 ),  as  Eq.  (3)  is  a  neces¬ 
sary  step  that  one  must  go  through  in  a  textbook  derivation 
of  the  Boltzmann  equation.5  In  such  a  derivation,  one  would 
normally  expand  the  left-hand  side  of  Eq.  (3)  to  first  order 
in  At,  then  take  the  limit  A/— 0  to  obtain  the  Boltzmann 
equation  [  Eq.  ( 1 )  ] .  In  the  velocity  space,  the  collision  term 
R(y,t)  in  Eq.  (3)  is  difficult  to  evaluate  in  Cartesian  coordi¬ 
nates,  but  straightforward  in  spherical  coordinates  since  the 
collision  cross  sections  depend  only  on  the  electron  impact 
speed  v.  Thus,  it  might  seem  that  spherical  coordinates 
would  be  the  natural  choice  for  evaluation  of  Eq.  (3).  In 
their  calculations,  Kitamori,  Tagashira,  and  Sakai4  used 
spherical  coordinates  to  evaluate  the  term  a-VJjv,/)  of  Eq. 
( 1 ).  This  involves  derivatives  of / with  respect  to  v  and  0 
(polar  coordinates  in  velocity  space)  which  must  be  evaluat¬ 
ed  numerically.  This  procedure  leads  to  strong  numerical 
instabilities  in  the  EVDF  and  is  computationally  expensive. 

On  the  other  hand,  as  we  will  now  show, 
/(v  +  Av,/  +  At)  of  Eq.(3)  is  extremely  easy  to  evaluate  in 
Cartesian  coordinates.  The  acceleration  of  electrons  along  E 
(which  we  choose  to  be  along  the  z  direction)  is  constant, 
and  is  zero  in  the  transverse  directions.  Thus,  in  Cartesian 
coordinates,  Eq.  (3)  is  written  as 

f[vx,vy,u,  -  (eE/m)At,t  +  At  ] 

=J(vx,v>,,v„t)  +R(\,t)At.  (4) 

Equation  (4)  is  very  well  suited  forevaluation  on  a  com¬ 
puter.  Since  there  is  axial  symmetry  along  nz,/nccd  only  be 
stored  as  a  function  of  v.  and  vx  (orvy)  in  such  a  way  that  the 
velocity  increments  Av  satisfy  An  =  (eE/m)At.  Evaluation 
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£/Tt  =  144  Td 


tyN  =  S66Td 


TKC  (rtf«c.)  TMC  (nice.) 

FIG  1  The  lime  evolution  of  various  electron  swarm  parameters  in  gaseous  F:IG  3,  The  time  evolution  for  vanous  electron  swarm  parameters  in  gase- 

neon  for  E /N  ■*  144  Td.  ous  neon  for  E /.V --  566  I'd. 


of /in  Eq.  (4)  then  merely  involves,  in  part,  a  shifting  of  the 
array /( vx,vx)  along  v,  at  each  time  interval  Af.  This  shifting 
procedure  accomplishes  all  the  acceleration  effects  of  the 
electrons  due  to  the  electric  field  without  incurring  the  relat¬ 
ed  problems  of  numerical  instabilities  which  arise  mainly 
from  the  evaluation  of  derivatives.  The  evaluation  of  the  col¬ 
lision  term4  in  Eq.  (4)  involves  an  integral  over  the  polar 
angle  0  (in  the  velocity  space)  and  requires  a  knowledge  of 
the  distribution  function  at  various  values  of  v  and  0,  that  is, 
at  various  v-0  grid  points  in  the  velocity  space.This  integra¬ 
tion  can  be  carried  out,  even  though  /is  known  only  as  a 
function  of  vx  and  v„  by  simply  interpolating^  d^.d, )  to  get 

Thus,  the  novel  algorithm  for  evaluating  the  time  evolu¬ 
tion  of  the  EVDF  is  as  follows:  Starting  from  a  distribution 
function  at  some  time  I  (a  Maxwellian  at  t  =  0,  for  example) 
which  is  stored  in  a  two-dimensional  array/to,  ,v, )  such  that 
Au  =  ( cE/m)At ,  the  collision  terms  R(vx,vx)  for  each  vx 
and  vx  are  calculated.  These  collision  terms  are  then  multi¬ 
plied  by  At  and  added  to  the  distribution  function  from 
which  they  were  obtained  in  accordance  with  Eq.  (4). 


tyN  =  283  Td 
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FIG  2  The  time  evolution  for  larious  electron  swarm  parameters  m  gase¬ 
ous  neon  for  li/S'  —  ?S3  Td. 

5602  J.  Appl.  Phys..  Vol.  63,  No.  1 1 , 1  June  1 988 


The  resulting  array  is  then  shifted  along  the  v1  index 
\f(ux,vt  )—f(vx,vt  +  Ad))  and  it  becomes  the  new  distribu¬ 
tion  function  at  the  later  time  /  -F  At.  This  procedure  is  re¬ 
peated  while  various  swarm  parameters  are  calculated  from 
each  new  distribution  function  corresponding  to  a  new  time 
I  +  At.  Equilibrium  is  obtained  when  the  swarm  parameters 
cease  to  change  in  time. 

We  have  used  this  method  to  ootain  the  time  evolution 
of  various  electron  swarm  parameters  for  the  electron-neon 
system  for  various  values  of  E  /N.  In  all  of  our  calculations 
we  assumed  an  initial  Maxwellian  velocity  distribution  at 
/  =  0  and  a  gas  density  N  of  3.54xl0+l6  cm'1  or 
1.32X  10_5amagat  (ITorrat  273  K).  The  velocity  steps  Ad 
ranged  from  2.2  x  107  to  3.4  x  107  cm/s  and  the  tunc  steps  A t 
ranged  from  0.064  to0.032  ns  as  E /N  was  varied  from  144  lo 
566  Td.  We  point  out  that,  in  practice,  A t  may  be  chosen  to 


E/N«566Td 


MtTlAL  MAXWELLIAN  EVDF 


final  EQUUBRajM  evdf 


FIG.  4  The  initial  Maxwellian  and  final  equilibrium  electron  velocity  .  i>- 
tnbution  function  in  neon  for  E/N  —  566  Td. 
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be  a  proper  fraction  of  mCut/eE  and  only  this  fraction  of  each 
element  of  the  EVDF  is  shifted  per  tune  step.  Such  a  proce¬ 
dure  is  observed  to  enhance  the  numerical  stability.  Another 
criterion  for  choosing  A  t  is  that  the  term  R  ( \,t )  At  of  Eq.  (4 ) 
be  smaller  than'fv,/)  for  all  values  of  v.  Figures  1—3  display 
the  calculated  time-dependent  behavior  of  various  electron 
swarm  parameters  in  gaseous  neon  for  three  different  values 
of  E/N.  The  relevant  collision  cross  sections  were  taken 
from  Ref.  4.  The  equilibrium  values  of  these  swarm  param¬ 
eters  are  in  very  good  agreement  with  those  calculated  by 
Kitamori,  Tagashira,  and  Sakai.4  We  have  also  noted  that 
the  final  equilibrium  values  of  the  swarm  parameters  are 
unaffected  by  the  average  energy  value  of  the  initial  velocity 
distribution,  although  the  transient  behavior  may  be  some¬ 
what  different.  For  example,  an  overshoot  in  the  drift  veloc¬ 
ity  is  observed  if  the  initial  average  energy  of  the  distribution 
function  is  less  than  the  final  equilibrium  value  of  the  energy. 


Figure  4  shows  the  initial  and  final  EVDF  for  the  case 
E/N =  566  Td. 

The  fortran  77  code  which  we  used  was  less  than  400 
lines  in  length  and  required  about  500  CPU  seconds  to  reach 
equilibrium  for  E /N  —  566  Td  on  an  Amdahl  470/V 6  main¬ 
frame  computer. 
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Differential  and  integrated  cross  sections  for  the  formation  of  antihydrogen  by  the  impact  of 
intermediate-energy  (20-500  keV)  antiprotons  on  positronium  arc  calculated  using  the  first  Born 
approximation.  The  calculations  are  carried  out  for  the  formation  of  antihydrogen  in  ground  and 
various  excited  electronic  states  (n  =  1-3)  when  positronium,  the  target  atom,  is  in  the  ground  state, 
and  for  the  formation  of  antihydrogen  in  the  ground  state  when  the  positronium  is  in  various  excit¬ 
ed  electronic  states  (n  =  1-2).  The  l/nJ  behavior  for  the  capture  cross  sections  is  used  to  calculate 
the  total  (that  is,  all  states  added  together)  integrated  cross  sections.  The  cross  sections  for  the  for¬ 
mation  of  antihydrogen  presented  here  are  obtained  from  those  for  the  formation  of  positronium  by 
the  impact  of  positrons  on  hydrogen  atoms  by  using  charge  invariance  and  the  principle  of  detailed 
balance. 


I.  INTRODUCTION 

Because  of  the  recent  availability  of  an  antiproton 
beam  at  the  Low-Energy  Antiproton  Ring  (LEAR)  facili¬ 
ty  at  Organisation  Europcenne  pour  la  Recherche 
Nuclcairc  (CERN),  experiments  for  the  formation  of  one 
of  thc_  simplest  atoms  of  antimatter,  namely,  antihydro¬ 
gen  (H),  are  being  proposed  and  planned.1'  4  One  reason 
for  the  strong  interest  in  the  formation  of  antihydrogen  is 
its  relative  stability  compared  to  other  exotic  atoms  such 
as  muonium,  protonium,  positronium,  etc.,  which  makes 
it  more  suitable  for  carrying  out  experiments  with  antihy¬ 
drogen  for  various  diagnostic  purposes.  In  the  present 
work  the  cross  sections  are  calculated  for  the  formation 
of  H  by  positron  (c+)  capture  during  the  collisions  of 
intermediate-energy  (20-500  keV)  antiprotons  (p)  with 
positronium  (Ps)  using  the  first  Born  approximation 
(FBA).  Schematically  the  processes  are 

p  +  Ps(  Is)— *H( «/mi+ e  “  ,  (1) 

where  nlm  =  \s0,  2 s0,  2 p0,  2 p±l,  3s0,  3 p0,  3p . 3 d0, 
3rf±p  and  3 rf±2>  and 

p  +  Ps( tilin )— H(ls)+c“  ,  (2) 

where  nlm  =  ls0,  2s0,  2 p0,  and  2p±, .  The  1  /n 5  behavior 
for  the  capture  cross  sections  in  the  Born  approximation 
is  used  to  estimate  the  differential  and  integrated  cross 
section  for  the  formation  of  higher  excited  states  of 
antihydrogen  starting  from  the  ground  state  of  Ps  as  well 
as  for  the  formation  of  ground  state  of  II  starting  from 
the  higher  excited  states  of  Ps. 

II.  THEORY 

Following  earlier  works,2,3  the  cross  sections  (o-j|)  for 
the  formation  of  II  by  the  impact  of  p  on  Ps  of  Eqs.  (1) 
and  (2)  are,  by  charge  invariance,  the  same  as  those  for 

thp  formntinn  nf  fl-tl  hy  tjio  Of  prCtOHS 

(p)  on  positronium,  that  is, 

<rn=tr(p+Ps-.II+e-)=a(p+Ps— .II+e  +  )  .  (3) 


Furthermore,  the  cross  sections  for  the  formation  of  H  by 
the  impact  of  protons  on  Ps  are  related,  by  the  principle 
of  detailed  balance,  to  the  cress  sections  for  the  forma¬ 
tion  of  Ps  by  the  impact  of  positrons  on  hydrogen. 
Hence  one  can  write 

,  k< 

Oji  =  tr(p  +  Ps— -H-l-e  +  )=— tr,,s  .  (4) 

kf 

ct,,s  is  the  cross  section  for  the  process 

e  +  +  H(«/m)-*Ps(  ls)+p  ,  (5a) 

or  for  the  process 

c  f +H(ls)— *Ps(n/m)+p  .  (5b) 

Jik,  and  ftYf  are  the  relative  momenta  of  the  positron  (in 
the  initial  channel)  and  the  positronium  (in  the  final  chan¬ 
nel)  of  Eqs.  (5).  In  the  present  work,  the  values  of  <7|>,  for 
the  processes  of  Eqs.  (5)  for  various  sets  of  nlm  arc  calcu¬ 
lated  first,  and  then  Eq.  (4)  is  used  to  obtain  the  values  of 
the  cross  sections  (Op)  for  the  formation  of  ar'.ihydrogen 
which  are  now  labeled  cr„/m. 

The  coordinate  system  used  for  calculating  oPs  of  Eqs. 
(5)  is  shown  in  Fig.  1.  It  is  clear  from  Fig.  1  that  r,  (r^) 
and  R,  (Ry)  arc  the  interna!  coordinate  and  the  coordi¬ 
nate  of  the  center  of  mass  of  the  bound  system  in  the  mi- 


FIG.  1.  The  coordinate  system  for  the  processes  of  Eqs.  (5). 
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tial  (final)  channel,  respectively.  R  locates  the  incoming 
positron  with  respect  to  the  proton.  Atomic  units  arc 
used  throughout  this  paper  unless  specified  otherwise. 
To  obtain  the  impact  energy  of  the  positron  for  the  pro¬ 
cesses  of  Eqs.  (5)  corresponding  to  an  impact  energy  of 
the  antiproton  for  the  processes  of  Eqs.  (1)  and  (2),  con¬ 
servation  of  momentum  is  used.  Let  E  ~k*/2mp  be  the 
laboratory  impact  energy  of  the  antiproton  of  mass  mp 
with  respect  to  the  positronium  at  rest.  Then  the  wave 
number  fcy  of  the  positronium  with  respect  to  the  station¬ 
ary  antiproton  is 

kf=mVtkf/mp  ,  (6) 

where  /»i !>s  is  the  mass  of  the  positronium.  Hence,  in  the 
frame  of  a  stationary  antiproton,  the  energy  of  the  posi¬ 
tronium  is  kj /{2m |is).  The  wave  number  k  of  the  posi¬ 
tron,  of  Eqs.  (5),  is  obtained  from  the  com*.,  .ation  of  en¬ 
ergy,  namely, 

k}/(2\\)=kf/(2\'f)—Zf+z,  ,  (7) 

where  c,  and  Cy  are  the  (positive)  binding  energies  of  II 
and  Ps,  respectively.  If  m  is  the  mass  of  an  electron  (or  a 
positron),  then  v,  =m (m  +  mp)/(2m  +»ip)  is  the  re¬ 
duced  mass  of  the  system  (e4  +11)  in  the  initial  channel 
and  vl=2mpmA2m+inp)  is  the  reduced  mass  of  the 
system  (Ps+p)  m  the  final  channel  of  Eqs.  (5).  In  the  first 
Born  approximation,  the  differential  cross  section  for  the 
processes  (5)  is  related  to  the  T-matrix  element  as 


(2a-)2  kj 


((0/  |  V,  |  0, )  m  the  prior  form  , 
7  ~  J'/ltfr/)  >»  the  post  form  . 


0,  and  t hf  are  the  wave  functions  for  the  system  (e  f,H) 
in  the  initial  channel  and  the  system  (Ps,p)  in  the  final 
channel,  respectively;  that  is, 

0,  =exp(;k(-R,  )0((r,) 


where  0,( r, )  and  <J/(ry)  are  the  bound-state  wave  func¬ 
tions  of  H  and  Ps,  respectively.  The  initial  and  the  final  , 
interaction  potentials  V,  and  Ky  are 

V,  =  V(R)-V(r;) ,  (1  la) 

V/  =  V(R)-V{rl)  ,  (lib) 

where  V(x)=l/x.  The  integrated  cross  section  is  given 
by 


sinOrfO  . 


0y=eXp(lky'Ry)0y(ry) 


Substituting  V,  and  Fy  from  Eqs.  (11)  in  the  T-matrix 
elements,  we  get 

<0y|  F(R)|0i>-<0y|  F(r/)|0,> 

=  A', +  A’,  (prioi  form)  ,  I* 3a) 
T=  <0/1  F(/?)|0,>-<0/|  F(r,)|0,> 

=  Ar,+iVj  (post  form)  .  (13b! 

These  '/’-matrix  terms,  Ar1,  A'2,  and  A'_,,  can  be  evaluated' 
easily  for  the  cases  when  only  the  spherically  symmetric 
states  are  involved.  Such  cases,  for  example,  arc  the  for¬ 
mation  of  Ps(l.v)  from  different  s  states  of  11  of  Eq.  (5a)  or 
the  formation  of  different  s  states  of  Ps  from  III  lx)  of  Eq. 
(5b).  Complexity  arises  when  the  initial  or  the  final 
bound-state  wave  functions  in  Eqs.  (5)  are  not  spherically 
symmetric  as  in  the  case  of  formation  of  Ps  in  the  p  or  d 
state  or  when  the  target  II  is  in  an  excited  p  state.  The 
complexity  due  to  the  nonspherical  wave  functions  can  be 
reduced  by  expressing  the  angular  dependence  of  the 
wave  functions  in  terms  of  derivatives  of  an  exponential 
factor  exp(i  A-r),  where  the  value  of  the  parameter  A  is 
eventually  set  to  zero.  In  particular,  for  the  processes 
(5a)  we  introduce  a  function  ,Y,(r)  as  follows: 

,V,(«/m,r)=p(«/m,r)exp(;  A-r)  ,  (14) 

where  p(nlm,r )  represents  R„t,  the  radial  part  of  the 
wave  function  of  hydrogen,  multiplied  by  the  constants  of 
Ylm,  the  angular  part  of  the  wave  function.  Then  the 
complete  bound-state  wave  functions  0,(r),  of  Eq.  (10), 
can  be  obtained  by  merely  taking  the  appropriate  deriva¬ 
tives  of  Xj(t)  with  respect  to  either  Ax  or  Ay  or  Az  and 
then  setting  A=C.  As  an  example,  for  the  3 d0  state  of 
hydrogen,  the  wave  function  is 


0/(3r/o,r)=R,2(r)  (3 cos2/?—  I ) 

1077 


=p(3d0,r)  —  -- 
r  3  ; 


r|exp(i  AT)o(3</?.r)/r2]— exp(/  AT)n(7r/c  r)/3 


-[.Vi(3(/0,r)/r2]-,V,(3d0,r)/3 
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•  where,  using  the  explicit  form  of  R}2, 

p(3d0,r)=  ^^^■/•iexp(-r/3)  .  (16) 

This  kind  of  trick  has  been  used  earlier  by  Sil  cl  al .6  in 
calculations  of  charge-transfer  cross  sections. 

Let  Q,-  and  Qy  be  the  momentum-transfer  vectors  in 
the  initial  and  final  channel,  respectively,  t.e., 

Q,=[»ip/(/m +m(>>]kl— ky,  Q/=k//2-k,  (17) 

and  let  p  be  the  reduced  mass  of  the  final  atom  positroni- 
um  of  the  processes  (5).  Then  the  parts  of  the  7-matrix 
elements  for  the  processes  (5a),  upon  replacing  <4,(r)  by 
,V,(«/m,r)  of  Eq.  (14),  can  be  written  as 

JV(  =(2tt)>/2 f  ^/(Qy  — t)D(t)p(  — t-Q,  —  AW3t  ,  (18a) 

q2 

N2  —  — (2ff)5  ■=-  +£/  0/(Q/)p(-Q,  —  A)  ,  (18b) 

Zfl 

Ni  =  -{2ir)3/2$}{Q/) f  P( t )p(  —  t  —  Q,  —  AWJl  ,  (18c) 


where  the  ovcrtildc  represents  the  Fourier  transform  of 
the  respective  function.  For  example,  the  Fourier  trans¬ 
form  /(t)  of  a  function  /(r)  is  defined  as 

/(t)=(2 it)' 3/2  f  rfVexpl— it-r)/(r) . 

Substituting  the  explicit  forms  for  the  potential  V,  and 
the  wave  functions  <4,  and  4>f,  the  term  incurs  an  in¬ 
tegral  of  the  form 

/,=  f  d3i - - - r - - ,  (19) 

J  (t-P )2+p2  <t-D)2+b2  t2+z 2 

which,  upon  using  the  Feynman  identity,  can  be  reduced 
to  a  one-dimensional  integral  as5 

/ 1  =  jr2  f 1 - “■ - j-  ,  (20) 

1  Jo  E[F2+(E+z)2] 

where  E2=x{  1  —  x)(P— B)2-fxp2+(  1  —  x)b2  and 
F=(P  —  B)x  -f-B.  N2  can  be  written  in  a  pure  analytic 
fashion.  The  term  N3  involves  integral  of  the  foim 


/3=JV 


‘[(t_p)l+p’]2*2+*2 
which  can  be  done  analytically  to  obtain 


p(P2-l-(p+z)2]  ‘ 


Integral  of  Eq.  (21)  easily  follows  from  integral  /,  of 
Eq.  (19)  by  setting  P  =  B  and  p  =b.  P,  B,  p,  b,  and  z  of 
Eqs.  (19)  and  (21)  are,  of  course,  related  to  Q,,Q/  and  the 
parameter  A  of  Eqs.  (18a)  and  (18c).  In  the  case  of  p- 
and  r/-s!ate  wave  functions,  A'),  N2,  and  N}  will  involve 
derivatives  with  respect  to  the  Cartesian  components  of 
A.  Note  that  since  <4„/m  is  a  complex  conjugate  of 
<5„;  and  |  7'  j 2  is  used  to  calculate  the  cross  section, 
the  cross  sections  with  both  the  wave  functions  (<4„; ±m) 
are  the  sam-.  This  technique  of  using  the  derivatives  of 
the  factor  exp(i  A^r)  for  representing  the  angular  depen¬ 
dence  of  the  wave  functions  can  also  be  used  for  the  for¬ 


mation  of  various  excited  states  of  Ps  from  H(ls)  of  Eq. 
(5b).  In  this  case  the  factor  exp(t  A-r)  is  introduced  in 
the  final  bound-state  wave  function  r)  rather  than  in 
the  initial  bound-state  wave  function  <4,(r),  Thus  after 
calculating  the  cross  sections  for  the  processes  of  Eqs.  (5), 
the  corresponding  cross  sections  for  the  formation  of  an¬ 
tihydrogen  arc  obtained  using  Eq.  (4). 

It  was  predicted  by  Oppenheimer7  that  the  cross  sec¬ 
tion  for  capture  into  any  ns  excited  hydrogenic  state  at 
high  incident  energies  falls  as  1/n3.  Later,  Omidvar8 
showed  that  the  cross  section  for  the  capture  into  s,  p, 
and  d  states  as  well  as  the  sum  of  the  cross  sections  over 
hn  states  falls  as  l//i3  in  the  first  Born  approximation. 
There  is  also  an  experimental  evidence  of  the  cross  sec¬ 
tion  for  electron  capture  into  excited  states  falling  as 
l/«  3  in  collisions  of  high-energy  fluorine  ions  with  ar¬ 
gon.9  Hence  using  the  l/«3  dependence  of  the  capture 
cross  sections  of  Eq.  (1)  for  the  formation  of  II  in  excited 
states  («  >4),  the  total  cross  section  could  be  written  as 


where 


^i  =  0,h0  > 

<7?=<72!o+"2?0  +  2^2/.|  -  <24> 

<r}  =  ‘7>!0  +  a}p0  +  ‘ia)?l+a!tl„  +  7-a3dl-'’  2 <7m2  . 

and  n  =3.  The  value  of  the  seiies  within  the  large  square 
brackets  of  Eq.  (23)  is  2.0805  for  n  =3.  In  the  case  of 
H(ls)  formation  from  various  states  of  Ps,  of  Eq.  (2),  the 
1  /n3  rule  is  used  for  positronium  states  n  >  3.  Hence  us¬ 
ing  n  =2  the  value  of  the  series  within  the  large  square 
brackets  of  Eq.  (23)  is  1.616. 

III.  RESULTS  AND  DISCUSSIONS 

The  present  work  on  charge-transfer  processes  has 
been  carried  out  using  both  the  prior  and  the  post  forms 
of  the  interaction.  The  difference  in  the  values  of  the 
cross  sections  using  these  two  forms  of  the  interaction 
(that  is,  the  post-prior  discrepancy)  is  quite  small  (a  few 
parts  in  a  thousand).  Hence  only  the  prior  form  of  the  re¬ 
sults  is  presented.  The  n  ~ 3  behavior  for  the  capture 
cross  sections  is  used  for  antihydrogen  formation  in  states 
n  >  4  from  Ps(ls)  [sec  Eq.  (1)]  and  for  H(ls)  formation 
from  Ps  in  stales  n  >  3  [see  Eq.  (2)].  The  differential  cross 
sections  (DCS’s)  for  the  formation  of  antihydrogen  in  all 
stales  tgiound  and  ah  possible  excited  states  added  to¬ 
gether)  by  the  impact  of  antiprotons  on  positronium  in 
ground  state  are  shown  in  Fig.  2  and  those  for  the  forma¬ 
tion  of  antiliydrogen  in  the  ground  state  from  all  possible 
slates  of  positronium  are  shown  in  Fig.  3.  From  Figs.  2 
and  3  we  see  that  the  nature  of  flic  DCS  curves  in  these 
two  cases  are  very  similar.  DCS  curves  in  Fig.  2  show 
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that  following  the  forward  peak  there  is  a  very  shallow 
minimum  at  lower  energies  that  deepens  with  increasing 
antiproton  impact  energies  while  its  position  stays  at 
roughly  the  same  angle  (~24°).  At  larger  scattering  an¬ 
gles,  the  values  of  DOS's  fall  smoothly  at  all  energies. 
Similar  is  the  case  with  the  curves  of  Fig.  3,  except  that 
the  minima  in  the  DCS  values  at  lower  impact  energies 
are  more  pronounced  than  those  of  Fig.  2.  The  minima 
in  the  DCS  curves  arise  from  the  opposite  nature  (attrac¬ 
tive  versus  repulsive)  of  antiproton-positron  and 
antiproton-elcctron  interactions;  essentially  the  N ,  and 
the  iV2  parts  of  the  T-matrix  element  [Eq.  (13)]  cancel  one 
another.  The  integrated  cross  sections  for  the  formation 
of  antihydrogen  in  individual  states  ls0,  3s0,  2 p0,  2 pi( 
3s0,  3Po>  3pi,  3 d0,  3rf|,  id 2,  as  well  as  in  all  possible 
states,  using  the  n  '1  behavior,  from  Psds),  arc  numeri¬ 
cally  presented  in  Table  I  and  arc  shown  in  Fig,  4.  The 
integrated  cross  sections  for  the  formation  of  H(  Is)  from 


FIG  2.  The  differential  cross  sections  for  (lie  formation  of 
amiliydrogcn  in  all  possible  slates  by  the  impact  of  anliprolons 
on  I’s(lj),  Eq,  (1).  The  DCS  curves  between  200  and  500  heV 
correspond  to  sequential  increase  of  impact  energy  at  an  inter¬ 
val  of  25  keV. 


Ps  in  individual  states  ls0,  2s0,  2 p0,  2 p,,  as  well  as  in  all 
possible  states,  using  the  n~3  behavior,  arc  numerically 
presented  in  Table  II  and  are  shown  in  Fig.  5.  All  the 
curves  for  integrated  cross  sections  in  Figs.  4  and  5  decay 
very  smoothly  with  increasing  p  impact  energies.  From 
Table  I  and  Fig.  4  we  see  that  formation  of  the  p  and  d 
states  of  antihydrogen  dominates  over  the  formation  of 
any  other  state  at  lower  antiproton  impact  energy  ( <  50 
keV).  However,  at  higher  impact  energies  the  formation 
of  H(  Is)  dominates  overjlie  formation  of  all  other  states. 
In  case  of  formation  of  H(  Is)  from  various  states  of  Ps,  as 
seen  in  Table  II  and  Fig.  5,  the  formation  from  Ps(ls) 
dominates  significantly  over  formation  from  any  other 
slate  of  positronium  almost  at  all  energies  considered. 

The  use  of  1  /n 3  behavior  of  the  capture  cross  sections 
for  the  entire  range  of  antiproton  impact  energies  to  cal¬ 
culate  the  cross  sections  for  states  n  ;>4  for  processes  of 
Eq.  (1)  and  for  states  n>3  for  piocesses  of  Eq.  (2)  is 


FIG.  3.  The  differential  cross  sections  for  the  formation  of 
antihydrogen  in  the  ground  stale  by  the  impact  of  anliprolons 
on  all  possible  stales  of  positronium,  Eq.  (2)  The  DCS  curves 
between  200  and  500  keV  correspond  to  sequential  increase  of 
impact  energy  at  an  interval  of  25  keV. 
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found  to  be  quite  reasonable.  For  example,  the  value  of 
the  ratio  (;i  +\)}a „ +  l/nia which  ideally  should  be  1  if 
the  l/;i3  behavior  is  valid,  is  1.54  at  20  kcV  and  1.07  at 
500  kcV  for  the  process  of  Eq.  (1)  when  n  =2.  The  same 
ratio  for  the  process  of  Eq.  (2)  is  1.26  at  20  kcV  and  1.C2 
at  500  keV  for  n  =  1.  These  numbers  suggest  that  the  use 
of  1/n3  rule  holds  relatively  more  accurately  at  higher 
impact  energies  than  at  lower  impact  energies. 

In  the  earlier  works,  Neuman  el  a/.1  have  carried  out  a 
laser-enhanced  electron-ion-capture  calculation  to  evalu¬ 
ate  the  rate  of  total  recombination  for  antihydrogen  for¬ 
mation.  Humberston  et  al.2  presented  the  values  of  cross 
sections  for  the  formation  of  H(ls)  from  l’s(lx)  at  lower 
impact  energies  using  the  p-,  and  d- wave  positronium- 


formation  phase  shifts  obtained  by  elaborate  variational 
calculations.  Later,  Darcwych3  calculated  the  cross  sec¬ 
tions  for  the  formation  of  II  in  states  1*0,  2*0,  2p  (all  ml, 
and  3s0  from  1’sUs)  at  low  antiproton  impact  energies  us¬ 
ing  the  first  Born  approximation  (FDA).  Since  the  FBA 
is  not  a  good  approximation  for  scattering  calculations  at 
low  projectile  energies,  the  present  work  is  carried  out  in 
an  intermediate  range  of  impact  energies.  However,  for 
computational  checking  purposes,  cross  sections  at  a  few 
low  impact  energies  were  calculated  and  it  was  observed 
that  both  Ref.  3  and  the  present  work  agree  veiy  well  for 
the  formation  of  II  in  *  states,  while  the  values  of  the 
cross  sections  for  2 p-  (all  m)  state  formation  in  the 
present  work  arc  much  larger  than  those  of  Ref.  3.  As  a 


TABLE  I  The  values  of  the  integrated  cross  sections  (un)  for  the  formation  of  11  in  ground  and  various  excited  slates  as  well  as  in 
all  possible  states  by  the  impact  of  p  on  l’s(ls)  in  the  incident-energy  range  of  £  =20-500  keV.  The  notation  a[b]  for  the  values  of 
the  cross  sections  means  a  X  10*. 


Op  (10" 16  cm2)  at  various  p  energies,  E  (keV) 


States 

£  =  20 

30 

40 

50 

75 

100 

125 

150 

l*o 

0.1048(1] 

0.4640 

0.2302 

0.1242 

0  3425(  -  1) 

0.12081-1] 

0  5025[-2) 

0,23521-2) 

0.2780 

0.1109 

0.53721-  1] 

0.28371-1) 

0.7344J  -  2) 

0.24381-2) 

0.9624[~3) 

043151-3] 

2/>o 

0.9554 

0  2794 

0.9847(-l] 

0  39941-1) 

0.63841-21 

0.15161-2) 

0.4643[-3] 

0.l697[-3] 

2pi 

0.4057 

0.1001 

0.31501-1] 

0.11761-1] 

0.1643[-2] 

0.35931-3) 

0.10401-3] 

0.3646[-4] 

3*o 

0.1149 

0.4015(-1] 

0.18621-1] 

0.96481-2) 

0.2439[-2] 

0  79731-3) 

0.31121-3) 

0.1383[-3] 

3p0 

0.3173 

0.1009 

0.3674(-l] 

0.15111-1) 

0.24291-2) 

0.5741[-3] 

0.17471-3) 

0.6345[-4) 

3p, 

0.1264 

0.3439[-l] 

0.I124(-I] 

0.42701-2] 

0.60421-3) 

0.13211-3] 

0.38131-4] 

0.13331-4) 

3d„ 

0.933 1  [  —  1  ] 

0.23051-1] 

0.66511-2] 

0,2241(-2] 

0.24511-3) 

0.43591-4] 

0.10611-4] 

0.32111-5] 

3d, 

0.5982[-l] 

0.1323[-1] 

0.35641-2] 

0.11441-2) 

0  11571-3) 

0  196l[-4] 

0.4613[-5] 

0.1360J-5) 

3d, 

0.1886[ -1] 

0.3340[-21 

0  79961-3] 

0.2379[-3] 

0.21491  —4] 

0  34l7[-5] 

O7709|-6] 

0.22061-6) 

All 

0.5039(1] 

0.1608(1] 

0.6393 

0.2957 

0.64991  -  1  ] 

0.2035[  - 1 ) 

0.7873(-  2) 

0  35141-2] 

£■  =  175 

200 

225 

250 

275 

300 

325 

350 

1*0 

0.1204[-2] 

0.66111-3] 

0.38431-3} 

0.2340(-3) 

0.14831-3) 

0.9710|-4] 

0  65461-4] 

0  45251-4] 

2*o 

0.2I31[  -3] 

0.11361-3] 

0.64381-4] 

0.3837(  -4) 

0.2385[-4) 

0.1537(-4] 

0.10211-4] 

0.69691-5] 

2po 

0.7069[  —4] 

0  3256[-4) 

O.I624{-4] 

0  8637(-5] 

0.48481-5] 

0.2847(-5) 

0.1737[-5] 

0.10961-5] 

2Pi 

0. 147 1  [  —  4  j 

0.66041-5) 

0.32251-5) 

O.I686[-5) 

0.93241-6) 

0.5405[-6] 

0.3262(-6] 

0  2038[-6] 

3s0 

0.6786[  —  4) 

0.35981-4] 

0.20301-4] 

0.12061-4] 

0.74781-5] 

0.48081  -5] 

0.31891-5] 

0.21731-5] 

3p0 

0  2629[—4] 

0.12051-4] 

0.5984[-5] 

0.31721-5] 

0  17751-5] 

0.10391-5) 

0.63251  -6) 

0.398l(-6] 

3Pi 

0.53601-5] 

0.2399[-5) 

0.11691-5] 

0.60961  —  6] 

0.3364(-6) 

0.1946[-6) 

0.11731-6] 

0.7315|~7] 

3d0 

0.1l41(-5] 

0.45831-6] 

0.20301-6] 

0.9738[-7] 

0.4990[ -7  J 

0.27041-7) 

0.1537[-7] 

0.91121-8] 

3d, 

0.47291-6] 

0.18641—6] 

0.81111-7] 

0.38221-7] 

O.I923[  —  7) 

0.1023[-7) 

0.57021-8) 

0.3309[-8] 

3‘  2 

0.7495[-7] 

0.290l(-7] 

0.12441  -7] 

0.5791[-8) 

0.2S84[  — 8] 

O.I520[-8] 

0  8408(-9) 

0  48461-9] 

All 

0.1740[-2) 

0.93231-3) 

0.53171-  3] 

0.31901-3] 

0.1997[-3) 

0.1295[-3] 

0.8656[-4] 

0.594l(-4) 

£=375 

400 

425 

450 

475 

500 

1*0 

0.3l98[-4] 

0.2304[-4] 

0.16901-4] 

0.1259[-4] 

095121-5] 

0.72801-5] 

2*o 

0.4869[-5] 

0.3473[-5] 

0.252d[  —  5] 

0.18641-5} 

0.13981-5) 

0.1062(-5] 

2po 

0.7l21(-6] 

0.47461-6) 

0.32361-6) 

0.225I[-6) 

0.1595[-6) 

0.11491-6) 

2pi 

0.13l2[-6) 

0.86771-7] 

0.58741-7] 

0.40601-7] 

0.28601-7) 

0.2049[-7) 

3*o 

O.I517[-5] 

0.10811-5) 

0.78441-6] 

0.5790[-6) 

0.43381-6] 

0.3294(-6) 

■>Po 

0.25811-6] 

0.17161-6] 

0.11681-6) 

0,81121  -7) 

0.57381 -7j 

0.4l27(-7] 

3 Pi 

0.4703[-7] 

0.3106(-7] 

0.21001-7] 

0  1450[ -7J 

O.I020[-7] 

0.73011-8) 

3  d0 

0.5604(-8] 

0.3560[-8] 

0.2328(-8] 

0.15621-8] 

0.10741-8) 

0.7534[-9] 

3d, 

O.I9891-8] 

0.1233[-8] 

0.78601-9] 

0.51331-9] 

0.3426[-9) 

0.2332[-9] 

3d2 

0.28961-9] 

0.1786(-9] 

0.11321-9] 

0.7361[-10] 

0.48931-10) 

0.3319(-10] 

All 

0.4173(-4] 

0.29911-4] 

0.2184[-4] 

0.1620[-4] 

0.12191-4) 

0.9303[-  5) 
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Antiproton  Impact  Energy  (keVI 

FIG.  4  The  integrated  cross  sections  for  the  formation  of  an- 
tihydrogcn  in  various  s,  p,  and  d  states  (summed  over  m  quan¬ 
tum  numbers)  by  the  impact  of  antiprotons  on  1’sOsl,  Eq.  tl). 
The  solid  curve  corresponds  to  the  total  integrated  u.c.,  sum  of 
all  possible  nl  states  of  aiitihydrogcn)  cross  sections. 


FIG.  5.  The  integrated  cross  sections  for  the  formation  of  an- 
lihydrogen  in  the  ground  state  by  the  impact  of  antiprotons  on 
po  tronium  in  various  s  and  p  states  (summed  over  m  quantum 
numbers),  Eq  (2).  The  solid  curve  corresponds  to  the  total  in¬ 
tegrated  (i.c.,  sum  of  all  possible  nl  states  of  positromum)  cross 
sections. 


TABLE  II.  The  values  of  the  integrated  cross  sections  (<7n)  for  the  formation  of  0(  Is)  by  the  impact 
of  p  on  ground  and  various  excited  states  as  well  as  on  all  possible  states  of  Ps  in  the  incident-energy 
range  of  E  =20-500  keV.  The  notation  a[b]  for  the  values  of  the  cross  sections  means  a  X  10*'. 

E  (keV) 

2s0 

Ps  state 
<7|j  (10" 16  cm2) 

2 Pi 

All 

20 

0.1048(1] 

0.1054 

0.39341-1) 

0.1020[-l] 

0.1315(1] 

30 

0.4640 

0.5235[-l) 

0.16591-1] 

0.4220[-2] 

0.5890 

40 

0.2302 

0.27571-  1] 

0.7753[-2] 

0.18771-2] 

0.2933 

50 

0.1242 

0.1539[-1] 

0.3896[-2) 

0.89781-3] 

0.1583 

75 

0.3425[  —  1  ] 

0.44081-2) 

0.88891-3) 

0.18541-3] 

0.4341(-1] 

100 

0  1208(-1  j 

0.1572[-2] 

0.26221-3] 

0.5093(  —  4] 

0.1521[-1) 

125 

0.5025(-2] 

0.6545[-3] 

0.92911-4) 

0.17U[-4] 

0.6288[-2) 

150 

0.2352[-2] 

0.3054(-3) 

0.37731-4) 

0.6665[-5] 

0  29281-2) 

175 

0.1204(-2] 

0  15551-3] 

0.1702[-4] 

0.2907(  —  5] 

0  14921—2] 

200 

0.66111-3) 

0.8490(-4) 

0.8342[-5] 

0.l387[-5] 

0.81631-3) 

225 

0.3843[-3] 

0.4903(-4) 

0.43751  -5) 

0.7106[-6] 

0  47291-3) 

250 

0.2340[-3] 

0.2967(-4] 

0  2427[-5] 

0.38641-6) 

0.28721-3) 

275 

0.1483[-3] 

0.1867[-4] 

0.14111-5) 

0.22081-6] 

0  1814[-3] 

300 

0.97101-4] 

0.12161-4) 

0  8544[-6) 

0.13161  -6] 

0.11861-3] 

325 

0  6546[-4) 

0  8146(-5] 

0.5355[~6) 

0.81321-7) 

0  79751  -4) 

^0 

04S?S( -4] 

0  55991-51 

0,3459[  —61 

0.51871-71 

0.55031-41 

375 

0  319si  —4) 

0.39361-5) 

0.22951 -6j 

0.3401(-7] 

0.38821  -4) 

400 

0.23041-4] 

0.28231-5] 

0.15581  -6) 

0.2285|-7) 

0.27931  -4) 

425 

0.16901-4] 

0.20601-5) 

0.10811-6) 

01569(-7] 

0.20461-4) 

450 

0.1259(-4] 

0.15281-5] 

0.7639[-7) 

0.10991  —7} 

0.I522[-4) 

475 

095121-5) 

0.11501  -5] 

0.54921-7] 

0.78291-8] 

0.1149[  -4] 

500 

0.72801-5] 

0.8768[-6] 

0  40101-7) 

0  5668(-8] 

0.8781(-5] 
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check  on  the  present  computer  code,  the  cross  sections10 
for  the  formation  of  H,  by  the  process  of  electron  capture 
during  proton-hydrogen  collisions,  in  states  Is  and  2 p 
were  reproduced. 

In  conclusion,  we  have  presented  a  simple  first  Born 
calculation  of  cross  sections  for  the  formation  of  antihy- 
drogen,  by  the  impact  of  intermediate-energy  (20-500 


ft.  Berger  cl  at.,  CERN  Report  No.  I’SCC  1985,  p.  86  (unpub¬ 
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keV)  antiprotons  on  positronium,  in  ground  and  in  vari¬ 
ous  excited  states. 
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Approximate  Green's  Functions  in  Electron 
(Positron)»  Atom  Collisions 


The  role  of  approximate  Green's  functions  in  electron  (positron)-atom  collisions 
is  discussed.  It  is  explicitly  demonstrated  that  the  scattering  amplitude  in  various 
approximation  methods,  like  the  Massey- Mohr,  simplified  second  Born,  plane 
wave  (adiabatic  as  well  as  nonadiabatic).  Glauber  (with  and  without  the  Wallace 
correction),  eikonal  Born  senes  (including  unitarization)  and  Schwinger  variational 
pnnciple.  is  obtained  by  mere  approximations  of  the  relevant  Green's  function. 
Some  of  the  shortcomings  of  different  methods  are  pointed  out.  It  is  shown  that 
the  Schwinger  variational  method  with  suitable  choice  of  tnal  wave  function  yields 
practically  the  same  values  of  the  differential  and  integrated  cross  sections  for  the 
elastic  scattenng  as  well  as  the  total  collisional  scattering  for  the  e--H  system  as 
obtained  by  the  use  of  the  unitarized  eikonal  Born  series  method. 

Key  Words:  electron-atom  scattering,  positron-atom  scattering,  Green's  functions, 
approximate  methods 


It  is  well  known  that  an  exact  evaluation  of  the  electron  (positron- 
atom  scattering  amplitude  is  not  possible  at  present.  A  number  of 
approximate  methods  starting  from  the  Lippmann-Schwinger  in¬ 
tegral  equation  have  been  developed  for  an  evaluation  of  the  direct 
scattering  amplitude.1  In  this  paper  we  will  explicitly  demonstrate 
that  most  of  these  approximate  methods  for  evaluating  the  scat¬ 
tering  amplitude  employ  different  approximate  forms  of  the  Green’s 
function.  We  will  also  compare  somewhat  in  detail  the  numerical 
results  of  the  elastic  and  total  scattering  of  electrons  and  positrons 
by  a  hydrogen  atom  obtained  recently  by  Byron  et  al.2  using  the 
unitarized  eikonal  Born  series  (UEBS)  method  and  by  Khare  and 
co-workers3-4  using  the  Schwinger  variational  (SV)  principle. 
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The  exact  direct  scattering  amplitude  for  the  e  =  -atom  scattering 
from  the  initial  target  state  i  to  the  final  state  /  is  given  by1 

/"(k/.k,)  =  -2~  <khf\viiG;Ur-%u)mi?B'  (1) 

n  "  1  n  ■  l 


with 


Here  k,  and  k;  are  initial  and  final  momenta  of  the  projectile  and 
U  is  the  reduced  interaction  energy,  that  is,  U  =  (2 mlhz)  V  where 
V  is  the  actual  total  Coulomb  interaction  between  the  projectile 
(electron  or  positron)  and  the  target  atom.  GJ  is  the  free  Green’s 
function  of  the  noninteracting  projectile  and  target.  The  Green’s 
function  GJ  is  given  by 


c<r 


|q,m)  (q,m| 

«,  -  r  +  *  q 


(2) 


where  |q)  and  | m)  are  the  intermediate  states  of  the  projectile  and 
target,  respectively,  and  A: *  =  k*  -  2 (Em  -  £,),  Em  and  £,  being 
the  intermediate  and  initial  target  state  energies.  (We  will  use 
atomic  units,  that  is,  h  =  l,m  -  lande  =  1  unless  stated  explicitly 
otherwise).  The  nth  term  of  series  (1),  /Bn,  is  referred  to  as  the 
nth  Born  term  and  fBp  is  the  pth  Born  approximation  for  the 
scattering  amplitude.  The  nonconvergence  of  the  Born  series  for 
rearrangement  collisions,  including  the  exchange  effects,  was  ex¬ 
plicitly  demonstrated  by  Aaron  et  al.5  We  may  remark  that  the 
investigations  of  Rosenberg6  strongly  suggest  that  the  Born  series 
for  the  direct  scattering  amplitude  for  the  e=-atom  system  may 
also  not  converge.  Furthermore,  the  evaluation  of  fBl/  for  n  >  1, 
is  extremely  difficult  due  to  the  presence  of  the  operator  (G0  + 
U)n~l.  The  first  Born  term,  fBU  which  does  not  contain  Go,  is 
easy  to  evaluate.  However,  it  completely  neglects  the  effects  due 
to  distortions  of  the  wave  functions  of  the  projectile  and  the  target. 
The  Green’s  function  contains  these  effects  through  the  interme¬ 


diate  states  jq)  and  j m).  Different  approximations,  which  have  been 
developed  to  obtain  ff,{ k/,k,)  using  an  integral  equation  approach, 
treat  the  Green’s  function  in  different  approximate  ways. 


For  fast  projectiles  the  most  important  correction  to  fB\  is  the 
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second  Born  term  fB2  given  by 

Jb2  (k/.k,)  =  -2n-  (k,.f\UG0'U[kJ).  (3) 

An  exact  evaluation,  either  analytical  or  numerical,  of  fB2  is  an 
involved  problem.  So  far  only  Ermoiaeve  and  Walters7  have  suc¬ 
ceeded  in  this  task.  Their  numerical  method  has  been  discussed 
briefly  by  Walters.8 

The  first  attempt  for  an  approximate  analytical  evaluation  of  JB2 
was  made  by  Massey  and  Mohr.9  They  replaced  k*  by  k ]  in  (2) 
for  all  values  of  m.  Then  using  the  closure  relation  one  obtains 
the  Green’s  function  Go(R.R':  x.x')  =  (R,x|Go|R',x'>  in  the 
Massey- Mohr  approximation  given  by 


Go+w(R.R';  x.x')  =  Gop(R.R')  8(x  -  x').  (4) 

where  R  and  x  are  the  coordinates  of  the  projectile  (electron  or 
positron)  and  the  bound  state  electrons  of  the  target,  respectively. 
G0p  is  the  free  particle  Green’s  function  of  the  projectile  only  and 
is  given,  in  operator  form,  by 


y*  -  M 

0p  J  k}  - 


|q>  <ql 


ql  +  it 


dq 


(5) 


Use  of  (4)  and  (5)  in  (3)  yields  an  integral  expression  for/^  which 
does  not  contain  intermediate  target  states  | m)  and  the  integration 
over  x'  becomes  trivial  due  to  the  presence  of  the  delta  function 
5(x  -  x').  Thus,  in  the  Massey-Mohr  approximation  we  have 


where  Jg^x)  is  the  second  Born  scattering  term  for  the  scattering 
of  the  projectile  by  a  target,  in  which  the  bound  electrons  are 
frozen  at  the  position  x,  and  is  given  by 


?«(x)  =  -  2tt2  (k^Go^Ulk,). 


(7) 


It  is  evident  from  (6)  and  (7)  that  the  Massey-Mohr  approxi¬ 
mation  completely  neglects  the  effects  due  to  the  distortion  of  the 
target  wave  function.  Since  this  approximation  assumes  the  exci¬ 
tation  energies  to  be  zero,  it  turns  out*  that  for  the  elastic  scattering 
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the  imaginary  part  of  the  second  Born  term.  Im  diverges  in 
the  forward  direction  and.  for  e--atom  elastic  scattering,  the  real 
part  in  the  forward  direction.  Re  7S2.W(0),  asymptotically  goes  as 
k,~-  instead  of  k,~l.  the  correct  dependence.10  A  similar  divergence 
in  the  forward  direction  is  also  observed  for  s-s  excitation,  for  the 
e  =  -atom  system,  which  arises  due  to  the  absence  of  the  interme¬ 
diate  p  states  in  the  evaluation  of  the  second  Bom  term.  To  a 
certain  extent  these  discrepancies  are  removed  by  taking  /c£  = 
k f  -  2A  with  A  as  the  mean  excitation  energy.  In  an  average  sense 
the  signature  of  the  intermediate  target  states  are  carried  to  JB2 
through  A.  The  choice  of  A  is  not  unique.  Ermolaeve  and  Walters7-8 
have  discussed  the  various  options.  To  obtain  better  values  of  the 
differential  cross  sections,  particularly  for  the  case  with  momentum 
transfer  |K|  large  in  which  the  higher  Born  terms  play  an  important 
role,  a  few  low  lying  excited  states  should  be  explicitly  included 
in  the  summation  of  Eq.  (2).  Based  on  these  ideas  Holt  and 
Moiseiwitsch10  proposed  to  include  in  Eq.  (2)  a  finite  number  of 
low  lying  target  states  exactly  and  the  rest  by  summing  up  through 
closure.  The  simplified  second  Born  approximation  (SSBA)  of 
Holt  and  Moiseiwitsch,  thus  obtained,  has  been  one  of  the  most 
popular  ways  for  evaluating  the  second  Born  term. 

An  approximation  just  opposite  to  that  of  Massey  and  Mohr 
was  attempted  by  Khare  and  Shobha11  when  they  evaluated/^  in 
the  plane  wave  approximation.  They  took  qz  =  k]  in  (2).  Thus 
Go(R,R';  x.x')  reduced  to  Go,(x,x')'8(R  -  R'),  where  Go,, 
the  target  Green’s  function,  is  given,  in  operator  form,  by 


G£  =  S' 


lm)  (m\ 
k2m  ~  kf 


(8) 


The  prime  excludes  the  final  and  initial  states  from  the  summation; 
hence,  these  states  do  not  contribute  to  the  polarization  of  the 
target  wave  functions.  Use  of  (8)  in  (3)  yields  as  the  second  Born 
term  in  the  plane  wave  approximation. 


WMO  =  -2tt2  (k^|£/ap|k,) 

where 

j,  Qi  <f\U\m)  < m\U\i ) 

3P  5  kzm  -  kj  ■ 


(9) 

(10) 
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It  is  evident  that  (9)  completely  neglects  the  effects  due  to  the  dis¬ 
tortion  of  the  projectile  wave  function  but  includes  those  due  to  the 
distortion  of  the  target  wave  function.  Since  the  projectile  is  being 
represented  by  a  plane  wave  in  (9).  this  approximation  is  referred 
to  as  the  plane  wave  approximation.  Uip  serves  as  the  interaction 
potential  between  the  projectile  and  the  distorted  (polarized)  tar¬ 
get.  Note  that  since  (10)  does  not  depend  upon  the  velocity  of  the 
projectile  the  polarization  potential  Uip  is  adiabatic.  The  approx¬ 
imate  Green's  function,  including  nonadiabatic  effects  up  to  the 
first  order  in  the  interaction,  can  be  obtained  by  taking  qz  =  -  v£ 
and  using 


1  1  1  +  k2 

&  +  ?*  ki-kz  +  {V-R  +  kz)-ki-k:  (ki-k-y 


(in 


in  (2).  We  note  that  (V^  +  k2)  when  operating  on  a  plane  wave 
yields  zero.  On  using  this  approximate  Green’s  function  in  Eq.  (3) 
we  obtain 


fxpfr/M  =  -2tt2  (kf\Ut p  +  UM p|k()  (12) 

where  the  nonadiabatic  potential  Vnip  in  its  Hermitian  form  is  given 
by12 


tf„.p(R)  =  S 


,  VH  (f\U\m)  •  VR  ( m\U\i ) 

(ki  ~  k2)2 


(13) 


It  is  easy  to  show12  that  for  elastic  scattering  the  asymptotic  form 
of  the  dynamic  interaction  potential  is  given  by 


a ,  (a„  -  6(3,) 

tyR)  +  £/fl,p(R)  =  -  ^  -  ■"  ~R(r  •  O4) 

where  ad  and  a,  aie  dipole  and  quadrupole  polarizabilities  and  (3, 
is  the  dipole  nonadiabatic  coefficient  of  the  target,  a  result  obtained 
earlier  by  Kleinman  et  al ,13  using  a  different  approach.  Although 
the  plane  wave  approximation  has  been  successful  in  explaining 
the  differential  cross  sections  for  the  elastic  scattering  of  fast  elec¬ 
trons  by  light  targets  like  hydrogen12  and  helium14  atoms  as  we!! 
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as  hydrogen  molecules. 15  it  suffers  from  the  defect  that  fB2p  is 
purely  real.  Thus,  although  this  approximation  can  be  used  to 
obtain  the  total  elastic  scattering  cross  section  by  integrating 
(/si  +  ?B2P)2  over  all  s°lid  angles,  it  cannot  be  used  to  obtain  the 
total  collisional  cross  section  (elastic  plus  inelastic)  which  is  related 
to  the  imaginary  part  of  the  forward  scattering  amplitude  via  the 
optical  theorem.  Furthermore,  even  at  very  high  energies  fBl  + 
fB2p  does  not  approach  fBl  for  any  value  of  the  momentum  trans¬ 
fer  |K|.  In  1959  Glauber16  employed  the  eikonal  approximation17 
of  potential  scattering  theory  to  discuss  the  many-particle  scatter¬ 
ing  process.  To  obtain  the  eikonal  wave  function  for  the  projectile 
one  introduces  a  new  variable  p  defined  by  p  =  k,  -  q  in  (5)  and 
approximately  writes 

1  ^  1  p 2 

k-  -  q2  +  it  2k,p  +  it  (2k,-p  +  it)2'  { 

If  we  put  only  the  first  term  of  (15)  in  (5)  we  get  the  linearized 
Green’s  function  of  the  projectile  leading  to  the  eikonal  wave 
function.  However,  the  inclusion  of  the  second  term  gives  the 
leading  Wallace  phase  correction.18  Use  of  the  linearized  Green's 
function  in  (1)  yields  the  scattering  amplitude  in  the  Glauber  ap¬ 
proximation,  given  by 

fc%X)  =  </!/e(x)|0  (16) 

where  /£(x)  is  the  scattering  amplitude  in  the  eikonal  approxi¬ 
mation  for  the  scattering  of  a  particle  by  a  target  atom  in  which 
the  bound  electrons  are  frozen  at  x.  Similarities  between  (11)  and 
(15)  as  well  as  between  (6)  and  (16)  are  to  be  noted.  If  we  include 
the  leading  Wallace  phase  correction. f^x)  changes  tofEW(x)  given 
by19 

/ew(X)  = 

—Z  I (Pb  exp(t'K-b)  {exp [t(£f  'xo(M)  +  fc,-3Xt(M))l  - 1}.  (17) 

K  is,  as  above,  the  momentum  transfer.  The  eikonal  phase  Xo(b’x) 
is  of  first  order  and  the  leading  Wallace  phase  correction  Xi(b-x) 
is19  of  second  order  in  the  interaction  potential  U.  Of  course,  one 
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gets  fE(\)  on  dropping  the  term  containing  Xi-  We  note  that  the 
Glauber  approximation  with  or  without  the  Wallace  correction  com¬ 
pletely  neglects  the  effect  due  to  the  distortion  of  the  target  wave 
function.  fc  as  well  as  fw.  defined  by 

fw(  k/.k,)  =  </|/£w(x)|/>.  (18) 

suffer  from  the  same  discrepancies  as  noticed  earlier  for  fBz.u-  The 
scattering  amplitudes  fG  and  fw  can  each  be  expanded,  in  a  manner 
analogous  to  the  Born  series,  in  powers  of  U  to  yield  the  Glauber 
series  and  the  Wallace  series,  respectively.  ~fGn  and  fwn  will  refer 
to  terms  of  order  n  in  the  Glauber  and  in  the  Wallace  series, 
respectively.  However,  whereas  the  Glauber  series  terms  are  al¬ 
ternately  real  and  imaginary,  the  terms  of  the  Wallace  series  are 
complex  (for  n  >  1).  Furthermore,  at  large  momentum  transfers, 
the  terms  of  the  Wallace  series  are  nearly  equal  to  the  correspond¬ 
ing  terms  of  the  Born  series.  Byron  et  al.1' 9  took  advantage  of  the 
above  property  and  proposed  a  unitarized  eikonal  Born  amplitude 
given  by 


/uebs  -  fw  ~  ?W2  +  ?B2  (19) 

where  is  the  second  term  of  the  Wallace  series.  However, 
difficulties  arise  in  the  evaluation  of  the  Wallace  terms  JWn  (with 
n  >  4)  of  the  e  =  -atom  scattering  amplitude.  This  led  Byron  et  al.1 
to  redefine  the  unitarized  eikonal  Born  series  in  the  following 
manner: 


/uebs  ~  fw  ~  ?W2  +  Ib2>  (20) 

where  fw  is  again  obtained  from  (18)  except  for  the  difference 
that  in  the  evaluation  of  /£VV(x)  the  phase  term  exp[/fc,_:,Xi(b-x)]  is 
replaced  by  1  +  ifc,_3X i(b,x).  We  note  that  the  eikonal  Born  series 
as  well  as  the  unitarized  eikonal  Born  series  include  the  effect  due 
to  the  distortion  of  the  target  wave  function  only  up  to  the  second 
order  in  the  interaction  potential  through  }&,  which  is  usually 
evaluated  following  the  procedure  of  Holt  and  Moisei witsch. 10 

Recently  Byron  et  al.2  utilized  the  UEBS  method  to  obtain  the 
cross  sections  for  the  elastic  and  inelastic  scattering  of  electrons 
and  positrons  by  atomic  hydrogen  at  intermediate  and  high  ener- 
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gies.  At  about  the  same  time  Khare  and  Prakash3  utilized  the 
Schwinger  variational  method  to  investigate  the  e  =  -H  scattering. 
We  now  explicitly  demonstrate  that  the  scattering  amplitude  in  the 
Schwinger  variational  approach  follows  by  approximating  the  Green's 
function  in  Eq.  (1).  On  splitting  the  infinite  sum  in  Eq.  (1)  at  n 
=  p.  the  scattering  amplitude  can  be  written  in  the  following 
alternate  form: 


//,(k/,kl)  =  fBp  -  2tt2  (kfJ\U(GZU)p'Z  (Ctft/y-'M.  (21) 


in  which  fBp  is  the  pth  Born  approximation  for  the  scattering  am¬ 
plitude.  The  Green’s  function  appears  explicitly  in  the  infinite  sums 
of  Eqs.  (1)  and  (21).  Now  we  make  the  approximation  of  truncating 
these  two  infinite  sums  to  finite  sums  of  m  terms  each  and  replacing 
the  Green’s  function  GJ  in  the  surviving  m  terms  by  a  multiple 
CGq  of  the  Green's  function,  with  the  complex  multiplying  factor 
C  to  be  determined  later.  On  using  this  approximation  for  the 
Green’s  function  in  Eqs.  (1)  and  (21).  we  get  two  approximate 
expressions  for  the  scattering  amplitude  which  we  label  as  f,  and 
f2,  respectively.  These  are 


fi  -  G  fBm 


and 


fl  ~  f Bp  +  G{ fflpTm  fBp)- 

The  multiplying  factor  C  is  now  determined  by  simply  equating  /, 
and  f2,  as  they  would  have  been  the  same  if  the  Green's  function 
was  not  approximated.  This  procedure  yields 


C  = 


/a? 


f  Bm  f  Bp  f  Bp  ■*  m 


(22) 


On  substituting  this  expression  for  C  in  either/!  or  in  f2  above, 
we  obtain  for  the  scattering  amplitude 


If  1  = 

l/pmj 


f Bmf Bp 


fBm  "b  /fin  f Bp  -f  m 


(23) 
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The  same  expression  tor  the  scattering  amplitude  is  obtained3  if 
one  takes  incoming  and  outgoing  scattered  vvaves  correct  to  (p  - 
l)thand(«  -  l)th  order  in  the  interaction  potential  in  the  Schwin¬ 
ger  variational  principle.  Khare  and  Prakash3  took  m  =  p  =  2, 
replaced  /S4  by  fCi.  and  noting  that  fai{ 0)  is  zero  in  the  closure 
approximation  for  e--H  scattering,-0  found 


l/d  - 


Qb 1  +  JbiY 

/si  +  )b2  ~  fc 4 


(24) 


for  the  scattering  amplitude  in  the  forward  direction.  The  use  of 
(24)  along  with  the  optica!  theorem  gave  slightly  different  values 
of  the  total  cress  sections  QT  for  electron  and  positron  scattering 
at  lower  impact  energies  even  when  exchange  was  excluded.  On 
the  other  hand  the  UEBS  method  of  Byron  e:  ri.1  gave  different 
values  of  Qr(e~)  only  when  exchange  was  included.  In  Table  I  we 
show  the  two  sets  of  values  which  are  in  good  agreement  with  each 
other,  'particularly  for  positron  scattering  where  exchange  does  not 
play  any  role.  For  the  electron  case  the  difference  between  the 
two  sets  of  values  seems  to  arise  mainly  due  to  the  exchange  con¬ 
tribution.  Such  an  agreement  is  not  surprising  because  [/^(O)]  and 
fv hbs(O)  agree  with  one  another  asymptotically  up  to  the  order  of 
k'3.  A  comparison  of  the  theoretical  values  of  the  cress  sections 
for  the  electrons  with  the  adopted  cross  sections  of  do  Heer  ei 


TABLE i 


Total  coilisional  cross  sections  (in  ai)  for  the  scattering  of  electrons  and 
positrons  by  a  hydrogen  atom 


£(eV) 

SV* 

UEBS1’ 

Hc 

e’ 

e" 

C  ' 

100 

6.77 

6.82 

7,04 

6.84 

6.85 

209 

4.18 

4.18 

4.23 

4.18 

4.18 

300 

2  93 

2.93 

J.10 

3  07' 

306 

400 

2,43 

2.43 

2.45 

2.44 

2. 45 

500 

2.03 

2.03 

- 

- 

*SV  are  the  results  of  Khare  and  Prakash  (Ref.  3'  who  employed  the  Schwinger 
variational  method,  Eq.  (23),  with  m  =  p  =  2. 

"UEBS  are  the  results  of  Byron  et  al.  (Refs.  2  and  19)  who  employed  the  um- 
tarized  eikonal  Bom  series  method. 

'H  are  the  adopted  cross  sections  of  de  Heer  -tt  al.  (Ref.  21). 
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al.  ,:i  presented  in  Table  I.  shows  that  the  agreement  between  them 
is  quite  satisfactory.  It  may  be  noted  that  de  Heer  et  al.  used 
theoretical  as  well  as  experimental  cross  sections  for  elastic  and 
inelastic  processes  from  various  different  sources  to  obtain  the  total 
collisional  cross  sections. 

For  a  better  comparison  of  the  two  methods — the  Schwinger 
variational  and  the  UEBS — one  should  compare  the  values  of  the 
differential  cross  sections  rather  than  the  total  collisional  cross 
sections.  However,  at  nonzero  scattering  angles  ~f8 3(0)  is  not  zero 
and  no  tractable  method  is  so  far  available  for  obtaining  Im  fB 3(0),. 
which  falls  as  k,~3.  Hence  Khare  and  Lata4  took  m  -  2  and  p  =  1 
in  (23),  replaced  Re/fl3(0)  by  /C3(0)  and  neglected  Im  /83(0).  Thus 
they  obtained  [/2i(0)j  which  agrees  with  AjebsW  asymptotically 
only  up  to  the  order  of  k~2.  Furthermore,  Khare  and  Lata  obtained 
the  exchange  contribution  through  the  Ochkur  approximation 
whereas  Byron  et  al.  obtained  the  exchange  contribution  through 
a  more  sophisticated  method.  Tables  II  and  III  compare  the  two 
sets  of  cross  sections  for  the  elastic  scattering  of  e-  by  hydrogen 
atoms  in  the  energy  region  of  100  to  400  eV.  At  the  highest  impact 
energy,  that  is,  400  eV.  the  agreement  between  the  two  sets  is  very 
good.  The  maximum  difference  is  only  3.6%.  However,  with  the 
decrease  of  the  impact  energy  the  difference  between  the  two  sets 
of  cross  sections  increases.  At  100  eV  the  maximum  difference  is 
about  25%.  Such  behavior  is  expected  because  the  scattering  am¬ 
plitudes  obtained  by  the  two  methods  agree  with  one  another  only 
up  through  terms  of  order  k~2.  In  general  the  differences  between 
the  two  sets  of  values  are  relatively  larger  at  higher  scattering 
angles  and  both  sets  of  values  are  smaller  than  the  experimental 
values  of  Williams22  and  van  Wingerden  et  alP  At  larger  angles 
though,  /uebsW  yields  closer  agreement  with  the  experimental 
data  of  Williams.22  However,  we  note  that  Byron  et  al.  have  com¬ 
mented  that  so  far  no  theoretical  method  has  given  excellent  agree¬ 
ment  with  the  experimental  data  of  Williams  at  large  scattering 
angles.  According  to  Kingston  and  Walters,24  the  data  of  Williams 
are  consistently  larger  at  larger  angles.  Perhaps  the  same  could  be 
said  for  the  data  of  van  Wingerden  et  al. 

In  summary,  various  approximate  forms  of  the  Green’s  functions 
have  been  shown  to  give  rise  to  practically  ail  the  theoretical  meth¬ 
ods  used  for  investigations  of  scattering  of  fast  projectiles.  In  their 
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latest  method  Byron  et  al.  have  put  emphasis  on  the  unitarization 
whereas  Khare  and  Prakash  based  their  method  on  the  Schwinger 
variational  principle.  Both  of  these  methods  seem  to  be  attractive 
and  their  application  to  different  processes  for  heavier  elements 
and  the  comparison  of  the  resulting  cross  sections  will  be  of  future 
interest. 
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Exact  time-dependent  evolution  of  electron-velocity  distribution  functions 
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A  numerical  technique,  starting  from  the  Boltzmann  equation,  for  obtaining  the  time-dependent 
behavior  of  the  electron-velocity  distribution  function  in  a  gas  is  presented.  A  unique  feature  of  this 
technique  is  that,  unlike  previously  used  procedures,  it  docs  not  make  use  of  any  expansion  of  the 
distribution  function.  This  allows  the  full  anisotropy  of  the  distribution  function  to  be  included  in 
the  solution.  Furthermore,  the  problem  associated  with  multiterm-expansion  techniques  of  choos¬ 
ing  a  sufficient  number  of  terms  for  convergence  is  completely  avoided.  The  distribution  function 
obtained  by  the  present  method  is  exact  and,  in  principle,  contains  all  of  the  expansion  terms  of  the 
previous  procedures.  Details  of  the  algorithm,  including  stability  conditions,  treatment  of  the  boun¬ 
daries,  and  evaluation  of  (he  collision  integrals,  are  presented.  This  technique  has  been  applied  for 
obtaining  the  time-dependent  behavior  of  electron  swarms  in  gaseous  argon  and  neon  for  various 
values  of  E/N  (the  ratio  of  the  applied  uniform  dc  field  to  the  gas  density),  and  the  corresponding 
results  arc  presented. 


I.  INTRODUCTION 

The  electron-velocity  distribution  function  (EVDF)  is 
fundamentally  important  in  virtually  all  aspects  of  gase¬ 
ous  electronics.  The  EVDF  provides  a  statistical  descrip¬ 
tion  of  the  motion  of  all  of  the  electrons  in  an  electron 
swarm.  The  motion  of  the  electrons  in  the  swarm  is 
affected  by  externally  applied  electric  and  magnetic  fields, 
and  by  collisions  of  the  electrons  with  the  particles  of  the 
ambient  gas.  These  external  forces  and  collisions  cause 
tune-dependent  changes  in  the  EVDF.  Stating  this  pro¬ 
cess  mathematically,  let  /(v,f)  represent  the  electron- 
velocity  distribution  function  at  a  velocity  v,  and  at  a 
particular  time  t.  Then,  at  some  later  tune  t  +  At,  the 
EVDF  can  be  described  very  simply  by 

/(v  +  Av,t  +  At)=/(v,0  +  R(v,l)At  .  (1) 

Here,  Av=aA /,  with  a  as  the  acceleration  of  the  electrons 
due  to  the  externally  applied  forces.  R  (v,t)At  is  the  col¬ 
lision  term  which  represents  the  net  change  in  /(v,t) 
during  the  time  increment  At  due  to  all  possible  collision 
processes  between  the  electrons  and  the  gas  particles. 
Now,  if  we  expand  Eq.  (1)  to  first  order  in  At,  and  then 
take  the  limit  as  At  goes  to  zero,  the  spatially  indepen¬ 
dent  Boltzmann  equation  is  immediately  obtained: 

+a ■V„/<v,U=K(v,/>  .  (2) 

It  is  thus  clear  that  the  physical  content  of  the 
Boltzmann  equation  (2)  is  entirely  equivalent  to  that  of 
the  difference  equation  (1).  Knowledge  of  the  EVDF  is 
usually  gained  by  solution  of  the  Boltzmann  equation 
and,  to  this  end,  many  techniques  for  its  solution  have 
been  developed. 

Traditionally,  the  techniques  used  for  solving  the 
Boltzmann  equation  for  an  equilibrium  electron-velocity 
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distribution  function  have  involved  expansion,  usually  in 
the  Legendre  polynomials,  of  the  distribution  function  as 
follows: 

/( v,r)=  2  fn(o,t)P„(0) .  (3) 

n  “0 

Often,  only  the  first  two  terms  in  the  expansion,  contain¬ 
ing  /0  and  /, ,  are  retained,  and  the  time  derivative  of  the 
distribution  function  in  the  Boltzmann  equation  is  set  to 
zero  to  correspond  to  the  equilibrium  situation.  The  re¬ 
sulting  coupled  time-independent  equations  can  then  be 
solved  for  /0  and  / (  using  standard  numerical  methods. 
The  two-term  expansion  method,  however,  breaks  down1 
under  situations  of  large  E/N  (the  ratio  of  the  applied 
uniform  dc  electric  field  to  the  gas  density)  or  for  cases  in 
which  the  inelastic  scattering  cross  sections  are  compara¬ 
ble  in  magnitude  to  the  elastic  cross  sections.  The 
shortcomings  of  the  two-term  expansion  can,  in  principle, 
be  overcome  by  retaining  more  terms  in  the  Legendre  ex¬ 
pansion  of  the  distribution  function.  These  multiterm 
methods,  however,  have  their  own  drawbacks.  As  more 
terms  in  the  expansion  are  kept,  the  computational  com¬ 
plexity  increases  rapidly.  Furthermore,  cross  sections 
which  have  been  adjusted  to  reproduce  experimental 
swarm  parameters  in  a  two-term  expansion  calculation 
do  not  yield  the  same  EVDF  or  the  corresponding  swarm 
parameters  when  used  in  a  multiterm  expansion  calcula¬ 
tion  and  vice  versa.  Clearly,  it  would  be  desirable  to  have 
a  procedure,  as  described  below,  which  can  provide  the 
equilibrium  EVDF  without  involving  any  expansion  of 
the  distribution  function.  Such  a  procedure,  which  in¬ 
corporates  a  finite-difference  technique,  was  developed  by 
Tagashira  and  co-workers.2  In  their  procedure,  the  dis¬ 
tribution  function  was  expanded  to  second  order  in  time 
using  a  standard  Taylor  series.  The  various  time  deriva¬ 
tives  of  /(v,r)  were  evaluated  by  direct  substitution  from 
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the  Boltzmann  equation,  namely, 

/(v,r  +A  <)=/(v,r)+^V,-At  +0((  At)2) 

at 

=  /<v,/)+[-aV„/(v,t)+K(v,r)]Ar 

+  0((Ar)2)  .  (4) 

Tagashira  and  co-workers  chose  to  evaluate  the  distribu¬ 
tion  function  of  Eq.  (4)  in  spherical  coordinates  in  the  ve¬ 
locity  space,  that  is,  /  was  stored  as  a  v-0  array,  simply 
because  the  evaluation  of  the  collision  term  R(\,t )  is 
most  convenient  in  spherical  coordinates.  Evaluation  of 
the  a-V„  term  then  involved  derivatives  of  /  with  respect 
to  t)  and  0.  These  derivatives,  which  had  to  be  taken  nu¬ 
merically  by  a  finite-difference  procedure,  are  prone  to  in¬ 
stability.  It  was  in  order  to  alleviate  this  instability  that 
Tagashira  and  co-workers  had  to  retain  some  of  the 
terms  proportional  to  (At)2  in  the  Taylor-series  expansion 
of  Eq.  (4).  In  a  previous  paper,3  we  briefly  described  a 
finite-difference  algorithm  for  determining  the  exact 
time-dependent  behavior  of  electron-velocity  distribution 
functions.  A  unique  feature  of  our  algorithm  is  that  it 
does  not  require  numerical  evaluation  of  any  derivatives, 
nor  does  it  make  use  of  any  term  expansions  of '!  '  distri¬ 
bution  function  in  terms  of  Legendre  functions.  The 
present  paper  will  provide  details  of  the  procedure,  such 
as  the  explicit  form  and  evaluation  of  the  collision  in¬ 
tegrals,  conditions  of  numerical  stability,  treatment  of  the 
_ I 


numerical  boundaries,  and  techniques  for  implementing 
the  conditions  of  numerical  stability.  We  will  also 
present  results  of  the  application  of  this  procedure  to 
electron  swarms  in  argon  and  in  neon  for  various  values 
of  E/N. 

In  the  present  solution,  the  following  finite-difference 
equation  in  Cartesian  coordinates  for  the  electron- 
velocity  distribution  function  is  evaluated 

/(u^.^.Uj+ADj,/  +  At) 

=/(ut,i)J,,oI,r)+/?(ut,uJ,,i»J,/)Af  .  (5) 

Equation  (5)  can  be  obtained  directly  from  Eq.  (1),  which 
is  equivalent,  in  its  physical  content,  to  the  Boltzmann 
equation  (2).  Until  the  collision  terms  R  (v,r)  are  known, 
however,  Eq.  (5)  is  of  little  practical  use.  So,  before 
proceeding  any  further,  we  will  explicitly  define  the  col¬ 
lision  terms  R  (v,t)  and  outline  procedures  for  their  eval¬ 
uation. 

II.  COLLISION  TERMS 

In  order  to  derive  a  very  general  expression  for  the  col¬ 
lision  term  we  ’will  assume  an  ambient  gas  of  constant 
density  interacting  with  a  spatially  homogeneous  swarm 
of  projectiles  of  arbitrary  mass  (for  example,  either  elec¬ 
trons,  positrons,  protons,  or  ions)  A  general  expression 
for  the  collision  term  in  Eq.  (5)  can  be  written4 


R(v,t)=  J  (N/v1)  J* Vpdvp  f^sinxpd tl>  f*”davpf(vp,t)crp{vr,’p)5{v  -gp(vp,ip))-Nvf(v,t)ar(v)  (6) 


where  op(vp,tl>)  is  the  differential  scattering  cross  section 
for  the  pth  ( p  ^elastic,  excitation,  ionization,  etc.)  col¬ 
lision  process,  and  crT(v)  is  the  integrated  total  cross  sec¬ 
tion  for  all  collision  processes.  The  function  gp(vp,<l>)  is 
defined  by  the  equation 


I 

the  convenience,  with  respect  to  either  z  or  v.  These  two 
possible  specifications  of  vp  lead  to  the  following  relation¬ 
ships  among  various  angles: 

v  =  v(u,0)  ,  (9a) 


v=gp(»p,'l>) ,  <7) 

which  relates,  via  the  energy-conserving  8  function,  the 
initial  speed  vp  to  the  final  speed  v  for  the  pth  collision 
process.  The  integral  terms  in  Eq.  (6)  represent  the  rate 
at  which  the  projectile  particles  are  scattered  into  a 
velocity-space-volume  element  dlv  located  about  v  due  to 
the  pth  scattering  process  and  will  be  denoted  by 
R  +(v,r).  The  last  term  in  Eq.  (6)  represents  the  rate  at 
which  the  projectile  particles  are  scattered  out  of  d3v 
about  v  due  to  all  possible  collision  processes,  and  will  be 
denoted  by  flf(v,/).  Thus 


K<v,f)=2tf/<v,f)-Kf<v,r)  ■ 


are  defined  in  Fig.  1  and  by  the  following  relations: 
d3v=v2dvsir.'9)d0d<f>  ,  (8a) 

d}vp=vpdvpsin(0p)d9pd<l>p  =  vpdvpsm(4i)dtlida  ,  (8b) 

so  that  the  direction  of  v  can  be  specified,  depending  on 


FIG.  I  Geometry  used  in  the  derivation  of  the  collision 
terms,  v  the  vp  denote  the  final  and  initial  velocities,  respective¬ 
ly,  and  ij/  is  the  scattering  angle. 
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cos(0p)=sin(0)sin(</')cos(«)+cos(0)cos(t/>)  ,  (9b) 

y P=y p(vt><Op)—'/p(Vp}0,'l>,a)  ,  (9c) 

and  finally, 

fl'/p,t)=:flvp,Op,t)=f{vp,0,il>,a,t) .  (9d) 

In  order  to  evaluate  Eq.  (6)  the  function  gp(up,^)  must  be 

specified  for  various  collision  processes.  For  clastic  col¬ 
lisions  (p  =e),  straightforward  kinematics  yields 

v  =veq(4l)=ge(v„tl>)  ,  (10a) 

where  q  (tp)  is  the  following  function: 

9(^)=[(l-/r2)l/2+/<cost/>]/(l+/r)  .  (10b) 


Here,  ft  is  the  ratio  of  the  mass  m  of  the  projectile  and 
the  mass  M  of  the  gas  particle,  that  is,  fi~m  /M.  For  the 
case  in  which  the  projectile  is  cither  an  electron  or  a  posi¬ 
tron  (/i « 1 ),  q(tp)  can  be  simplified  to 

q(i/')=[l—  /t(l~COS0)]  .  (10c) 

For  inelastic  collisions  (p  =/)  the  energy-conserving  ex¬ 
pression  for  gp  is 

v  =(t>2— 2£/m)1/2=g,(i), ) ,  (11) 

which  does  not  depend  on  the  angle  i)i.  §  is  the  energy 
loss  associated  with  the  inelastic  process.  Performing  the 
radial  vp  integrations  to  eliminate  the  Dirac  8  functions 
in  Eq.  (6),  one  is  left  with  a  two-dimensional  angular  in¬ 
tegral  over  tp  and  a.  Because  of  the  8  function,  vp  in  the 
integrand,  which  includes  the  distribution  function  /  and 
the  collision  cross  section  ap,  is  replaced  by  11/9(1/1)  (for 
the  elastic  part)  or  by  (w24- 2f/m)l/2  (for  the  inelastic 
part).  For  brevity  in  writing  and  in  accordance  with  Eqs. 
(10!  and  (11),  we  replace,  111  the  resulting  angular  in¬ 
tegrals,  v/q(tp)  and  (ir+2f/;n)l/2  with  and  v,,  re¬ 
spectively,  so  that  the  final  expression  for  the  collision 
term  looks  like 

y  4 

R  (v,r)=  N  f-Ara^v^tpmvpdtp  f2“/(vr,t)da 
J  0  J  0 

u2 

(12a) 

So  far,  the  collision  integral  is  valid  for  any  type  of  pro¬ 
jectile  except  that  the  term  corresponding  to  the  ioniza¬ 
tion  process  has  to  be  treated  slightly  differently  when  the 
projectile  is  an  electron.  When  ionization  is  considered, 
the  final  energy  of  both  the  incident  electron  and  the  free 
electron  that  is  produced  via  the  ionization  process  must 
be  properly  accounted  for.  To  this  end,  an  electron- 
energy  partition  ratio  A /( I  —  A)  is  used,  which  denotes 
the  ratio  of  the  available  energy  that  goes  to  each  of  the 
two  electrons  (labeled  1  and  2  below).  The  integral  that 
represents  the  rate  at  which  electrons  scatter  due  to  the 
ionization  process,  into  the  velocity  space  element  d1  v 


about  (u ,0),  then  has  two  parts: 

R*nA(v,0)+R^(v,O) 

~  ^f*-^°,o'Sv\>'l’)s'''’!>dtl>f2"f(vut)da 

(12b) 

where  U|  =(2g/m  +p2/A)i/2  and  i>j  =[2f /m  -t-u2/(  1 
—  A))l/2.  These  two  terms  correspond  to  two  Afferent 
ionization  events  in  which  electrons  of  initial  speeds  t)| 
and  t>2  ionize  gas  particles.  One  of  the  two  electrons 
[having  energy  ratio  A/(l—  A)]  resulting  from  each  of 
these  two  ionization  events  has  speed  u. 

In  evaluating  Eqs.  (6)  or  (12),  proper  account  must  be 
taken  of  the  vector  nature  of  v  and  vp.  Although  the  dis¬ 
tribution  function  is  symmetric  about  the  axis  parallel  to 
the  electric  field  (the  z  axis),  the  polar  axis  of  the  integrals 
over  a  solid  angle  ( \p,a )  is  tilted  with  respect  to  z  and 
cannot  take  advantage  of  the  symmetry. 

The  surface  over  which  the  distribution  function  in  the 
elastic  part  of  Eq.  (12a)  is  evaluated  for  angular  integra¬ 
tion  is  represented  graphically  in  Fig.  2.  Note  that,  be¬ 
cause  Df  depends  on  tp,  the  surface  is  not  a  perfect  sphere 
but  an  "egg-shaped”  surface  with  azimuthal  symmetry 
about  v.  More  specifically,  such  a  surface  is  realized  by 
the  tracing  of  the  lip  of  a  vector  whose  length  increases 
continuously  as  the  polar  angle  is  vg-ied  from  one  pole 
( 0=0)  to  the  other  pole  (t/t=ir),  but  the  length  of  the 
vector  remains  fixed  as  the  azimuthal  angle  a  is  varied, 
for  a  given  ip,  from  0  to  2ir.  The  deviation  of  this  surface 
from  a  perfect  sphere  is  related,  from  Eq.  (10),  to  the 
mass  ratio  ft.  If  the  projectile  particles  are  electrons  or 
positrons  the  surface  is  very  nearly  a  sphere,  and  very  lit¬ 
tle  difference  was  found  on  simply  replacing  vc  by  u  in  Eq. 
(12a).  However,  in  order  to  keep  the  analysis  more  gen¬ 
eral  and  not  limit  the  distribution  functions  to  only  those 
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FIG.  2.  Surface  over  which  the  elastic  component  of  the  col¬ 
lision  term  15  evaluated.  0  is  the  polar  angle  about  z,  and  $  1S 
the  polar  angle  about  v. 
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of  electrons  or  positrons,  we  will  not  make  any  assump¬ 
tions  about  the  relative  masses  of  the  gas  molecules  and 
the  projectiles.  Nevertheless,  Eq.  (12a)  remains  valid  in 
all  possible  cases  including  the  ones  in  which  the  projec¬ 
tile  particles  are  more  massive,  such  as  protons  or  ions, 
for  which  the  distortion  and  tilt  of  the  surface  would  be¬ 
come  very  important. 

Often,  the  differential  scattering  cross  sections  are  not 
sufficiently  available  to  cover  all  angles  and  energies  for 
each  scattering  process  in  the  system  of  interest.  In¬ 
tegrated  cross  sections,  however,  can  usually  be  found 
and  may  be  the  only  alternative.  In  this  case,  or  in  the 
cases  in  which  the  scattering  is  not  strongly  dependent  on 
angle,  we  may  make  the  approximation  of  isotropic 
scattering  (a  docs  not  depend  on  ip).  For  isotropic  elastic 
scattering  we  replace  the  differential  elastic  scattering 
cross  section  or(vr,\ !’)  with  of(vt )/4tt,  where  o,(t>,)  is 
the  integrated  elastic  scattering  cross  section.  Doing  this, 
the  elastic  part  of  Eq.  (12)  can  be  rewritten  as 

R*iv,t)=~^-r  f"v*ae(vc)s\ni!)dxlifi’!f(vr,i)da  . 

4iru3  jo  jo 

03) 

If  we  assume  isotropic  inelastic  scattering,  then  the 
differential  scattering  cross  section  for  the  ith  scattering 
process  (r,(t),,(M  can  be  replaced  by  a,(v,  )/4rr,  where 
<r,( v, )  is  the  integrated  cross  section  for  the  ith  inelastic 
process.  Also,  since  v,  does  not  depend  on  the  scattering 
angle,  the  surface  of  integration  is  spherical  with  a  con¬ 
stant  radius  t’,,  and  the  terms  containing  e,  can  be  taken 
outside  of  the  angular  integral.  We  can  then  write 


A’d,2 

R^iv,!)-- — 0,(0,; 
4 17V 

I  fjsi  n  ifi  dil>  /J  V<  v, ,  t  Ida 

(14a) 

Ah),2 

=  ~ - (7,(l),, 

4 771) 

)fjsm0,d0:  f^f(v„t)d<pi 

(14b> 

Au,2 

=  —<7,(1), 

)f*f(o„0„t)s\n6ld0l  . 

(14c) 

In  going  from  Eq.  (14a)  to  (14b)  we  have  made  use  of  the 
fact  that  the  integral  is  independent  of  the  choice  of  polar 
axis  and  we  have  chosen  the  z  axis  (which  is  the  axis  of 
symmetry  for  the  EVDF)  as  the  polar  axis,  so  that  the  0, 
integration  becomes  trivial  [see  Eq.  (8b)]. 

Another  problem  that  must  be  considered  in  the  evalu¬ 
ation  of  Eq.  (12a)  is  that  /  is  stored  in  a  rectangular  array 
Gn  the  x-z  plane)  and  has  rectangular  boundaries,  but  the 
surfaces  of  integration  in  the  R  +  terms  are  nonrcctangu- 
lar  and  sometimes  lie  outside  of  the  region  in  which  /  is 
known.  To  handle  this  problem,  an  extrapolation  pro¬ 
cedure  has  to  be  devised.  Regrettably,  extrapolations  are 
a  risky  business,  and  one  can  only  hope  that  the  relevant 
numerical  errors  will  be  small.  What  wc  did  was  to  set 
up  the  initial  ( t  =  0)  Maxwellian  distribution  of  projec¬ 
tiles  (electrons  in  our  actual  calculations)  so  that  the 
values  of  the  distribution  function  near  the  boundaries 
were  less  than  10-}  of  the  peak  value  so  that  the  bound¬ 
ary  contributions  would  be  small.  Then  we  assumed  that 
the  behavior  in  the  high-velocity  regions  (near  the  boun¬ 


daries  and  beyond)  remained  essentially  Maxwellian  at  all 
later  times.  Based  on  the  shape  of  the  distribution  func¬ 
tion  near  but  inside  the  boundary,  wc  extrapolated  the 
distribution  function  using  the  simple  recursion  relation 
f(v  +  2Au)  =  [/2(i)  +  Au)//(u)]C,  where  C  is  a  constant 
that  depends  on  A v.  This  recursion  relation  follows  from 
the  assumed  Maxwellian  form  of  the  distribution  function 
near  the  boundary.  The  case  C  =  I  corresponds  to  the 
logarithmic  approximation  for  extrapolation  used  by  pre¬ 
vious  investigators.2  Several  tests  were  made,  by  chang¬ 
ing  the  boundaries,  showing  that  the  high-velocity  tail  of 
the  numerically  obtained  distribution  function  conformed 
to  this  Maxwellian  behavior,  and  that  deviation  was  less 
than  1%. 

III.  METHOD  OI'  SOI.UTION 

For  the  situation  in  winch  the  external  force  on  the 
projectiles  (of  charge  q)  is  provided  by  a  uniform  dc  elec¬ 
tric  field  aligned  in  the  z  direction,  the  constant 
Ar. =(qE/in)At,  and  Eq.  (5)  can  be  rewritten  as 

/( n  < ,  v.  +  ( qE  /m )  A/,  t  +  A  t ) 

—/(u,, t>}  vx,t)+R (vxvv,v.,t)At  .  (15) 

Equation  (15)  is  well  sailed  for  evaluation  on  a  comput¬ 
er.  Since  there  is  axial  symmetry  around  the  z  axis,  / 
need  only  be  stored  as  a  function  of  v.  and  vx  (or  vy )  in 
such  a  way  that  the  velocity  increments  Ai>  satisfy  the  re¬ 
lation  Au  =  (qE/m) Ar.  Evaluation  of  /  in  Eq.  (15)  then 
merely  involves  a  shifting  of  the  two-dimensional  array 
/( t\,Dj)  along  Uj  at  each  time  interval  At,  and  then  add¬ 
ing  to  each  array  element  the  corresponding  collision 
term  R  (vx,vx,t)At.  This  shifting  procedure  accomplishes 
all  the  acceleration  effects  of  the  projectiles  due  to  the 
electric  field  and  is  inhcicntly  immune  to  round-off  error. 
Carrying  out  this  procedure  will  require  knowledge  of  the 
collision  integrals  for  each  vx  and  i).  at  time  t .  Evalua¬ 
tion  of  these  integrals  was  described  in  Sec.  II,  and  in¬ 
volves  an  integral  over  the  polar  angle  0  (in  velocity 
space)  and  thus  requires  a  knowledge  of  the  distribution 
function  at  various  values  of  v  and  0.  This  integration 
can  be  carried  out,  even  though  /  is  known  only  as  a 
function  ofu,  and  vx,  by  simply  interpolating /(t),, vzl  to 
get  f(v,0). 

The  algorithm  for  obtaining  the  time  evolution  of  the 
velocity  distribution  function  is  as  follows. 

(i)  Store  an  initial  distribution  function  (for  instance,  a 
Maxwellian  at  t  ~0)  in  a  two-dimensional  array  f(vx,vz) 
such  that  Au  =(qE/m)At. 

(ii)  From  the  existing  distribution  function,  evaluate 
the  collision  terms  R  ( vx ,  v, )  for  each  vx  and  v,. 

(iii)  Multiply  each  of  the  collision  terms  by  At  and  add 
to  the  corresponding  distribution  function  array  element 
[f(vx,vx)—rf(vx,v.)+R  (Uj.UjlAt]. 

(iv)  Shift  the  resulting  array  along  the  t>.  index 
[f(vx, vz)-rf{vx,vz  +  Avz)]  to  obtain  a  new  distribution 
function  which  corresponds  to  time  t  4-  At. 

(v)  Go  to  step  (ii). 

The  procedure  outlined  in  steps  (ii)-(v)  are  repeated 
while  various  swarm  parameters  are  calculated  from  the 
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distribution  function  obtained  in  step  (iv)  for  each  time 
cycle.  Equilibrium  is  obtained  when  the  swarm  parame¬ 
ters  cease  to  change  in  time. 

IV.  STABILITY  CONDITIONS 

In  any  numerical  solution  to  a  differential  equation, 
stability  is  always  a  major  consideration.  The  procedure 
of  shifting  the  array  elements  along  the  vz  index  accom¬ 
plishes  all  of  the  acceleration  effects  due  to  the  electric 
field  and  avoids  the  need  to  evaluate  any  derivatives  nu¬ 
merically.  This  procedure  in  itself  greatly  enhances  the 
numerical  stability  of  the  calculation;  however,  it  is  not  in 
itself  sufficient.  One  of  the  conditions  of  stability  is  a  re¬ 
statement  of  the  so-called  Courant-Friedrichs-Lcwys 
(CFL)  condition  and  is  simply 

(nArl/Au  <  1  ,  (16) 

where  a=qE/m  is,  in  this  case,  the  acceleration  of  the 
projectiles  due  to  the  electric  field  E.  Since  /  is  stored  so 
that  its  velocity  increments  satisfy  An  =  (qE/m)At,  Eq. 
(16)  is  always  minimally  satisfied  as  an  equality.  The  cal¬ 
culation  can  satisfy  Eq.  (16)  more  strongly  by  choosing  a 
smaller  time  increment  At',  which  is  a  proper  fraction  of 
At  such  that 

At'-Al/n ,  Av'-Av/it  ,  (17) 

while  leaving  a  and  Av  unchanged.  Then  Eq  (5)  can  be 
rewritten 

/(v  +  Av,t  +Af')=/(v  +  Av-Av',t) 


FIG.  3  Schematic  representation  of  the  acceleration  pro¬ 
cedure  described  in  Eq  119),  sh  awing  the  implementation  of  the 
Courant-Fricdrichs-I.cwy  stability  condition 

velocities  in  order  to  still  be  satisfied  for  high  velocities. 
If  the  adiabatic  condition,  Eq.  (20),  is  not  met,  an  obvious 
consequence  is  that,  if  an  array  clement  of  the  collision 
term  is  negative,  the  corresponding  array  element  of  the 
resulting  distribution  function  will  become  negative, 
yielding  an  unphysical  result. 


+R  (v  +  Av  —  Av',t)Al' 


V.  RESULTS  AND  DISCUSSION 


=  /(v  +  Av(l  — l/»),l) 

+  R  (v  + Av(  1  —  1 /«),( )Af'  .  (18) 

A  consequence  of  using  a  time  step  At'  smaller  than 
Al—Avm/(qE )  is  that  the  right-hand  side  of  Eq.  (18) 
calls  for  values  of  /  and  R  from  velocity-space  locations 
which  are  not  explicitly  stored  m  the  array.  We  can  ap¬ 
proximate  Eq.  (18)  into  a  usable  form  with  some  simple 
linear  interpolation  of  the  /  and  R  arrays  from  their 
stored  values  to  obtain 

/(v  +  Av,<  +At') 

=  [(«  —  l)/«][/(v  +  Av,I)  +  f((v  +  Av,r)A('J 


The  algorithm  described  in  Sec.  Ill  was  used  to  obtain 
the  transient  behavior  of  the  EVDF  and  of  various 
electron-swarm  parameters  for  electrons  in  gaseous  neon 
and  gaseous  argon.  The  swarin  parameters  under  investi¬ 
gation  were  Vd,  ( e ),  and  /(,,  which  are  the  drift  velocity, 
average  energy  of  the  electrons,  and  the  ionization  rate  of 
the  gas  atoms,  respectively.  The  time  dependence  of 
these  quantities  was  calculated  from  the  normalized  ve¬ 
locity  distribution  function,  which  was  obtained  at  each 
time  step.  The  normalized  distribution  function  F(v,rp,t) 
is  defined  as  follows: 


f{v,rl’,t) 

r 


(21) 


-Hl/«)[/(v,r)  +  R(v,t)Ar']  .  (19) 

The  above  procedure  for  rigorously  implementing  the 
CFL  stability  condition  is  represented  graphically  in  Fig. 

3. 

Another  condition  for  the  numerical  stability  requires 
that  each  of  the  array  elements  of  the  distribution  func¬ 
tion  /be  larger  than  the  corresponding  array  elements  of 
the  collision  term.  Restating  algebraically, 


where 

A  (t)  —  2-f^°v2du  JJf(v,ifi,l)siiu/idi/i  (22) 

and 

A  (0)=  1  .  (23) 

The  vmious  electron-swarm  parameters  are  defined  as 
follows: 


/(v, t)/[R  (v,f)Af]>  1  .  (20)  /<  —2ttN  f  "v2dv  f’Tucr,m(v)FU>,ifi,t)siiufrdifi ,  (24) 

J  o  Jo 

This  condition  can  be  met  by  selecting  At  such  that  Eq, 

(20)  is  satisfied  for  all  v  Experience  has  shown  that  the  where  <710„Ui)  is  the  ionization  cross  section  of  the  atoms 
ratio  in  Eq.  (20)  must  be  of  an  order  of  102  or  103  for  low  by  electron  impact: 
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Vj  =  2tt  f™v2dv  co&\fiF(v, tp, t )sint/; tl tji  ,  (25) 

and 

(e)=2irfo’eu2dvfJjinu2F(v,ifi,t)S:!v!)di/i .  (26) 

By  using  a  normalized  distribution  function,  the  total 
number  of  electrons  in  the  distribution  function  is  kept 
constant;  this  situation  is  analogous  to  a  steady-state 
Townsend  (SST)  experiment.  A  pulsed  Townsend  (PT) 
experiment  can  be  simulated  by  using  an  unnormalized 
time-evolved  distribution  function,  since  in  such  an  ex¬ 
periment  the  number  of  electrons  docs  not  stay  constant. 
In  a  PT  calculation,  the  final  equilibrium  values  of  the 
swarm  parameters,  obtained  from  the  unnormalized  final 
equilibrium  distribution  function,  will  depend  on  the  ini¬ 
tial  conditions.  On  the  other  hand,  the  equilibrium 
values  of  the  swarm  parameters  in  an  SST  calculation  are 
independent  of  the  initial  conditions.  For  this  reason,  we 
have  chosen  to  present  only  the  normalized  SST  results 
for  various  swarm  parameters.  In  each  case  a  Maxwelli¬ 
an  velocity  distribution  of  electrons  was  assumed  at  time 
t  =0and  a  gas  density  N of  3.54 X  10l6cm"3  (1.32X  10~3 
amagat  or,  equivalently,  1  Torr  at  273  K)  was  used.  Iso¬ 
tropic  scattering  was  assumed  in  all  eases. 

The  convergence  of  the  swarm  parameters  to  their  final 
equilibrium  values  occurs  more  quickly  as  the  value  of 
E/N  becomes  larger.  For  cases  in  which  the  E/N  values 
are  small,  the  time  required  to  reach  equilibrium  can  be¬ 
come  large,  forcing  the  calculation  to  consume  more 
computing  time.  In  cases  of  small  E/N  and  when  only 
the  final  equilibrium  values  of  various  swarm  parameters 
are  of  interest,  the  time-independent  two-term  expansion 
methods  may  be  computationally  more  efficient  (although 
the  time  dependence  will  be  lost).  In  cases  of  large  E/N, 
however,  convergence  is  fast  enough  so  that  the  time- 
dependent  calculations  described  in  this  paper  become 
practical  with  very  modest  computing  resources 

A.  Electrons  in  neon 

For  neon,  the  velocity  steps  Av  ranged  from  2.05  X 107 
to  3.40X  107  cm/sec,  and  the  time  steps  At  ranged  from 
0.1  to  0.02  nsec  as  the  E/N  ratio  was  varied  from  35  to 
566  Td.  The  relevant  scattering  cross  sections  used  in  the 
calculations  were  taken  from  Ref.  2.  Figure  4  displays 
the  calculated  time-dependent  behavior  of  the  electron- 
swarm  parameters  for  E/N  =566  Td  and  with  initial 
average  energies  of  44  and  20  eV.  From  this  figure,  it  is 
evident  that,  although  the  final  equilibrium  values  of  the 
swarm  parameters  are  unaffected  by  the  average  energy 
value  of  the  initial  velocity  distribution  function,  the 
transient  behavior  may  be  considerably  different.  For  ex¬ 
ample,  an  overshoot  in  the  drift  velocity  is  observed  if  the 
initial  average  energy  of  the  EVDF  is  less  than  the  final 
equilibrium  value,  but  the  overshoot  docs  not  appear  if 
the  initial  average  energy  is  somewhat  higher  than  the 
final  value. 

Figure  5  displays  the  time-dependent  behavior  of  the 
electron-swarm  parameters  in  gaseous  neon  with  E/N  ra¬ 
tios  of  35  and  72  Td.  The  initial  value  of  (t)  for  the 
cases  of  35  and  72  Td  are  12  and  16  eV,  respectively.  For 


FIG.  4.  Time  dependence  of  various  elect ron-swarm  parame¬ 
ters  in  gaseous  neon  for  E/N  =  566  Td.  The  data  correspond¬ 
ing  to  the  solid  curves  and  the  dashed  curves  have  initial  aver¬ 
age  electron  energy  of  44  and  20  eV,  respectively 

the  case  in  which  E/N  =35  Td.  an  overshoot  of  the  drift 
velocity  is  present  and  slight  undershoots  of  both  the 
average  energy  and  ionization  rate  can  also  be  seen.  The 
undei shoots  arc  not  seen  in  the  72-Td  data;  however,  a 


TIME  (nsec) 


FIG.  5  Time  dependence  of  various  electron-swarm  parame¬ 
ters  in  gaseous  neon  for  E/N  =  72  Td  (solid  curves)  and  35  Td 
(dashed  curves).  The  upper  and  lower  time  scales  correspond  to 
the  35-  and  72-Td  data,  respectively. 
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FIG  6.  Equilibrium  values  of  various  electron-swarm  param¬ 
eters  in  gaseous  neon.  The  solid  circles,  open  circles,  and  trian¬ 
gles  are,  respectively,  the  equilibrium  values  of  <c>,  R„  and  Vd 
from  the  calculations  of  Kitamori  cl  a!,  (Ref.  2).  The  dashed 
lines  for  the  values  of  E/N  below  35  Td  correspond  to  extrapo¬ 
lated  values  of  various  parameters. 


slight  overshoot  of  the  drift  velocity  is  present.  Whether 
or  not  these  undershoots  or  overshoots  occur  depends  on 
the  initial  value  of  the  average  energy  of  the  swarm  com¬ 
pared  to  that  of  the  final  equilibrium  value.  Recently  ob¬ 
served6,7  current  overshoots  in  Ar-Hg  and  Nc-Hg 
discharges,  and  ionization  rate  overshoots  in  N2 
discharges,  are  presumably  related  to  the  initial  condi¬ 
tions  of  the  electron-energy  distribution  function. 

The  final  equilibrium  values  of  the  electron-swarm  pa¬ 
rameters  as  a  function  of  E/N  are  depicted  in  Fig.  6,  and 
arc  in  very  r  J  agreement  with  the  values  calculated  by 
Kitamori,  Tagashira,  and  Sakai.2  The  zero-field 
(E/N=  0)  values  of  various  parameters  can  be  extrapo¬ 
lated  from  the  curve  of  Fig.  6.  These  values  are 

Krf~»O.Ocm/sec,  R,->0.0  sec-1,  (e)— <-8.7  eV  . 

The  variation  in  slope  of  the  drift  velocity,  particularly 
near  the  lower  E/N  values,  suggests  that  the  electron 
mobility  in  neon,  which  is  related  to  8Krf/3 (.E/N),  is 
shghtly  dependent  upon  E/N. 

B.  Electrons  in  argon 

The  cross  sections  for  the  elastic  scattering  of  electrons 
with  argon  were  taken  from  Massey  and  Burliop,8  and 
the  ionization  cross  sections  were  from  Rapp  and 
Englander-Golden.9  The  total  excitation  cross  sections 


FIG.  7  Time  dependence  of  various  electron-swarm  parame¬ 
ters  in  gaseous  argon  for  E/N  =  565  Td  (solid  curves)  and  283 
Td  (dashed  curves) 

were  adapted  from  Sakai  ei  ul. 10  by  assuming  a  constant 
energy  loss  of  1 1  5  eV  for  all  excitation  processes. 

Figures  7  and  8  display  the  time-dependent  behavior  of 
the  electron-swarm  parameters  for  E/N  ratios  ranging 
from  72  to  565  Td.  For  a  given  value  of  E/N,  the  veloci¬ 
ty  step  An  and  time  step  At  m  the  present  case  are  com¬ 
parable  to  those  of  the  neon  case.  In  each  of  these 
figures,  overshoots  in  the  drift  velocities  are  observed  and 
they  arc  most  dramatic  when  the  initial  average  energy  is 


FIG.  8.  Time  dependence  of  various  electron-swarm  parame¬ 
ters  in  gaseous  argon  for  fc'/IV=72Td. 
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FIG.  9.  Equilibrium  values  of  various  election-swarm  param¬ 
eters  in  gaseous  argon.  The  solid  circles  and  the  open  circles 
are,  respectively,  the  equilibrium  values  of  <e>  and  R,  from  the 
SST  calculations  of  Sakai  a  s'  (Ref.  iO). 

closest  to  the  final  equilibrium  average  energy.  That  the 
size  of  the  overshoot  seems  to  vary  inversely  with  the 
E/N  ratio  is  merely  an  artifact  resulting  from  the  fact 
that  the  initial  average  energy  values  just  happened  to  be 


chosen  closer  to  the  equilibrium  values  as  the  E/N  ratios 
were  reduced.  In  fact,  the  overshoot  can  be  enhanced  or 
completely  eliminated  at  any  E/N  ratio  by  merely  adjust¬ 
ing  the  value  of  the  initial  average  energy. 

Some  preliminary  analysis  of  the  transient  behavior  of 
the  presently  calculated  electron-swarm  parameters  (for 
both  neon  and  argon)  suggests  that  their  time-dependent 
behavior  can  be  accurately  fitted  by  a  sum  of  two  ex¬ 
ponentials.  In  fact,  exponential  behavior  of  the  time 
dependence  of  various  swarm  parameters  can  be  analyti¬ 
cally  justified, 11,12  especially  for  low  values  of  E/N,  by 
assuming  a  constant  collision  frequency  which  leads  to  a 
very  simplified  collision  term. 

The  equilibrium  values  of  the  electron-swarm  parame¬ 
ters  in  argon,  plotted  as  a  function  of  E/N,  are  shown  in 
Fig.  9.  As  was  done  in  the  case  of  neon,  the  zero-field 
values  of  various  swam,  parameters  can  be  obtained  by 
extrapolation  from  this  figure;  this  extrapolation  pro¬ 
cedure  yields  the  following  values: 

Id— >0.0  cm/sec,  R,— >0.0  sec-1,  (e)— >5.2  eV  . 

The  equilibrium  values  in  Fig.  9  can  be  compared  with 
the  values  calculated  by  Sakai  el  a/,10  In  their  paper, 
various  expressions  used  to  define  the  drift  velocity  Vd 
were  different  from  the  expression  used  in  the  piesent  cal¬ 
culations  [Eq  (25)].  Thus  meaningful  comparisons  could 
not  be  made  for  that  parameter.  The  expressions  for  R, 
and  (e)  used  in  the  SST  condition  calculations  of  Ref.  10 
were  equivalent  to  the  expressions  used  in  the  present  pa¬ 
per  so  that  comparisons  among  these  parameters  are 
feasible.  Figure  9  shows  that  the  agreement  between  the 
two  calculations  is  excellent. 

The  unnormahzed  equilibrium  EVDF  in  argon  is 
shown  for  two  different  values  of  E  /N  in  Fig.  10.  In  both 
cases  the  initial  distribution  function  (at  I  =0)  is  a  splieri- 


e'/ ARGON 


(a) 


(b) 


FIG.  10.  Equilibrium  velocity  distribution  function  of  electrons  in  argon  for  (a)  E/N -IS  Td  and  (b)  E/N  =  424  Td. 
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cally  symmetric  Maxwellian  with  (e)=5  eV  for 
E/N  —25  Td  and  <e>  =  30  cV  for  E/N=  424  Td.  At 
equilibrium,  the  distribution  function  retains  much  of  its 
spherical  symmetry  for  low  E/N.  For  large  values  of 
E/N,  however,  the  equilibrium  EVDF  becomes  highly 
asymmetric,  suggesting  that  the  cfiect  of  the  electric  field 
on  the  distribution  function  dominates  over  the  effects  of 
collisions.  This  clearly  indicates  that  the  two-term  ex¬ 
pansion  procedure,  which  retains  only  first-order  devia¬ 
tions  from  spherical  symmetry,  for  obtaining  the  equilib¬ 
rium  EVDF,  would  be  valid  only  for  small  values  of 
E/N  The  present  procedure  for  obtaining  the  EVDF  is 
valid  lor  any  value  of  E  /N,  small  or  large. 

The  “valley"  near  the  origin  (u  =0)  of  the  distribution 
function,  which  becomes  very  pronounced  for  large 
values  of  E  /N,  is  probably  due  to  the  fact  that  very-low- 
energy  electrons  have  a  very  small  collision  probability 
and  are  quickly  accelerated  by  the  electric  field  to  a 
higher  velocity  where  they  become  more  likely  to  have 
collisions.  On  the  other  hand,  the  “upstream”  electrons 
are  not  as  efficiently  accelerated  into  the  origin  (where 
they  would  replace  those  that  have  been  accelerated  out), 
because  their  velocities  are  already  large  enough  so  that 
they  are  inhibited  by  collisions. 

VI.  CONCLUSIONS 

A  very  simple  numerical  algorithm  has  been  described 
which  obtains  the  time-dependent  behavior  of  an 
electron-velocity  distribution  function  in  a  gas.  Aside 
from  its  simplicity,  the  algorithm  has  many  unique  and 
valuable  features.  Unlike  many  other  methods  of  solu¬ 
tion  of  the  Boltzmann  equation,  the  present  method  does 
not  make  use  of  any  term  expansions  of  the  distribution 


function,  and  in  this  respect  the  solution  is  exact.  The 
need  for  numerical  evaluation  of  derivatives  has  been 
completely  eliminated,  which  allows  for  a  much  more 
stable  solution.  The  computational  algorithm  itself  only 
involves  summing  and  shifting  of  various  array  clemt  ts, 
which  can  be  done  without  incurring  any  round-off  error; 
this  fact  enhances  the  stability  of  the  numerical  pro¬ 
cedure.  Because  of  the  simplicity  of  the  procedure,  the 
calculation  can  be  penormed  with  very  modest  comput¬ 
ing  resources.  11ns  is  especially  true  in  cases  of  lngli- 
E/N  values  for  which  convergence  to  equilibrium  is 
much  faster  than  for  low  E  /N. 

Although  the  calculations  that  are  presented  in  this  pa¬ 
per  are  for  electron  swarms  in  a  pure  gas,  like  neon  or  ar¬ 
gon,  subjected  to  a  constant  electric  field,  other  more 
complicated  situations  can  be  very  easily  adapted  to  the 
present  procedure.  For  instance,  the  present  algorithm 
could  be  easily  adapted  to  the  case  in  which  the  electron 
swarm  interacts  with  a  gas  mixture.  Other  situations  of 
interest  include  the  cases  in  which  the  projectile  particles 
are  more  massive  than  electrons  such  as  muons,  protons, 
and  heavy  ions.  Furthermore,  the  present  procedure  can 
easily  accommodate  the  case  in  which  the  external  elec¬ 
tric  field  varies  with  time,  such  as  an  rf  field.  These  other 
applications  of  the  present  procedure  are  under  current 
investigation. 

ACKNOWLEDGMENTS 

It  is  a  pleasure  to  thank  Professor  A.  Garscadden  for 
introducing  us  to  this  research  problem  and  Professor  H. 
II.  Denman  for  valuable  conversaMons.  The  support  of 
the  U  S.  Air  Force  Office  of  Scientific  Research  through 
Grant  No.  AFOSR-87-0342  is  gratefully  acknowledged. 


'A  V  Phelps  and  L  C  Pitchford,  Phys  Rev  A  31,  2932  (1985) 
2K.  Kitamori,  II.  Tagashira,  and  Y  Sakai,  }.  Pliys  D  11,  283 
(1977). 

’P.  J.  Drallos  and  J.  M.  Wadehra,  J  Appl.  phys.  63,  5601 
(1988) 

4T.  Holstein,  Pliys  Rev.  70,  367  (1946). 

5R  Courant,  K  Friedrichs,  and  H.  Lcwy,  Math.  Ann  100,  32 
(1928) 

6M.  E.  Duffy  and  J.  II.  Ingold,  Bull.  Am.  Phys  Soc  34,  307 
(1989). 

7J  T  Verdeyen,  L.  C  Pitchford,  Y  M  Li,  J  B.  Gerardo,  and 
G  N.  Hays,  Bull.  Am.  Phys  Soc.  34,  322  (1989). 

*11.  S.  W  Massey  and  E  11.  S  Burhop ,  Electronic  and  Ionic  Im¬ 


pact  Phenomena,  2-id  ed.  (Clarendon,  Oxford.  England.  1969), 
Vol  1 

9D  Rapp  and  P.  Englander-Golden,  J.  Cliem  Phys  43,  1464 
(1965). 

I0Y  Sakai,  H.  Tagashira,  and  S.  Sakamoto,  J.  Phys  D  10,  1035 
(1977). 

nB.  Shizgal  and  D.  R.  A  McMahon,  Phys.  Rev.  A  32,  3669 
(1985);  D.  R.  A  McMahon,  K.  Ness,  and  B  Shizgal,  J.  Pliys. 
B  19,2759(1986). 

I2I„  C.  Pitchford,  Air  Force  Wright  Aeronautical  Laoralorv 
Technical  Report  No.  AFWAL-TR-85-2106,  1985,  G.  N 
Hays,  L.  C  Pitchford,  J  B.  Gerardo,  J.  T  Verdeyen,  and  Y. 
M.  Li,  Phys  Rev.  A  36,  2031  (1987) 


INTERNATIONAL  JOURNAL  OF  QUANTUM  CHEMISTRY,  VOL  XXXVII.  797-S09  (1990) 


Analytical  Expressions  and  Recursion  Relations  of 
Two-center  Harmonic  Oscillator  Matrix  Elements  of 
Arbitrary  Functions 

P.  J.  DRALLOS  AND  J.  M.  WADEHRA 

Department  of  Pnysics  and  Astronomy,  Wayne  State  University,  Detroit.  Michigan  48202 


Abstract 

The  matrix  elements  of  various  analytical  functions  / (AO,  X  being  the  intcmuclear  separation,  are 
required  for  the  description  of  transition  probabilities  and  other  molecular  properties.  These  matrix  ele¬ 
ments  can  be  conveniently  estimated  by  assuming  vibrational  wave  functions  of  two  relatively  displaced 
linear  harmonic  oscillators  of  arbitrary  frequencies  to  represent  the  vibrational  levels  of  two  electronic 
states  of  a  molecule.  Using  this  assumption,  analytical  expressions  for  the  matrix  elements  of  an  arbitrary 
analytical  function  f(X)  are  obtained.  Useful  recursion  relations  among  these  matrix  elements  are  derived 
and  an  elegant  graphical  representation  of  the  recursion  relations  is  obtained.  These  graphical  representa¬ 
tions  are  utilized  to  obtain  new  more  general  recursion  relations  among  matrix  elements  of  the  arbitrary 
function  f(X). 


I.  Introduction 

In  describing  various  physical  properties,  like  the  oscillator  strengths,  the  transition 
rates  etc.,  of  a  molecule,  the  matrix  elements  of  various  functions  of  the  intemuclear 
separation  X  appear  naturally  in  the  formulation.  A  part  of  these  matrix  elements  is 
normally  to  be  evaluated  between  vibrational  wave  functions  belonging  to  two  differ¬ 
ent  electronic  states  of  the  molecule.  For  small  vibrational  quantum  numbers  it  is  rea¬ 
sonably  accurate  to  replace  the  actual  potential  curves  of  the  relevant  electronic  states 
by  those  of  linear  harmonic  oscillators.  For  the  cases  of  large  vibrational  quantum 
numbers,  where  the  effects  of  anharmonicity  are  relatively  significant,  one  could  use 
a  Morse-type  potential  to  represent  the  potential  curves.  If  one  were  then  to  use  the 
standard  perturbative  techniques,  with  the  linear  harmonic  potential  as  the  zeroth 
order  approximation,  the  harmonic  oscillator  matrix  elements  of  various  functions  of 
X  would  appear  in  the  correction  terms.  In  a  recent  paper  (1)  we  obtained  closed  form 
expressions  and  recursion  relations  among  two-center  harmonic  oscillator  matrix  ele¬ 
ments  of  some  definite  functions  of  X.  The  two  linear  harmonic  oscillators  (lho) 
were  assumed  to  have  arbitrary  frequencies  and  equilibrium  positions.  The  functions 
of  X  considered  earlier  (1)  were  the  exponential,  Gaussian,  and  powers  of  X.  In  the 
present  paper  we  will  generalize  the  earlier  results  to  include  any  arbitrary  analytical 
function  f(X)  which  could  be  expanded  as  a  power  series  in  X. 

The  use  of  linear  harmonic  oscillator  wave  functions  in  obtaining  the  Franck-Condnn 
overlap  integral  dates  back  [2]  to  1930;  however,  there  has  been  some  recent  interest 
(3-5)  in  obtaining  analytical  expressions  as  well  as  recursion  relations  among  the 
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two-center  harmonic  oscillator  integrals  of  arbitrary  functions  of  X  because  of  their 
possible  applications  in  molecular  physics  and  vibrational  spectroscopy.  Two  impor¬ 
tant  conclusions  reached  previously  [))  will  be  of  interest  in  the  present  discussion. 
First,  the  two-center  LIIO  matrix  elements  of  X-dcpcndent  functions  like  the  exponen¬ 
tial  (exp(-aX)]  and  the  Gaussian  (exp(— /3X*)]  could  be  obtained  from  the  corre¬ 
sponding  Franck-Condon  overlap  integrals  by  a  simple  scaling  procedure  Second, 
the  two-center  mo  matrix  elements  of  powers  of  X  (the  Franck-Condon  overlap  inte¬ 
gral  is  a  special  case  of  these  corresponding  to  X°)  satisfy  four-  or  five-term  recursion 
relations  which  can  be  utilized  for  rapid  evaluation  of  various  matrix  elements.  We 
will  show  below  that  the  two-center  lho  matrix  elements  of  an  arbitrary  analytical 
function  f(X)  also  satisfy  some  simple  recursion  relations  and.  furthermore,  these  re¬ 
cursion  relations  can  be  represented  in  an  elegant  graphical  manner  The  graphical 
representation  of  these  generalized  recursion  relations  will  be  used  as  an  aid  in  ob¬ 
taining  new  general  recursion  relations  which  are  also  valid  for  any  analytical  func¬ 
tion  of  X. 


II.  Arbitrary  Function  f(X) 

In  the  present  notation  (X  j  m)  is  the  wave  function  of  the  mth  level  of  the  harmonic 
oscillator  associated  with  the  potential  V,  =  /u  v2/2,  and  «X|n»  is  the  wave  func¬ 
tion  of  the  nth  level  of  the  harmonic  oscillator  associated  with  the  potential 
V,  =  fia i2(X  -  A"0)'/2 .  X0  is  the  separation  between  the  two  oscillators  and.  for  con¬ 
venience,  define  w0  =  h/(txXl).  The  wave  functions  (X  \  m)  and  ((X  \  n»  of  the  two 
harmonic  oscillators  arc  written  :n  terms  of  standard  Hcrmitc  polynomials  Hm  and  //„ 
Note  the  analogy  of  the  present  single  and  double  ket  notation  (for  example.  |»i)  and 
| n»)  with  the  standard  single  and  double  prime  notation  which  is  used  to  distinguish 
between  the  vibrational  eigenfunctions  belonging  to  the  two  different  electronic  states 
of  a  molecule.  For  convenience,  define  «  =  |2<o  ,/(<«,  +  co2)j''2,  b  =  [2(u;/(W|  +  to,))1'2 
and  c  =  ((oi|  +  <u,)/2oj0)  1 :  Using  the  lho  wave  functions,  the  matrix  element  of  the 
arbitrary  function /(X)  can  be  written  as 


(m|/(X)|«»  = 


£4r",r“(i5fi),‘ 


exp 


w,X2 


2w0XoJ 
W;(X  ~  X„)2 


2w0X2 


H  „ 


=  N„ir -,,2f  </« 

x /[*(«)]  cxp[-(« 


hV 

2  i 


\H„(au)HJibUt  -  c)] . 


(1) 


where,  in  the  last  step,  the  integration  variable  has  been  changed  from  X  to  dimen¬ 
sionless  it  as 


«  =  [(«»,  +  ft>2)/(2o>0)f2  •  X/X„  =  c  •  X/X0 
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and 

N„„  =  ((0,(11,)'  *[(<«,  +  exp{-<u,a>]/[2cDo(<i>i  +  <»,]} 

=  [ab/(2m"'m\  n!)]1'1  cxp{— (a/>c/2)2} 

Applying  the  recursion  relation  for  Hermite  polynomials,  Eq.  (A  I),  to  either  Hm(au) 
or//,|fi(«  -  <-)]  in  Eq.  (I)  yields,  respectively, 

2oc(m|X/(X)|n»  =  [2(m  +  l)]w*0<m  +  l|/(,Y)|n» 

+  (2 m)"%(m  -  1|/(X)|«»  (2a) 

and 

2bc{m\(X  -  X0)f(X)\n))  =  [2 (n  +  l)]'%<m|/«)|n  +  1» 

+  (2n)'’%{m\f(X)\n  -  l»  (2b) 

Using,  in  the  integrands  of  Eqs.  (2), 

ii  exp [-(ii  -  y):]  =  -1/2 (d/du)  exp[-(i<  -  y)!]  +  y  cxp[-(u  -  y)2] 

with  y  =  b!c/ 2.  then  performing  integration  by  parts  and  the  necessary  derivatives 
(using  Eq.  (A. 2)),  we  obtain  two  independent  recursion  relations  for  the  matrix  ele¬ 
ments  of  f(X): 

X0(m\df(X)/dX\n))  =  [2(m  -r  l)Hc/«)<m  +  1|/M|«» 

+  (2m)'a(c/a)(\  -  a2)(m  -  1|/(X)|«» 

-  (2n)'abc{m\f(X)\n  -  1»  -  fc2c2<m|/(X)|n»  (3) 

and 

X0(m\df(X)/dX\n))  =  (2(«  +  \)}'n(c/b){m\f(X)\n  +  1» 

+  (2n)'a(c/b){\  -  b2)(m\f(X)\n  -  I)) 

-  (2 m)',2ac{m  -  l|/(X)|n»  +  <i2c2(w|/(X)|n» .  (4) 

Equations  (3)  and  (4)  are  independent  and  can  be  combined  to  obtain  two  additional 
recursion  relations  (which  are  not  independent  from  (3)  and  (4)].  These  four  relations 
are  generalizations  of  the  recursion  relations  for  two-center  lho  Franck-Condon  over¬ 
lap  integrals  obtained  previously  by  Ansbacher  (6)  and  Manneback  (7).  The  relations 
of  Ansbacher  and  Manneback  are  obtained  by  setting /(X)  =  constant  in  (3)  and 
(4)  above. 

An  arbitrary  analytical  function  /(X)  can  always  be  expanded  in  a  Taylor  senes, 
so  that 

m  =  £  (5) 

no  L 

where  fm( 0)  denotes  the  /th  derivative  of  f(X)  evaluated  at  X  =  0.  Operating  on 
Eq.  (5;  with  (»ij  trom  the  left  and  |n»  from  the  right,  a  simple  expression  for  the  ma¬ 
trix  element  of  f(X)  in  terms  of  the  matr'-  dements  of  X‘  is  obtained: 
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*  /'m(0l 

<m|/(X)|«»  =  £  1  „  (m|X'|n))  ■  (6) 

1-0  * 

On  using  the  closed  form  expression  for  the  matrix  elements  of  X1,  given  by  Eq.  (A. 3), 
the  matrix  elements  of  f(X)  can  finally  be  written  as 


•  yl  fn.ll  fi.l-i>1  f  nl  l1 

=  v  f '"to)—  y  y  — — 

yj  '(2c)'  £  ,.o  L<*  "  P)«"  -  9)!j 


{m  -  p\n  -  q)) . 


p\q\(l  -  p  -  q)\ 

The  notation  \a,b  1  refers  to  the  smaller  of  a  and  b.  The  overlap  integral 
(m  -  p\n  -  q))  can  be  easily  obtained  from  Eq.  (A. 4). 

III.  Graphical  Representation 

Two  different  kinds  of  recursion  relations  arc  obtained  above  for  the  matrix 
elements  of  f(X).  One  type  of  recursion  relation  relates  a  single  matrix  element 
(mpf/(,>0|/i»  (m  and  n  refer  to  the  vibrational  quantum  numbers  of  the  two  oscilla¬ 
tors)  to  a  linear  combination  of  matrix  elements  of f(X).  Equations  (2a)  and  (2b)  em¬ 
body  such  recursion  relations.  In  the  other  kind  of  recursion  relation,  which  are  oi 
interest  to  us  here,  a  single  matrix  element  of  df(X)/dX  is  expressed  as  a  linear  com¬ 
bination  of  matrix  elements  of  f(X).  Equations  (3)  and  (4),  which  are  CAamples  of 
this  kind  of  recursion  relation,  can  be  further  generalized  by  writing  the  arbitrary 
function  f(X)  as  the  /tli  derivative  of  another  arbitrary  analytical  function  g(X),  that 
is,/(X)  =  g‘\X): 

Xo</Hs,HI’(X)|n»  =  [2 (m  ¥  1) ]'%:/«)</*  +  l|«<n(X)|«» 

+  (2m)''\c/a)(l  -  a2)  {m  -  l^Wk) 

-  (2n),,2bc(m\g^X)\n  -  !»  -  bV(m\gm(> Oi«.»  (8) 


X0Hg,"l,C>0|/i»  =  [2  (n  +  l)]l/3(c/h)(m!g("(X)|n  +  1» 

+  (2/iHc/f>)(l  -  b2) <m|g(°(20|n  -  1)> 

-  (2m)w«c(m  -  I|gw(2f)|n»+a2cJ(m|gw(A()|n» .  (9) 

Generalized  recursion  relations  (8)  and  (9)  are  shown  graphically  in  Fig.  la  and  b, 
respectively.  Each  solid  circle  in  this  figure  represents  a  matrix  element  (depending 
upon  m,  n,  and  /)  in  the  recursion  relation.  Here,  l  is  the  order  of  derivative  of  g(X) 
with  respect  to  X.  A  similar  graphical  representation  of  the  recursion  relations  for  the 
Franck-Condon  overlap  integrals  [which  are  obtained  by  setting  gm(X)  =  constantj 
was  presented  by  Manneback  [7],  Note  that  in  Eqs.  (8)  and  (9)  matrix  elements  of 
different  orders  of  derivative  of  g(X)  are  mixed  and,  therefore,  m  reality  the  various 
terms  in  these  relations  represent  matrix  elements  of  different  functions  of  X.  A  re¬ 
cursion  relation  that  does  not  mix  matrix  elements  of  different  derivatives  of  a  func- 
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m 

(e) 

Figure  1  The  recursion  relations  of  Eqs  (8), (9),  (II).  (12).  and  (10)  arc  represented 
graphically  in  a,  b,  c.  d.  and  c,  respectively  m  and  n  arc  the  vibrational  quantum  numbers 
of  the  two  linear  harmonic  oscillators,  and  t  is  the  order  of  derivative  of  the  arbitrary  func¬ 
tion  g(X)  with  respect  toX  Each  solid  circle  represents  a  single  matrix  element  appearing  in 
the  recursion  relation. 


non  of  X,  making  it  truly  valid  for  matrix  elements  of  any  analytical  function  of  X,  is 
easily  obtained  by  combining  Fig.  la  and  b  and  eliminating  points  labelled  1  in  the 
figures.  The  resulting  five-term  recursion  relation  is  represented  graphically  in 
Fig.  le,  and  it  corresponds  to  the  following  equation: 

(2)l,2aic(m|g <0(X )|n»  =  b  •  {[(#»  +  1)]%.  +  l|gw(X)|«» 

-  a  •  {[(«  +  DJ'W’WI'.  +  1» 

+  (n)l,2Hg!"(X)|/i  -  1»} . 


(10) 
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In  (his  relation,  both  of  the  indices  (m  and  it)  of  the  matrix  elements  vary.  The  rela¬ 
tion  (10)  is  also  a  generalization  of  the  so-called  "diamond"  relation  [8]  for  one-center 
Ltto  matrix  elements  of  arbitrary  functions.  We  note,  in  passing,  that  relation  (10)  is 
also  obtained  by  combining  Eqs.  (2a)  and  (2b). 

Sometimes  one  needs  (9j  to  have  a  recursion  relation,  such  as  relation  (2),  in 
which  only  one  of  the  indices  (either  m  or  n)  vanes  and  the  other  index  remains  con¬ 
stant.  To  derive  alternate  recursion  relations  of  this  kind,  we  first  obtain  Fig.  lc  and  d 
by  combining  Fig.  la  and  b  and  eliminating  the  points  labelled  2  and  3,  respectively. 
The  two  relations  obtained  in  this  manner  are  merely  alternative  forms  of  the  recur¬ 
sion  relations  (8)  and  (9)  and  do  not  contain  any  information  that  is  not  already  in  (8) 
and  (9).  The  algebraic  expressions  for  these  alternative  recursion  relations  corre¬ 
sponding  to  Fig.  lc  and  d  are,  respectively, 

X0<m|g('*"(X)|n»  =  [2(m  +  l)]'V(m  +  l|g‘"(X)|/i»  -  nJc  W(X)|n» 

+  [2(n  +  l)Hc/f>)(l  -  «:)Hg,0tt)|n  +  1» 

-  (2n)m(c/b)  (m|g(ft(X)|n  -  1»  (11) 

and 

Xo<H8,'H,(X)|"»  =  [2  (n  +  l)]l'Jfcc<m|g(',(X)|n  +  1»  +  *>JcW’(X)|n» 

+  [2(m  +  l)J1,2(c/ a) ( 1  -  b2)(m  +  l|gH’(X)|n» 

-  (2 m)''\c/a){m  -  l|g(,,(X)|n» .  (12) 

Next  the  points  labelled  1,  2.  3,  and  4  of  Fig.  lc  are  eliminated  using  point  2  of 
Fig.  lb  successively.  Figure  2  depicts  this  elimination  procedure.  The  resulting  fig- 


ngurc  2  The  graphical  aid  used  for  deriving  the  recursion  iciauoa  de¬ 

merits  of  the  arbitrary  function  g(X)  The  shaded  portion  encompasses  the  matrix  elements 
that  appear  in  the  denved  recursion  relation 
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urc.  which  is  in  the  form  of  a  triangle  lying  completely  in  the  n-l  plane  with  in  fixed, 
represents  a  nine-term  recursion  relation  in  which  the  matrix  elements  of  different 
orders  of  derivative  arc  still  mixed.  The  resulting  algebraic  expression  for  this  rela¬ 
tion  is 

b :(X9/ c):(»t|g (<* !l|n)>  +  b\X0/c) (2«),'2<m|g',*"|/i  -  l»  - 
fr'(X0/c)[2(n  -r  l)]l,:<m|g'(‘"|n  +  1»  + 

2(1  -  ,r)[n(n  -  l)],'V!g,'l|n  "  2»  -  a4W2n),'2<m|g,',|n  -  l»  + 

{-«4(>V  -  2(2 «  +  1)  +  2 aV(m  +  n  +  l)}(ni|g(',|/i»  - 

(i‘bc[2(n  +  l)]l'2(m|g<,in  +  1)>  + 

2(1  -  «-)[(«  +  2)(n  +  l)]l,2<m|gl'l|n  +  2»  =  0,  (13) 

The  same  procedure  of  eliminating  four  terms  of  Fig.  Id  using  a  single  point  of 
Fig.  la  repeatedly  can  be  used  again  to  obtain  the  following  nine- term  recursion  rela¬ 
tion  with  n  fixed  and  m  and  /  varying: 

+«2(X„/ c)!(m|g"* 2l|n»  +  a\XJc)(2m)'%n  -  !|g'""|n»  - 
a*(XJc)  [2(m  +  l)f2<m  +  l|g"‘"|/t»  + 

2(1  -  -  l)]',2(m  -  2|gll,|n)>  +  ab'c(2m)'%n  -  l|g,,,|n»  + 

{-«Vc2  -  2(2 m  +  1)  +  2 a:b\m  +  n  +  l)}(m|g',l|n»  + 
ab*c[2(m  +  !)]%i  +  l|g,,,j«»  + 

2(1  -  (>2)[(m  +  2) (hi  +  !)]“<«  +  2|gH,|«)>  =  0.  (14) 
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^(m|g,'*"(X)|»»  =  (2  (n  +  l)]l,Jc(m|g(,,(X)|M  +  1» 

-  (2m)'l2c(m  -  l|gw(X)|n»  +  c-'Ms'VOIn)) .  (16b) 

The  recursion  relations  of  Eqs.  (13)  and  (14)  become 

+(Xo/c)JH*,"2,|n»  +  (X0/c)(2«),'J<m|«t't,’l"  -  D>  - 

(. Xjc){2(n  +  l)]l'!(m|g<'*ll|«  +  1»  -  c(2n)''2(m|g(,l|/t  “  1»  + 

{2 (m  -  n)  -  c2}<m|g(»>  -  cf2(«  +  l)]w<»n|j“’|/«  +  1')  =  0  (17a) 

id 

+(V4:<m|g,'tJI|«»  +  (X0/c)  (2m)w(m  -  l|g'"»>  - 
OCJcMm  +  1)}%I  +  l|g<,*,l|n»  +  c(2m)ll2(m  -  l|g«|«»  + 

{2(n  -  m)  -  c!}{m\gm\n))  +  c[2(m  +  1 )] ,,2</m  +  l|gw|fl»  =  0.  (17b) 


Case  2 


When  cu,  is  not  equal  to  <v2  but  X0  =  0  (that  is,  (XJcf  =  2A/{/x(o>,  +  tu2)}  and 
c  =  0],  the  two-center  harmonic  oscillator  matrix  element  of  the  arbitrary  function 
f(X)  becomes, 


M/WI«»  = 

1*0 


(a  \  n  fo.i-p] 

2  m(o»,  +  <o2))  ^  ^ 


2 "'“mlM 


(m  -  p)\  (n  -  q)\ 


p* 0  ?«0 

aW 


J  p! 


■(m  -  pin  -  q)) 


with 


Pm./il 

(m  |  n))  =  Nm  2 
*•0 


m  -  k\.  n  -  k 


W 


•(2  ab)1 


(b2  -  l)(^)(o2  -  l)(=-:‘). 


\k\ 


(18) 


Only  those  terms  in  the  triple  sum  [over  /,  p,  and  q]  for  which  (l  -  p  -  q)  is  an  even 
integer  contribute  to  the  sum.  Also  only  those  terms  in  the  k  sum  for  which  (n  -  k) 
and  (m  -  k)  are  even  integers  contribute  to  the  overlap  integral  (m  \  n)).  The  recur¬ 
sion  relations  of  Eqs.  (8)  and  (9)  reduce  to 

a(Vc)<m|*('H’(X)|n»  =  [2(m  +  DJ'V  +  l|g<0(X)|n» 

+  (2m),l2(\  -  a2)(m  -  l|gl0(X)|«» 

-  (2nyaab{m\gll)(X)\n  -  1»  (19a) 


and 


b(X0lc){m\gu"‘{X)\n))  =  [2(n  +  l)]''W’(X)|/«  +  1)) 

+  (2n)l,2(l  -  b2)(m\gm(X)\n  -  1» 
-  (2 m)',2ab{m  -  l|g"'(X)|«»  . 


(19b) 
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The  recursion  relations  of  Eqs  (13)  and  (14)  reduce  to 

-*■  b HXJc):(m\g +  bHXJcXln)'  3<m|g''‘"|/r  -  1»  - 
/>Wc){2(n  +  l>]w<m|s'""|n  +  1»  + 

2(1  -  <t:)[n(n  -  !)]'3<m|gu'|n  -  2»  + 

{-2(2«  +  I)  +  2a:b'(m  +  n  +  I )}(»i|gl,l|u))  ~ 

2(1  -  «:)((/i  +  2) (/i  +  l)f:<m|g,'’|/i  +  2»  =  0  (20a) 

and 

+ a :(A’0/ c)‘(m|j> 11 ' 2l|/t»  a\XJc)(2m)'2<m  -  l|g“*"|«»  - 

a\XJc)[2(m  +  l)]l,J</n  +  l|g'""|/t»  + 

2(1  -  b:)[m(m  -  l)]l!(m  -  2|g',l|»)>  + 

{-2(2m  +  1)  +  2<rV(m  *■  n  +  l)}(»tjg',l|rt»  + 

2(1  -  b:)[(m  +  2 )(/»  +  l)],,J(m  +  2|s,w|n»  =  0  (20b) 


Case  3 


In  the  case  in  which  w,  =  ca,  =  w  as  well  as  X0  =  0  the  two  wave  functions  (X|/n) 
and  (( X  | «))  belong  to  the  same  LUO,  and  the  single-center  matrix  element  of  the  arbi¬ 
trary  function  f(X)  simplifies  to 


*  /  * l  \  /  2  fm  n  n,  /-pi 

<m|/(X)|«»  =  Z/m(0)  J-  I  1 

ImQ  \4/X<0/  0 

]V1 _ 8, 

<!)' 


I  VI — /  , 

f  2T'mln' 
[(>«  ~  />)!<»  - 


(21) 


where  the  Kronecker  delta  represents  the  overlap  integral  (m  -  p  |  n  -  r/»  Analyti¬ 
cal  expressions  for  single-center  harmonic  oscillator  matrix  elements  for  a  few  defi¬ 
nite  functions  of  X  have  been  obtained  previously  (10).  The  reduction  of  Eq.  (21)  to 
those  expressions  has  provided  a  useful  check  on  the  present  work  The  recursion  re¬ 
lations  (8)  and  (9)  also  simplify  to  the  following  three-term  relations: 


(Vc)<tn|s<'“’(X)|n»  =  [2(m  +  UHm  +  l|g,n(X)|«» 

-  (2»t)'  J</n|g,n(20|n  -  l»  (22a) 

and 

{Xa/c)  (m|gH*"(.V)|«»  =  [2(/t  +  l)]l:<m|g,%Y)|n  +  1» 

-  (2m)'%n  -  l|g'"(X)|n»  (22b) 


The  recursion  relations  (13)  and  (14)  simplify  to  the  following  four-term  relations 

+(Vc)W"2’l«»  +  (AVc-)(2n)l  !(m|g',,"!u  -  0)  - 

(XJc)[2(n  +  l)]l'3(m|g,,*,l|n  +  l»  +  2(m  -  ;i)(m|gl',|n»  =  0  (23a) 


806  DRALLOS  AND  WADEHRA 

and 

+(X0/d:<m|j>l'*:'|n»  +  (Xq/c)  (2m)i  :(m  -  l|g"*"|«»  - 
(X„/c)[2(in  +  I  |J'  '{m  +  l|j>u,"|«»  +  2(-mi  +  n)  (m|j>l"!/i»  =  0.  (23b) 

Case  4 

Consider  the  ease  in  which  all  three  <o,.  oo.  and  X0  are  nonzeio  such  that  to,.  to.  >  a>0 
and  b :  >  Unjin-,).  It  follows  from  Eq.  (7),  on  taking  the  limit  b~  t>  (to0/i o,).  that 

<m|/(X)|/i» - *  W«»  ■/(!>%/ 2) .  (24) 

Note  that  this  scaling  property  is  valid  for  any  arbitrary'  analytical  function  / 

V.  Discussion  and  Conclusions 

In  this  paper  we  have  derived  a  closed  form  expression,  Eq.  (7),  for  the  matrix  ele¬ 
ment  of  an  arbitrary  analytical  function  f{X)  between  the  vibrational  states  belonging 
to  two  distinctly  separated  linear  harmonic  oscillators  of  different  frequencies.  We 
have  explicitly  verified  that  for  the  cases  where  f(X)  is  a  definite  specified  function 
like  an  exponential,  a  Gaussian  or  powers  of  X.  the  expression  (7)  reduces  to  the  pre¬ 
viously  derived  analytical  expressions  ( 1 .  1 1 )  It  essentially  provides  a  useful  check 
on  the  correctness  of  the  present  work.  For  the  case  of  diatomic  molecules,  various 
matrix  elements  obtained  above  by  using  one-dimensional  harmonic  oscillators  could 
be  used  in  a  straightforward  manner.  For  the  cases  of  polyatomic  molecules,  on  the 
other  hand,  one  will  have  more  degrees  of  freedom  and  it  would  be  necessary  to  use 
two  or  more  dimensional  oscillators.  Even  in  such  cases,  as  long  as  the  relevant  func¬ 
tions  /  arc  separable  in  Cartesian  coordinates,  the  matrix  elements  of  /  can  be  evalu¬ 
ated  using  the  techniques  presented  above.  In  particular,  for  a  function  /  which  is 
separable  in  Cartesian  coordinates  x.  v,  and  :.  a  general  matrix  clement  of  /  will  be¬ 
come  a  product  of  three  one-dimensional  matrix  elements,  one  each  for  x,  y,  and 
The  recursion  relations  outlined  above  will,  then,  be  valid  for  any  one  of  these  one¬ 
dimensional  matrix  elements. 

A  number  of  recursion  relations  |Eqs  (2),  (8)-(l4)|  among  the  matrix  elements  of 
the  arbitrary  function  }(X)  and  its  derivatives  arc  presented.  It  should  be  noted  that 
relation  (10)  can  be  used  for  obtaining  the  matrix  clement  (m|/(X)|/i»  for  any  in  and 
n  from  a  knowledge  of  merely  two  numbers  (0|/(A')|0))  and  (0|/(X)|1))  (or 
(0|/(A)|0»  and  (1|/(A')|0))).  Indeed  the  numerical  merit  of  this  procedure  is  that  one 
needs  to  evaluate  only  the  two  lowest-order  (and.  therefore,  the  simplest)  matrix  ele¬ 
ments  using  the  time-consuming  summation  form  (7)  One  can  then  utilize  the  recur¬ 
sion  relations  for  rapid  evaluation  of  any  higher-order  matrix  clement  (m\f  (,Y )!?;)}. 
Furthermore,  the  relations  (2),  (8).  (9).  (II),  and  (12)  can  be  used  judiciously  for 
quickly  obtaining  the  matrix  elements  of  other  functions  of  X  once  the  complete  ma¬ 
trix  (wi|/(X)|«»  has  been  determined  for  a  certain  definite  function  /  The  recursion 
relations  (2),  (13).  and  (14)  are  useful  in  situations  wh<*re  one  needs  a  recursion  rela¬ 
tion  in  which  one  of  the  vibrational  quantum  numbers,  either  m  or  n,  stays  fixed. 
Such  situations  arise,  for  example,  in  investigations  of  vibrational  spectroscopy  |9). 
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As  an  application  of  the  recursion  relations  derived  in  this  paper,  we  have  calcu¬ 
lated  (9|  (lie  resonant  vibrational  excitation  (v,  — »  ly)  of  molecules  Lij  and  N-  by  im¬ 
pact  of  low  energy  electrons.  In  the  resonance  model,  the  projectile  electron  forms  an 
intermediate  autodctaclung  molecular  anion  state  which  decays  into  a  vibrationally 
excited  molecule  after  the  detachment  of  the  electron.  In  the  simple  model  that  we 
used  19].  the  potential  curves  of  the  neutral  molecule  and  of  the  intermediate  resonant 
anion  state  were  replaced  by  those  of  two  one-dimensional  harmonic  oscillators  of  ar¬ 
bitrary  frequencies  and  equilibrium  separations.  The  vibrational  excitation  amplitude, 
then,  involved  matrix  elements  of  the  resonance  width  function  between  vibrational 
wave  functions  of  the  two  oscillators.  Using  the  results  obtained  in  the  present  paper, 
we  obtained  recursion  relations  among  vibrational  excitation  amplitudes  such  that 
when  the  amplitudes  for  the  transitions  (0  — *  0)  and  (0  — >  I)  were  known,  then  exci¬ 
tation  amplitude  for  any  other  transition,  inelastic  or  superelastic,  were  easily  ob¬ 
tained  by  simple  use  of  recursion  relations.  This  simple  model,  involving  the 
recursion  relations  among  the  linear  harmonic  oscillator  matrix  elements,  was  able  to 
account  successfully  the  experimentally  observed  ( 13]  spectacular  peaks  in  the  cross 
sections  for  vibrational  excitation  of  N,  by  low-energy  electron  impact. 

Finally,  we  comment  briefly  on  the  possible  usefulness  of  the  present  results.  As 
mentioned  in  the  Introduction,  the  matrix  elements  of  various  functions  of  inter- 
nuclear  separation  X  are  relevant  in  the  discussions  of  molecular  spectroscopy.  For 
the  cases  of  those  diatomic  molecules  whose  potential  curves  arc  accurately  known, 
these  matrix  elements  could  be  computed  numerically  using  the  actual  potential 
curves.  However,  for  those  molecules  for  which  the  potential  curves  arc  not  accu¬ 
rately  known  —  this  includes  many  diatomic  molecules  and  the  majority  of  poly¬ 
atomic  molecules  —  the  replacement  of  actual  potential  curves  by  those  of  linear 
harmonic  oscillators  or  Morse  oscillators  is  a  reasonable  approximation  (14],  Two- 
center  linear  harmonic  oscillator  matrix  elements  can  then  reasonably  approximate 
the  required  matrix  elements  for  low-lying  vibrational  levels  and  can  provide  impor¬ 
tant  correction  terms  for  high-lying  vibrational  levels.  It  is  appropriate  to  remark  that 
if  one  were  to  obtain  the  anharmonic  correction  terms  for  the  linear  harmonic  oscilla¬ 
tor  by  simply  expanding  the  rotationicss  Morse  potential,  the  integration  coordinate 
in  the  matrix  elements  of  the  correction  terms  would  be  radial  rather  than  Cartesian, 
and  the  limits  of  integration  would  be  from  zero  to  infinity.  In  that  case,  the  matrix 
elements  <m|X'|n)>  would  not  be  explicitly  given  by  Eqs.  (A. 3)  and  (A. 4).  even 
though  the  recursion  relations  of  Eqs.  (2)— (4)  and  (8)— ( 14)  will  still  be  valid.  Use  of 
analytical  harmonic  oscillator  matrix  elements  along  with  the  recursion  relations,  as 
outlined  in  this  paper,  has  the  computational  advantage,  especially  for  polyatomic 
molecules,  over  the  direct  numerical  evaluation  of  these  matrix  elements  in  terms  of 
little  memory  and  time  requirements.  Furthermore,  analytical  results  can  be  applied 
in  ways  that  arc  simply  inaccessible  to  numerical  solutions  since  analytical  results, 
unlike  the  numerical  solutions,  explicitly  show  the  dependence  of  matrix  elements  on 
various  parameters  like  the  mass,  frequency,  etc.,  of  the  oscillator.  As  a  specific 
example,  the  numerical  observation  of  Fraser  1 151  and  of  Nicholls  and  Jarmnin  lift] 
that  the  matrix  elements  (m\X'\n))  for  the  first  positive  system  of  Nj  satisfy,  under 
certain  conditions, 
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{mpf:!»)>  _  (m[X'|»)) 

(rnlX'l^))  </n|X°|n» 

could  be  easily  explained  as  a  special  case  of  the  scaling  law  depicted  analytically  in 
Eq.  (24)  above. 
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We  have  investigated  the  role  played  by  initial  rovibrational  excitation  of  Li;  on  the  cross  sections 
and  rates  for  dissociative  electron  attachment  to  the  molecule.  For  a  given  internal  energy,  ihc  vi¬ 
brational  excitation  enhances  the  attachment  cross  section  more  than  the  rotational  excitation.  The 
attachment  cross  sections  and  the  attachment  rates  reach  their  maximum  values  when  the  process 
of  dissociative  attachment  to  rovibrationally  excited  molecules  is  still  cndoergic  and,  furtheimore, 
these  quantities  stay  close  to  their  maximum  values  even  when  the  process  changes  from  being  en- 
doergic  to  exocrgic.  The  tipper  bounds  on  the  cross  sections  and  the  rates  for  dissociative  electron 
attachment  to  Li2  are  12.8  A  and  1.25X  10' 8  cm's1.  At  a  fixed  electron  temperature,  the  l.metic 
energy  of  the  negative  ion  formed  by  this  process  increases  as  the  vibrational  quantum  number  of 
the  initial  neutral  molecule  increases;  the  maximum  kinetic  energy  of  the  Li  ~  ion  formed  by  attach¬ 
ment  to  the  u  =  1 2  level  of  Li 2  is  0. 1 53  eV 


I.  INTRODUCTION 

It  has  been  amply  demonstrated,  both  theoretically1  as 
well  as  experimentally,2  that  the  cross  sections  for  II- 
formation  via  the  process  of  dissociative  electron  attach¬ 
ment  to  molecular  hydrogen  are  significantly  enhanced  if 
the  molecule  H2  is  initially  rovibrationally  excited.  In  or¬ 
der  to  assess  the  analogous  degree  of  enhancement  for 
other  molecular  systems,  we  have  presently  investigated 
in  detail  the  efTect  of  initial  rovibrational  excitation  on 
the  rate  of  production  of  Li ~  via  the  process  of  dissocia¬ 
tive  electron  attachment  to  diatomic  lithium  molecules. 
Investigations  of  electron  attachment  to  lithium  mole¬ 
cules  are  especially  appropriate  at  this  time  due  to  several 
reasons.  First,  recent  experimental  observations’  reveal 
that  the  rate  of  Li-  formation  by  the  impact  of  thermal 
electrons  on  highly  vibrationally  excited  Li2  is 
(2d  I ) X  10s  cm's  '.  One  of  the  principal  aims  of  the 
present  work  is  to  confirm  this  expet internal  observation 
of  attachment  rate  by  explicit  calculations  as  well  as  to 
provide  detailed  cross  sections  for  electron  attachment  to 
Li2.  Second,  the  Li"  ions  could  possibly  play,  m  the  fu¬ 
ture,  the  same  roles  as  have  been  played  by  II-  ions  for 
neutral  beam  formation.  Third,  since  lithium  dimers  are 
isovalent  with  H2,  an  investigation  of  the  dependence  of 
the  electron  attachment  to  Li2  on  the  initial  rovibrational 
excitation  of  the  molecule  would  be  similar  to  the  previ¬ 
ous  detailed  study4  on  H2. 

Similarities  between  Li2  and  II2  suggest  that  theoretical 
methods  used  successfully  in  the  past  for  obtaining  the 
cross  sections  and  rates  of  electron  attachment  to  H2  can 
be  employed  for  analogous  investigations  for  Li2.  In  par¬ 
ticular,  the  process  of  dissociative  electron  attachment  to 
Li2  is  understood  to  proceed  through  the  formation  of  an 
intermediate  resonant  anion  state  Li2  which,  on  dissoci¬ 
ation,  leads  to  Li-.  The  fact  that  both  the  lithium  dimer 
molecules  and  the  hydrogen  molecules  are  isovalent  pro- 
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vides  similarities  between  *!>e  electronic  configurations  of 
the  two  molecules.  For  example,  the  lowest  electronic 
states  of  the  negative  molecular  ions  with  configurations 
( los  )2(  la„  )2(2ag  )2(2ou)  for  Li2-  and  (\ag  )2(  lou)  for 
II,  have  similar  symmetry,  namely,  :2„‘ .  However, 
compared  to  the  hydrogen  molecule,  the  lithium  mole¬ 
cule  possesses  a  large  polarizability  and  a  weak  bond 
strength  winch  makes  the  ground  state  of  Li2-  a  true 
bound  state.  In  the  ease  of  H2-,  on  the  other  hand,  the 
22  4  slate  is  a  true  bound  state  only  for  mtermiclear  sepa¬ 
rations  R  larger  than  2.9  a.u  and  is  an  autodetaching 
state  for  smaller  values  of  R.  The  first  excited  state  of  the 
negative  molecular  ions  with  symmetry  2lg  and 
configurations  ( ]<r? )2(  ltru  )2(2</g  )(2ou  )2  for  Li2‘  and 
( \<xg )( Io„  )2  for  II,  is  partly  Feshbach  and  partly  shape 
resonance  in  nature  for  both.  This  state  is  the  essential 
channel  for  dissociative  attachment  of  low-energy  elec¬ 
trons  to  lithium  molecules.  The  process  of  dissociative 
electron  attachment  to  molecular  lithium,  then,  is: 

e  -f-LijLY'Z4)— Li,  (A  %'  )-->Li  +  Li-  . 

II.  CALCULATIONS 

Fortunately,  a  number  of  accurate  calculations5  9  of 
the  potential  curves  of  the  X  state  of  Lt2  and  the 
A  2Z4  state  of  Li,  are  available.  The  potential  curves 
that  we  utilize  in  the  present  work  were  obtained1’  by  an 
ab  initio  calculation  using  optimized  configuration- 
interaction  (Cl)  wave  functions  built  from  orthonormal 
Slater-type  orbitals  (STO's).  Because  of  its  nature  (name¬ 
ly  Feshbach)  the  2£4  resonance  of  Li,  is  expected  to 
have  a  small  width  and  a  long  lifetime.  Potential  curves 
of  the  electronic  states  of  Li-  and  I  i-  relevant  to  the  at¬ 
tachment  process  are  shown  in  Fig  1 .  The  A  21*  elec¬ 
tronic  stale  of  Li2  exhibits,  due  to  its  autodetaching  na¬ 
ture.  a  complex  potential-energy  curve  whose  real  part 
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Intemuclear  Separation  (a  u.)  Intcrnutlcar  Separation  (a  u  ) 


FIG.  1.  The  potential  curves  of  U2(X  '2/1  and  Li2‘  ( A  :2^).  FIG.  2.  The  width  of  the  A  *2/  state  of  Li2~ 


F'(R),  along  with  the  potential  curve  of  the  X  '%+  state 
of  Li2,  is  shown.  The  potential  minima  of  the  X  and  the 
A  curves  are  at  5.05  and  5.91  a.u.,  respectively.  The  two 
curves  cross  at  R  =  R,  =  6.51  a.u.  and,  therefore,  the  pos¬ 
sibility  of  autodctachment  of  the  resonant  molecular 
anion  exists  only  for  intemuclear  separations  smaller 
than  R,,  the  so-called  stabilization  radius.  The  width  for, 
equivalently,  the  lifetime)  HR)  of  this  resonance  is  relat¬ 
ed  to  the  imaginary  part  of  the  complex  potential-energy 
curve  of  the  A  22*  state  of  Li2~.  In  the  present  autode¬ 
tachment  process  22st'-*l2g!  +e~,  the  lowest  contribut¬ 
ing  partial  wave  is  an  s  wave  and,  thus,  by  Wigncr’s 
threshold  law10  the  width  of  this  state  is  given  by 
r (R)  =  c*(R),  where  k(R )  is  the  wave  number  of  the 
electron  emitted  at  intemuclear  separation  R  and  c  is  a 
constant.  Tilts  constant  c  is  obtained  as  follows.  The  au¬ 
todetachment  width  T  is  related,  by  Fermi's  golden  rule, 
to  the  matrix  element  V(R)  coupling  the  discicte  reso¬ 
nant  state  with  the  continuum  state  of  the  electron- 
molecule  system.11  The  wave  function  of  the  discrete 
A  state  was  estimated  by  smoothly  extrapolating  the 
fully  optimized  exponents  of  the  Cl  wave  functions  front 
the  vanationally  stable  region  (R  SR,)  into  the  autode- 
tachmg  region  (R  <R,).  The  continuum  state  was  ap¬ 
proximated  by  extrapolation  of  a  series  of  wave  functions 
constructed  by  adding  an  electron,  in  a  series  of  diffuse 
STO’s,  to  the  unperturbed  ground-state  wave  function  for 
Li2(X  'Ig).  The  coupling  matrix  clement  K(R )  was  then 
Calculated  by  using  these  two  wave  functions.  A  compar¬ 
ison  of  the  width  calculated  by  using  the  golden  rule  with 
the  threshold  law  expression,  given  above,  yielded 
c=0.0i43  a.u.  The  width  of  the  A  2Z*  state  of  Li2“  as 
a  function  of  the  intemuclear  separation  R  is  shown  in 
Fig.  2  and  is  given,  in  atomic  units,  by  f(R) 


=  0.0143*  (R). 

In  the  traditional  resonance  theory,  within  the  local 
approximation,  the  radial  nuclear  wave  function  £(R)  of 
the  resonant  anion  state  satisfies12 


fi2  £ i1  n2j,ut  + 1) 

-r  + - — - 1 -VAR) 

2 M  dR 2  2 MR2 


—  yir(R)— £ p'(R)=  —  V(R)t,v  J(R)  ,  (1) 

where  f„_,(R)  is  the  nuclear  wave  function  of  the  initial 
rovibrational  state  of  the  neutral  molecule  and  M  is  the 
reduced  mass  of  the  nuclei.  E  is  the  total  energy  of  the 
system  and  its  conservation  in  the  initial  channel  (elec¬ 
tron  plus  molecule)  and  the  final  channel  (atom  plus  nega¬ 
tive  ion)  provides  an  expression  foi  the  threshold  for  dis¬ 
sociative-  election  attachment  (DA)  to  the  molecule  that 
is  rovibrationally  excited  to  a  particular  level, 


£t>A_j^  jV,:a  t’lJ  ^r.A- 

lh  [0,  otherwise. 


Here  D  is  the  dissociation  energy  of  the  lithium  molecule 
(1.0372  eV),  is  the  electron  affinity  of  the  lithium 
atom  (0.6182  eV)  and  EvJ  is  the  excitation  energy  of  the 
initial  rovibrational  level  of  the  molecule. 

The  nuclear  wave  function  f(R)  is  obtained  by  numeri¬ 
cally  solving  Eq.  (1)  subject  to  the  boundary  conditions 


i(R  =0)=0  , 

§(R  — > oo )  -  ►RR/ijll(RR )  . 

Here  ii1K2/2M  is  the  relative  kinetic  energy  of  the  lon- 
atom  pair  after  dissociative  electron  attachment  and  hjU 
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is  the  spherical  Hankel  function  of  the  first  kind. 

If,  after  the  formation  of  the  temporary  resonant  state, 
the  nuclei  separate  to  R  >  Rs  without  autodetachment 
having  occurred,  the  detachment  of  the  election  becomes 
energetically  impossible  and  dissociative  attachment  can 
result.  Thus  for  R  <RS,  the  envelope  of  l|(f?)l2  de¬ 
creases  with  R  because  of  the  possibility  of  autodetach¬ 
ment  and  it  is  the  asymptotic  (R  —  oo)  value  of  |£(,R)|2 
which  determines  the  probability  of  dissociative  electron 
attachment  to  the  molecule.  In  fact,  the  integrated  cross 
section  for  dissociative  attachment  of  an  electron  with 
energy  ftik?/2m  to  the  molecule  is  given  by12  (assuming 
momentum-normalized  functions) 

<7da=4it2-“  lim  \$(R)\2  .  0) 

M  K,  R  -•  oo 


In  order  to  convert  the  attachment  cross  sections  into 
attachment  rates  one  needs  the  energy  distribution  of 
electrons.  In  the  present  work  we  take  this  distribution 
to  be  Maxwellian,  namely, 

27 E  ”/2 


/(£)= 


2-rrE 


exp 


3 E 
2  E 


(4) 


where  the  aveiage  energy  E  is  related  to  the  electron  tem¬ 
perature  7'via  E=2kBT/2.  The  attachment  rate  ft (C)  is 
merely  a  convolution  of  (2£/m)1/2<7DA(£),  that  is, 
r  U/2 


k(E)= 


111 


J*Eu2crDA(E)f(E)dE  . 


(5) 


III.  RESULTS  AND  DISCUSSION 

In  the  present  work  we  have  calculated  the  cross  sec¬ 
tions  and  the  rates  for  dissociative  electron  attachment  to 
molecular  lithium  when  the  molecule  is  either  in  one  of 
the  vibrational  levels  u  —0  to  12,  all  rotationless,  or  m 
one  of  the  rotational  levels  J  =0  to  25,  all  with  v  =0. 
Figure  3  shows  the  cross  sections,  as  a  function  of  the  in¬ 
cident  electron  energy,  for  election  attachment  to  Li2  in 
various  iv,J)  levels.  Besides  the  lowest  level  (0,0),  the 
other  two  levels  shown  in  Fig.  3  (1,0)  and  (0,22)  have  ap¬ 
proximately  the  same  internal  energy.  The  attachment 
cross  sections  exhibit  a  rapid  increase  at  the  threshold, 
attaining  a  peak  value  <7^,  followed  by  a  uniform  de¬ 
crease  as  the  energy  of  the  incident  electron  increases. 
This  almost  vertical  onset  of  the  cross  section  is  attribut¬ 
ed  to  the  attractive  nature  of  the  potential  curve  of  the 
resonant  anion  state.13  The  attachment  rate  is  essentially 
determined  by  the  peak  attachment  cross  section.  The 
results  in  Fig.  3  clearly  show  that  the  attachment  cross 
section  is  enhanced  if  the  molecular  lithium  is  initially  ro- 
vibrationally  excited.  A  part  of  the  enhancement  of  the 
cross  section  occurs  due  to  the  lowering  of  the  threshold 
for  attachment  as  the  molecule  is  iuvibi.uiuu.iiiy  exciied. 
Furthermore,  for  a  fixed  internal  energy,  vibrational  exci¬ 
tation  of  the  molecule  is  more  effective  in  enhancing  the 
cross  section  than  rotational  excitation.  For  example, 
when  the  molecule  is  initially  provided  with  an  internal 
energy  of  0.04  eV,  then  the  peak  attachment  cross  section 


FIG.  3  The  cross  sections  for  dissociative  attachment  of 
low-energy  electrons  to  molecular  lithium  in  various  rovibra- 
tional  tv,J)  levels 


is  increased  by  a  nctor  of  6  8  if  this  internal  energy  is 
purely  vibrational  (from  v  =0  to  1)  and  by  a  factor  of  2.5 
if  the  internal  energy  is  purely  rotational  (from  J  =0  to 
22).  Tabic  I  provides  the  energetics  as  well  as  the  peak 
cross  sections  for  dissociative  electron  attachment  to  the 
ground  electronic  state  of  molecular  lithium  in  various 


TABLE  I  Internal  rovibrationa!  energy  threshold  for 
dissociative  attachment  (E [J1'')  ;">d  peak  attachment  cross  sec¬ 
tion  (o^)  for  various  rovibrational  levels  of  the  ground  elec¬ 
tronic state  ofLi]  Asterisk  denotes  exoergic. _ 


V 

J 

Ej  (eV) 

AT  (eV) 

Cfpcak  <A2) 

0 

0 

00 

0  4190 

0.368 

i 

0 

0  04292 

0.3761 

2  50 

2 

0 

008514 

0  3339 

3  96 

3 

0 

0.1267 

0.2923 

3.69 

4 

0 

0.1676 

0.2514 

4  52 

5 

0 

0  2078 

0.2112 

4.02 

6 

0 

0.2472 

0.1718 

9.26 

7 

0 

0.2860 

0.1330 

12.8 

8 

0 

0  3240 

00950 

128 

<> 

0 

0.3612 

0.0578 

8.94 

0 

0  3976 

0.0214 

10  3 

II 

0 

O.S333 

* 

9.43 

0 

1 

000017 

04188 

0.369 

0 

1 

0.00030 

0.4 » 55 

0.312 

0 

5 

0.00249 

04165 

0.389 

0 

10 

0.00910 

0.4099 

0.450 

0 

20 

0.03461 

0,3844 

0  809 

0 

22 

0.04164 

0  3774 

0.930 

0 

25 

005338 

0.3656 

1.16 

3610 
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rovibrational  levels.  It  is  interesting  that,  unlike  molecu¬ 
lar  hydrogen,  the  peak  attachment  cross  section  here 
does  not  increase  uniformly  as  the  internal  vibrational  en¬ 
ergy  is  increased.  Whether  this  is  due  to  the  local  nature 
of  our  present  calculations  or  to  the  nature  of  the  poten¬ 
tial  curves  of  the  lithium  molecule  remains  to  be  ascer¬ 
tained  in  the  future.  However,  to  maintain  numerical 
consistency  of  our  calculations  the  following  equality  was 
satisfied,  as  a  check,  to  within  a  few  parts  in  104  for  all 
values  of  the  incident  electron  energy  and  for  all  rovibra¬ 
tional  levels  of  the  molecule, 

hrn  \$(R)\2+  J*r(RmR)\2dR 

=  2/o"lm[£*m>fl,Ji(fl)]rfR  .  (6) 

This  relationship  is  obtained  by  multiplying  Eq.  (1)  by 
g*( /? ),  subtracting  the  resulting  equation  from  its  com¬ 
plex  conjugate  and  then  integrating  over  R  from  0  to  co. 

Unlike  the  case  of  molecular  hydrogen,  the  width  of 
the  resonance  responsible  for  attachment  to  Li2  is  quite 
small  (or,  equivalently,  the  lifetime  is  large)  so  that  the 
resonance  model  is  very  reasonable.  The  enhancement  of 
the  attachment  cross  section  is  directly  related  to  the  in¬ 
crement  of  the  range  of  intcrnuclcar  separations,  due  to 
internal  rovibrational  excitation,  over  which  the  electron 
can  be  captured.  This  range  is  increased  due  to  an  in¬ 
creased  vibrational  amplitude  during  vibrational  excita¬ 
tion  and  centrifugal  stretching  during  rotational  excita¬ 
tion.  The  probability  of  an  electron  capture,  to  form  the 
resonant  anion  state,  is  maximum  at  an  internuclear  sepa¬ 
ration  at  which  the  energy  difference  between  the  poten¬ 
tial  curves  of  Li2  and  Li2 '  is  equal  to  the  energy  of  the  in¬ 
cident  electron.  This  internuclear  separation  is  referred 
to  the  capture  radius,  Rc.  As  the  nuclei  separate  from  Rc 
to  R,  the  autodetachmcnt  of  the  electron  from  the  anion 
state,  leaving  behind  a  rovibrationally  excited  neutral 
molecule,  is  a  distinct  possibility.  Dissociative  attach¬ 
ment,  of  course,  results  when  the  internuclear  separation 
far  exceeds  the  stabilization  radius  Rs.  For  Li2  molecules 
in  levels  v  >  7  with  J  =  0  the  range  of  intermolecular  sep¬ 
arations  over  which  the  electron  could  be  captured  in¬ 
cludes  Rs.  The  energetic  threshold  for  dissociative  elec¬ 
tron  attachment  to  Li2  is  nonzero  for  vibrational  levels 
ti  5  10.  If  the  molecule  is  initially  in  vibrational  levels 
v  £  1 1,  the  attachment  process  is  exoergic. 

Figure  4  shows  the  effect  of  initial  vibrational  versus 
rotational  excitation  of  the  molecule  on  the  rate  of  elec¬ 
tron  attachment  as  a  function  of  electron  temperature. 
As  seen  in  this  figure,  the  initial  vibrational  excitation  to 
the  t>  ~  1  ( J  =0)  level  enhances  the  attachment  rale  more 
dramatically  than  the  initial  rotational  excitation  to  the 
J  =22  (v  =0)  level  which  has  roughly  the  same  internal 
energy.  The  rate  of  electron  attachment  to  l,i2  in  various 
vibrational  levels  is  shown  in  Fig.  5.  We  note  the  in¬ 
teresting  feature  that  the  attachment  rate  saturates  at  the 
value  1.25 X 10— 8  cm3s~'  which  is  reached  when  the 
molecule  Li2  is  initially  in  vibrational  levels  u  — 8-12. 
Furthermore,  this  rate  is  achieved  for  electrons  with  an 
average  energy  of  —0.2  eV.  These  results  are  consistent 


ii  -  (3  knTc/2)  (eV) 


FIG.  4  Rates  of  dissociative  electron  attachment  to  Li,  in 
various  rovibrational  U \J)  levels  as  a  function  of  the  electron 
temperature.  The  levels  ll.O)  and  10,22)  have  approximately  the 
same  excitation  energy 

with  the  recent  experimental  observations’  of  the  rate  of 
attachment  of  thermal  electrons  to  molecular  lithium. 

A  quantity  which  often  is  of  interest  to  experimental¬ 
ists  and  is  useful  for  plasma  diagnostic  purposes  is  tile  en- 


0  1  2  3  4  S 

£-(UbTc/2)  <eV| 


FIG.  5.  Rates  of  dissociative  electron  attachment  to  Li2  in 
various  vibrational  levels  as  a  function  of  the  electron  tempera¬ 
ture. 
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crgy  of  the  negative  ion  formed  during  the  process  of  dis¬ 
sociative  electron  attachment.  The  energy  of  this  ion  de¬ 
pends  on  the  energy  of  the  incident  electron  E,  as  well  as 
on  the  excitation  energy  E„j  of  the  initial  rovibrational 
level  of  the  molecule.  In  case  of  a  homonuclear  diatomic 
molecule  like  Li2  the  atomic  anion  carries  one-half  of  the 
relative  kinetic  energy  of  the  ion-atom  pair.  The  kinetic 
energy  of  the  anion  is 

£-=i(£,-D+A'EA-H^).  (7) 

Assuming  a  Maxwellian  distribution  of  energies  for  the 
electrons,  the  average  kinetic  energy  of  the  negative  ions 
is  given  by 

<E-)  =  J”E-oDA(E)f(E)dE/ J\DA(E)f(E)dE  . 

(8) 

Figure  6  shows,  as  a  function  of  the  electron  temperature, 
the  average  energy  of  the  Li~  ions  formed  by  dissociative 
electron  attachment  to  Li2  in  various  vibrational  levels. 
Since  the  attachment  process  becomes  exoergic  when  the 
molecule  is  initially  in  vibrational  levels  v  >  1 1,  the  rela¬ 
tive  kinetic  energy  of  the  ion-atom  pair,  then,  is  more 
than  the  energy  of  the  incident  electron  and  therefore  the 
average  energy  of  Li'  ions  formed  by  dissociative  elec¬ 
tron  attachment  to  molecular  lithium  continues  to  in¬ 
crease  as  the  internal  vibrational  energy  increases. 
Among  the  vibrational  levels  that  we  have  investigated 
(tr  — 0— 12)  the  maximum  kinetic  energy  of  the  Li '  ions  is 
0  153  eV. 

In  an  actual  plasma,  though,  the  energy  distribution  of 
electrons  is  non-Maxwcllian  and  only  a  detailed  solution 
of  the  Boltzmann  equation1,1  including  all  possible  col¬ 
lision  mechanisms  of  electrons  provides  the  real  distribu¬ 
tion  If  available,  this  real  energy  distribution  of  elec¬ 
trons  should  be  used  for  calculating  the  attachment  rates 
and  the  average  energy  of  the  negative  ion. 

We  conclude  by  noting  that  the  limited  experimental 
information3  on  the  rates  of  dissociative  electron  attach¬ 
ment  to  molecular  lithium  that  is  cuirently  available  is 
consistent  with  our  theoretical  results.  Detailed  experi¬ 
ments  investigating  the  effect  of  initial  rovibrational  cxci- 


0  I  2*  3  4  5 

P  -  (3  XBTCZ2)  (eV) 


FIG  6  Average  energy  of  Li  ions  formed  by  dissociative 
electron  attachment  to  Li2  in  various  vibiational  levels  as  a 
function  of  electron  temperature. 


tation  of  Li2  on  the  cross  sections  foi  formation  of  L> 
would  be  worthwhile  m  the  future. 
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ABSTRACT 

Formation  of  beams  of  negative  ions  by  the  process  of  dissociative  electron  attachment  to  a  molecule 
AB  (e"  +  AB  —  A  +  B")  is  investigated.  It  is  demonstrated  that  the  initial  rovibrational  excitation  of  the 
molecule  AB  plays  a  significant  role  in  enhancing  the  rate  of  production  of  negative  ions.  A  simple 
t  physical  picture,  involving  the  formation  of  a  transient  resonant  state  of  the  molecular  anion,  is  used  to 
explain  this  enhancement. 


1 . . ,  INTRODUCTION 

Production  of  high  intensity  neutrai  particle  beams  by  neutralization  of  negative  ion  beams  is  now  well 

established1.  The  negative  ion  beams  are  accelerated  to  high  energies  and  then  neutralized  by  electron 
detachment.  One  efficient  way  of  producing  beams  of  negative  ions  is  by  the  process  of  dissociative 
electron  attachment  to  molecules.  Schematically,  this  process  for  attachment  to  a  molecule  AB  is:  e'  +  AB 

—  A  +  B*.  It  has  been  extensively  demonstrated,  both  experimentally2  as  well  as  theoretically3,  that  the 
cross  sections  and  the  rates  of  production  of  negative  ions  by  this  process  are  dependent  upon  the  initial 
rovibrational  level  of  the  attaching  molecule.  In  the  cases  of  molecular  hydrogen  and  molecular  lithium, 
the  rates  of  production  of  negative  ions  H'  (as  well  as  D'  and  T")  and  Li"  by  dissociative  electron 
attachment  to  FU  (and  its  five  isotopes)  and  Lfy  respectively,  are  known  to  be  strongly  dependent  upon 
the  amount  of  internal  rovibrational  energy  of  the  neutral  molecule.  For  example,  the  rates  for  negative  ion 

beam  production  via  electron  attachment  to  molecular  hydrogen  (and  its  isotopes)  are  enhanced4  by  orders 
of  magnitude  if  H2  is  initially  vibrationally  pumped.  In  the  present  paper  only  the  results  of  our 
investigations  on  the  effect  of  initial  rovibrational  excitation  on  the  cross  sections  and  rates  of  dissociative 
electron  attachment  to  molecular  hydrogen  and  its  five  isotopes  will  be  discussed. 

2.  THE  RESONANCE  MODEL 

A  traditional  way  of  describing  the  physics  of  the  process  of  dissociative  electron  attachment  to  a 
molecule  AB  is5  via  the  formation  of  a  temporary  compound  state  of  the  electron-molecule  system.  In  this 
molecular  anion  state  AB*  (also  called  the  resonance  state)  the  electron  can  autodetach  with  a  finite  lifetime 
(related  to  the  width,  I',  of  the  resonance),  leaving  behind  a  vibrationally  excited  neutral  molecule.  On  the 
other  hand,  if  the  iifetime  of  the  resonance  is  long  enough,  the  anion  AB"  can  dissociate  into  A  +  B", 
corresponding  to  the  process  of  dissociative  electron  attachment. 


A  possible  scenario  of  the  resonance  model  is  depicted  in  Figure  1.  Shown  schematically  in  this 
Figure  are  the  potential  curves  of  the  neutral  molecule  AB  (labelled  V0)  and  of  the  resonant  state  AB" 
(labelled  V").  The  two  potential  curves  cross  at  an  intemuclear  separation  R  =  Rs  such  that,  for  R  2  Rs, 

Q  t  Ifrvjotnr'hmonr  />f  »Ua  « n  AvtsJ  *hn  mkOAMOMAA  ‘i 

-  »-kvMV»MViu»ivii4  Vi  U1V  V1VVUV1I  1J  Ulll>lgVUWUUJ  IIV/l  pVlllUUWU  U11U  U1V  IVJVIHUIW  iUtllJ  4 

state.  Rs  is  referred  to  as  the  stabilization  radius.  Before  electron  capture,  the  nuclei  are  rovibrating 
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INTERNUCLEAR  SEPARATION  R 


Figure  1.  A  schematic  representation  of  the  processes  of  dissociative  electron  attachment  and  vibrational 
excitation  of  a  molecule  AB  via  resonance  formation. 

in  the  level  (vj.Jj)  under  the  influence  of  the  potential  V0(R).  After  electron  capture,  the  nuclei  of  the  anion 

move  under  the  influence  of  V’(R).  The  probability  of  electron  capture  to  form  the  resonant  molecular 
anion  state  depends  on  the  intemuciear  separation  and  this  probability  is  maximum  at  an  intemuciear 
separation  (labelled  Rc  in  the  Figure  and  referred  to  as  the  capture  radius)  at  which  the  energy  separation 
between  the  two  potential  curves  is  equal  to  the  energy  of  the  incident  electron.  If  the  potential  curve  V'  is 
repulsive  in  nature,  the  nuclei  in  the  anion  state  begin  to  separate  such  that  the  electronic  potential  energy  is 
convened  into  nuclear  kinetic  energy.  Now  if  the  autodetachment  of  the  electron  occurs  at  some  specific 
intemuciear  separation,  labelled  R  in  the  Figure,  the  neutral  molecule  is  left  in  a  rovibrationally  excited 
le>  cl  due  to  the  gain  in  the  nuclear  kinetic  energy  (indicated  by  a  vertical  dotted  line  in  the  Figure).  The 
exact  rovibrationally  excited  level  (vf,Jf)  achieved  by  the  molecule  depends  on  the  gain  in  the  kinetic 
energy  of  the  nuclei  as  well  as  on  the  relevant  selection  rules.  Depending  upon  the  lifetime  of  the 
resonance  the  nuclei  in  the  anion  state  may  separate  to  an  intemuciear  separation  larger  than  Rs  beyond 
which  the  autodetachment  of  the  electron  is  energetically  not  possible  and  dissociative  attachment  may 
occur  resulting  in  the  formation  of  a  stable  negative  ion. 

3.  RESULTS  AND  DISCUSSION 

The  relative  masses  of  the  six  isotopes  H2,  HD,  HT,  D2,  DT  and  T2  are  1.00,  1.33,  1.50, 2.00, 2.40 
and  3.00,  respectively.  We  have  done  fully  quantum  mechanical  calculations  of  cross  sections  and  rates  of 

dissociative  electron  attachment  to  these  six  isotopes.  In  particular,  we  have  investigated6  the  effect  of 
initial  rovibrational  excitation  in  enhancing  the  rate  of  negative  ion  production.  The  qualitative  behavior  of 
the  attachment  cross  sections  is  the  same  for  all  isotopes ,  namely,  a  sharp  threshold  peak  followed  by  a 
rapid  decrease  on  increasing  the  impact  electron  energy.  The  numericai  resuits  are  provided  111 1  able  1 
which  summarizes  three  enhancement  factors: 

(a)  For  a  given  rovibrational  level,  the  factor  Rj  by  which  the  peak  attachment  cross  section  0^ 

is  altered  on  replacing  H2  by  one  of  its  isotopes  X, 
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Rt  =  ^pcak(H  ,  v  -  0,  J  -  0)/CTpCaj;(X  ,  v  =  0,  J  —  0). 

(b)  For  a  given  isotope  X,  the  factor  Rv  by  which  the  peak  attachment  cross  section  Cp^  is 
altered  on  exciting  the  molecule  vibrationally  from  v  =  0  to  v  =  1 , 

Rv  =  cpcak(X~.  v  =  1,  J  =  0)/apcak(X‘,  v  =  0,  J  =  0). 

(c)  For  a  given  isotope  X,  the  factor  Rj  by  which  the  peak  attachment  cross  section  o!)cak  is  altered 
on  exciting  the  molecule  rotationaliy  from  J  =  0  to  J  =  10, 


Rj  =  '-fycakC^  ,  v  =  0,  J  -  10)/CTpcak(X  ,  v  —  0,  J  =  0). 


The  cross  sections  for  electron  attachment  to  various  isotopes  of  molecular  hydrogen  are  shown,  as  a 
function  of  electron  energy,  in  Figure  2.  For  each  isotope  three  cross  sections  are  shown  in  this  Figure. 
The  first  set  of  six  curves  provides  the  cross  sections  for  electron  attachment  to  the  six  isotopes  in  the 
lowest  rovibrational  levels,  namely,  (v,J)  =  (0,0).  The  next  set  of  six  curves  in  Figure  2  shows  the  cross 
sections  for  attachment  to  the  vibrationaily  excited  isotopes  in  the  levels  (1,0).  The  third  set  of  six  curves 
depicts  cross  sections  for  attachment  to  rotationaliy  (0,J)  excited  molecular  isotopes;  the  value  of  rotauonal 
quantum  number  J  is  chosen,  for  each  isotope,  in  such  a  manner  that  the  rotational  excitation  energy  is 
roughly  equal  to  the  vibrational  excitation  energy  of  the  (1,0)  level.  For  each  of  the  isotopes,  the 
rovibrational  excitation  is  seen  to  enhance  the  attachment  cross  section.  Furthermore,  vibrational  excitation 
is  more  effective  in  enhancing  the  cross  section  than  rotational  excitation.  The  fundamental  reason  for  the 
enhancement  of  the  electron  attachment  cross  sections  is  an  increase  in  the  range  of  intemucicar  separations 
over  which  the  electron  can  be  captured  for  anion  formation  once  the  molecule  is  internally  excited.  This 
increase  in  range  occurs  due  to  an  increased  vibrational  amplitude  during  vibrational  excitation  and  due  to 
centrifugal  stretching  during  rotational  excitation. 

Table  1.  Various  factors  indicating  the  enhancement  of  the  peak  cross  section  for  attachment  to 
rovibrationally  excited  H2  and  its  isotopes. 


Isotope 

R1 

Rv 

R] 

h2 

1.0 

33.8 

15.9 

HD 

10.3 

39.4 

11.8 

HT 

29.6 

42.1 

10.5 

D2 

527 

48.9 

8.04 

DT 

4064 

54.2 

6.89 

To 

65217 

60.9 

5.72 

Various  enhancement  factors,  Rj,  Ry  and  Rj.  depend  upon  the  mass  (M)  of  the  isotope  and 
empirically  it  is  determined  that  the  mass  dependences  of  these  factors  are: 

Rj  exp  (- constant  *M 1/2)  (la) 
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Electron  Energy  (eV) 

Figure  2.  Cross  sections  for  dissociative  electron  attachment  to  rovibrauonally  excited  H2  and  its  isotopes. 

The  change  in  the  peak  attachment  cross  section  on  isotope  substitution,  namely  RI(  is  referred  to  as 
the  isotope  effect.  The  isotope  effect  is  most  dramatic  for  lighter  molecules  like  hydrogen  because  of  the 
greater  relative  change  in  the  reduced  mass  of  the  nuclei  on  isotope  substitution.  The  effect  could  be 

completely  understood  in  terms  of  a  semiclassical  analysis7,8.  In  such  an  analysis  the  cross  section  for 
dissociative  electron  attachment,  OoA,  can  be  written  as  a  product  of  two  factors: 

Th»  r  f  CTDA  =  °cap 

1  ne  nrst  ractor,  ocap,  is  interpreted  as  the  cross  section  for  the  formation  of  the  resonant  anion  state 
b^ the  capture  of  the  incident  electron.  The  second  factor,  S,  is  interpreted  as  the  probability  that  the 
separation  of  the  nuclei  in  the  anion  state  increases  from  the  capture  radius  Rc  to  the  stabilization  radius  Rs 
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without  electron  autodetachment  to  assure  the  process  of  dissociative  electron  attachment  to  occur.  The 
factor  S  is  referred  to  as  the  survival  probability  and.  in  the  semiclassical  analysis,  is  given  by: 


where  F(R)  is  the  width  of  the  resonant  state  and  v(R)  is  the  speed  with  which  the  nuclei  arc  moving  apart 
at  the  intemuclear  separation  R.  Now,  in  the  semiclassical  expression,  Eq.  (2),  the  survival  probability  S 
is  a  strongly  mass-dependent  quandty.  It  can  roughly  be  approximated  by  exp  (-  Pt/fi)  where  x,  the  time 
taken  by  nuclei  to  separate  from  the  capture  radius  Re  to  the  stabilization  radius  Rs,  is  directly  proportional, 

by  simple  kinematical  considerations,  to  M1/2.  Thus  nuclei  of  heavier  isotopes,  taking  a  longer  time  than 
'  nuclei  of  Hi  to  separate  out  to  Rs,  experience  stiffer  competition  from  the  process  of  electron 
autodciachment  which  reduces  the  probability  of  dissociative  attachment. 

As  an  example  of  the  isotope  effect  one  should  note,  using  the  values  of  Rj  in  Table  1,  that  the  value 
of  the  ratio 

ln  (~q peak  (TJ  /  °penk  (H2)~| 

•n[OpeakP2)/V<>k(H2)l  (3) 


is  in(65217)  /  ln(527)  =  1.7688.  The  values  of  the  peak  attachment  cross  sections  are  obtained  by  fully 
quantum  mechanical  calculations.  Now,  if  the  isotope  effect  of  the  semiclassical  Eq.  (la)  is  valid,  this 
ratio  should  be 


mTj 


m^2-m 


V?  - 1 
fz- 1 


(4) 


whose  value  is  1.7673.  The  two  values  are  quite  close,  any  difference  being  due  to  (he  mass  dependence 
of  the  vibrational  frequencies  of  the  various  isotopes. 

For  the  purpose  of  having  a  panicle  beam  with  high  energy  deposition  characteristics,  one  looks  for 
beam  panicles  as  massive  as  possible.  To  investigate  die  production  of  massive  negative  ions  we  have 
obtained  the  factors  by  which  the  cross  sections  (and,  therefore,  the  rates)  for  the  production  of  Li*  from 
Li2  as  well  as  of  Cl'  from  HC1  are  enhanced  on  vibrationaUy  exciting  the  neutral  molecules.  These 

enhancement  factors  are  approximately  8  for  Li*  production9  and  40  for  Cl*  production10,11  when  the 
molecules  Li2  and  HCl  are  initially  excited  fror, .  die  v=  0  to  the  v  =  1  level.  The  corresponding 
enhancement  factor  for  H*  production  from  H2  is  calculated  to  be  approximately  34.  Furthermore,  the 
maximum  rate  of  Li*  production  from  Li2  is  calculated9  (as  well  as  observed12)  to  be  about  10"8  cm3  s*1, 
which  is  same  as  the  rate  of  production6  of  H*  from  H2.  This  is  not  surprising  since  Li2  is  isovalent  with 
H2. 


4.  ru'TURE  POSSlbiLITlhS 

Calculations  of  the  rate  of  negative  ion  production  require  a  knowledge  of  the  electron  energy 
distribution  in  the  ion  source.  In  our  previous  calculations4,6  a  simple  Maxwellian  distribution  for  electron 
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energies  was  assumed .  Using  this  electron  energy  distribution  the  maximum  possible  rate  for  dissociative 
electron  attachment  to  Hj  was  calculated  to  be  about  10'8  cm3  s'1.  Using  the  same  distribution  function 
the  maximum  kinetic  energy  carried  by  the  negative  ion  H'  formed  by  dissociative  electron  attachment  to 
H2  was  estimated  to  be  less  than  0.5  eV.  We  have  recently  developed13  a  novel  procedure  for  obtaining 
the  time-dependent  behavior  of  the  electron  swarms  in  a  gas  mixture.  This  numerically  stable  algorithm 
generates  an  exact,  time-dependent  solution  of  the  Boltzmann  equation  for  charged  particle  swarms  in  a 
dilute  gas  and  uniform  electric  field  and,  thereby,  provides  a  temporal  evolution  of  the  electron  energy 
distribution  in  the  gas.  We  have  already  tested  this  procedure  for  generating  the  time-dependent  behavior 
of  various  parameters  like  the  average  electron  energy,  average  drift  velocity,  ionization  rate  etc.  for 
electron  swarms  in  gaseous  neon  and  argon.  We  are  currently  utilizing  the  same  numerical  procedure  for 
obtaining  the  time-dependent  behavior  of  parameters  used  for  diagnostics  purposes  in  a  source  containing 
a  mixture  of  atomic  and  molecular  hydrogen. 
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ABSTRACT 

Cross  sections  for  the  vibranonal  excitation  of  molecular  hydrogen  and  its  five 
heavier  isotopes  by  the  impact  of  low  energy  electrons  are  calculated  using  both  the 
local  and  the  nonlocal  versions  of  the  resonance  model.  It  is  demonstrated  that,  for  a 
given  incident  electron  energy,  the  cross  sections  for  the  vibrational  excitanon  of 
heavier  isotopes  can  be  obtained  from  those  of  molecular  hydrogen  by  a  simple  scaling 
procedure.  It  is  seen  that  the  scaling  law  holds  for  cross  sections  with  values  tanging 
over  eight  orders  of  magnitude. 

The  vibrational  excitanon  of  a  molecule  (such  as  H2 )  occurs  via  the  formanon 
of  an  intermediate  resonant  state  (H2‘ )  which  can  cither  autodctach  the  temporarily 
bound  electron  (and,  then,  leaves  behind  a  vibiauonally  excited  molecule)  or  can 
dissociate  into  H  +  IT  (that  is.  leads  to  dissociative  electron  attachment). 
Schemancally,  these  two  complementary  processes  are: 

r~  e’  +  Hj(vf,  Jf) 

e‘  +  H2  (v, ,  J, )  -  H2-  ( %*  or  2Zg+ )  -l 

V-  H  +  IT. 

The  nuclear  wave  function  q  of  the  resonant  state,  in  the  most  general  thcorcncal 
description,  satisfies1  a  iniegrodiffcremial  equation  with  nonlocal  potennals: 

[Tn(R)  +  V  (R).EK(R)=.  V(E-EvR)x,i(R) 

•  JiR'Zz'WXvfR)  {  P/dt  V'(E  Te^T ~  - 1  «  V»(E  -  Ev .  K)  V(E  •  Ev .  R)  }  \  (R7  . 

(1) 

Here  xv  arc  the  vibrational  bound  state  wave  functions  of  the  neutral  molecule,  TN  is 

the  nuclear  hnenc  energy  operator  and  V (R)  is  the  potential  curve  of  the  resonant 
anion  state.  The  sum  on  the  right  hand  side  includes  all  energetically  open  channels 
The  matrix  elements  V(e,R)  couple  the  discrete  resonant  state  with  the  background 
continuum;  this  coupling  leads  to  an  energy-shift  of  the  resonance  (the  principal  pan 
integral  on  the  right  hand  side)  and  provides  a  width  to  the  resonance  (the  imaginary 
pan  on  the  right  hand  side)  which  determines  the  lifetime  of  the  resonance.  The  total 
cross  section  for  vibrational  excitation,  in  atomic  units,  is1  (assuming  momentum 
normalized  continuum  functions), 

4 

,  ..  16 r  k'f  ,,.2 

<j(v  —  v  )  =  — - — j  die  IT(v  —  v  )l  . 
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threshold.  In  any  case  the  nonlocal  effects  are  found  to  be  small.  All  the  excitation 
cross  sections  sansfy  a  useful  isotope  scaling  law  such  that  the  vibrational  cxcitauon 
cross  sections  for  the  heavier  isotopes  of  molecular  hydrogen  can  be  obtained,  within 
a  factor  of  two,  from  the  corresponding  cross  secuons  for  Hj. 
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Tabic  J  Thermodynamics  of  LixHy  clusters 


Fig.  1  UH/LiH“  potential  energy  curves. 


Trigonal 

pyramid 

(C3v) 


_ _ i 

2.58  U 

•table 


unstable  as  a 
C3V  structure 


Kite 

structure 

(C2v) 


s»7i 


oy  i 
Y»i 
31 


Fig.4  LijH/LijH-  structures. 


unbound  by  0.24  eV 
relative  to  neutral  (C2v ) 


Linear 

(Ccov) 


Rhomboid 

(°2h) 


V  -  structure 
(C2v) 


u^L8H!±5u!flH 

M-  =  13.7  debyo 
sacfdte  pt.:  1  imag.  freq. 
,  M 


stable  sructure 


unstable  U2  +  H2 


Lll!_0H!l3u!flH 

unbound  by  0.040  eV 
relative  to  neutral  (Co>v) 

.LI 


unhnunH  hu  o  nc  Atr 

- —  —  —  «mV\t  W  T 

relative  to  neutral  (D2h) 


u-522-ucdl0' 

unbound  by  1.50  eV 
relative  to  neutral  (D2h) 


Fig.S  LiiHi/Li2^h~  structures. 


150  Gas  Phase  At  kali -Hydrogen  Interactions 


"Relativistic  approach  for  e  sca“c"n.fq?^m  arg0n  ’ 


0X33  Relativistic  Approach  for  e"  Scattering  froa 
Argon.  SOUfARTtr:  RaUAR,  Georgia^  Shtt  State  U.,  andJT 
KT^AOEHRA,  Vayne  State^T^--  ClTTerenTTal  and 
Integrated  crS3s  s«7I3ns  and  various  polarization 
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parameters,  such  os  polarization  (P)  or  the  the  Sheraon 
function  (S)  and  parameters  T  and  U,  for  electron  and 
positron  scattering  from  argon  are  calculated  using  the 
relativistic  Dirac  equation.  It  contains  the  spin-orbit 
interaction.  The  real  part  of  the  projectile- garget 
interaction  is  represented  by  n  aodel  potential  that 
includes  static  potential  (repulsive)  and  the  Buckinghaa 
type  eodel  polarization  potential  (attractive)  for  the 
positron  scattering  and  static  potential  (attractive), 
the  sa*e  polarization  potential  (attractive)  and  electron 
exchange  potential  for  electron  scattering  froa  argon. 
The  phase  shifts  vilh  large  angular  nonenta  it  are 
calculated  by  using  the  Born  approxination  The 
polarization  paraaeter  P  for  electron  scattering  is  found 
to  be  in  good  agreeaent  vith  the  available  calculated  and 
oeasured  values.  A  fev  different  aodels  of  the  absorp¬ 
tion  potential  for  the  inelastic  processes  are  used  to 
calculate  the  elastic  differential  and  integrated  as  veil 
as  the  total  integrated  cross  sections  for  positron  scat¬ 
tering  fron  argon.  It  is  noticed  that  even  though  the 
calculated  total  integrated  cross  sections  agree  reason¬ 
ably  veil  vith  the  aeasjred  values,  the  differential 
cross  section  curves  shov  features  different  fron  those 
aeasured  for  the  positron  scattering  froa  argon. 

*  Support  of  AFOSR  is  gratefully  acknovledged. 

•S.N.  Nahor  and  J.H.  Vadehra,  Phys.  Rev.  A  35,  2051(1987). 


"Formation  of  ground  and  excited  states  of  antihydrogen", 

Bull.  Am.  Phys.  Soc.  21,  991  (1988) 

presented  at  the  1988  annual  meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics, 
Baltimore,  Maryland,  April  18-20, 1988. 


EX  5  E££iiii£5,l.STN.si§!!^NAmE  sfltM  Stfirtr.  aK3 

antihydroRen. 5ULTXNA  _»•  • - lrr(.r»ntlri  and 

r-R^-tteHRA,  SMSS  «as  th“Ut£ot..,lon  of 

integrated  orcaa  sections  ^for  lntecB#ilate  energ* 

(lO^SOO^keVl^antlprotons  on^prsltronlua^ere^calculate^ 
carr?ed^eoufrStor°the*for«ation  of ^ntlhydrogen  In  ground 

&S2SL. 

for  the  forsation  of  ‘  exclted  electronic 

„h=n  the  pos  tronluo  i  In  «|ot  the  capture 

u::r  » «ic«i.t«  «■..  *.«i  (that  i5> 
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all  states  added  together)  integrated  cross  sections. 
The  cross  sections  for  the  fornatlon  of  antihydrogen 
presented  here  ore  obtained  troo  those  for  the  formation 
of  positronlua  by  the  iapact  of  positrons  on  hydrogen 
ntonj  by  using  the  charge  invariance  and  the  principle  oi 
hnlnnce. 


•Support  of  AFOSR  is  gratefully  acknowledge. 


"An  exact  numerical  solution  of  the  time  dependent  Boltzmann  equation", 
Bull.  Am.  Phys.  Soc.  M,  295  (1989) 

presented  at  the  41st  Annual  Gaseous  Electronics  Conference,  Minneapolis,  Minnesota,  October 

18-21,  1988. 


£-2  An  Exact  .Numerical  Solution  of  the  Tlgc_Dcpendent 
WnlrytcAnn  Equation . *  P.J.  DRALLQS  And  J.M.  UADEHRA,  Dept . 
of  Physics.  Wavne  State  Unlv.--  An  exact,  time  dependent 
numerical  solution  of  the  Boltzmann  equation  for  charged 
particle  swarms  in  a  dilute  gas  and  unifortr  electric 
field  is  presented  in  detail.  The  method  incorporates 
the  full  anisotropy  of  both  the  velocity  distribution 
function  and  the  collision  cross  sections  as  it  does  not 
involve  any  tern  expansions.  An  exact  analysis  of  the 
collision  terms  is  described  and  conditions  for  numerical 
stability  of  the  solution  are  discussed.  Results  are 
presented  for  electron  and  positron  swarms  in  gaseous 
Neon  and  Argon  and  in  some  model  gasos  at  various  values 
of  E/N. 

*This  work  was  supported  by  the  Air  Force  Office  of 
Scientific  Research  through  Grant  No.  AF0SR*87*03&2. 
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"Energetics  of  negative  ion  formation  via  dissociative  attachment  of  e  +  LiH", 
Bull.  Am.  Phys.  Soc.  M.  1401  (1989) 

presented  at  the  1989  annual  meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics, 

Windsor,  Ontario,  May  17-19, 1989. 


2:00 

FC6  Energetics  of  Negative  Ion  Formation  vil-Dissociativc 
Attachtrvm  of  c  +  l.iH*.  H.  H.  Mi  chela,  UTRC.  and  J.  M-Wadehra, 
Wavnr  Starr  U. —  The  tarnation  of  H"  and  Li"  by  dissociative 
attachment  (D A)  of  e  +  H2  and  e  +  Ltj,  respectively,  ij  now  well 
characterized  both  experimentally  and  theoretically.  the  role  of  UH 
(or  CsH),  which  could  be  tamed  from  seeding  an  alkali  into  a 
hydrogen  plasma,  is  presently  not  well  understood,  but  the  addition  of 
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an  alkali  appears  to  enhance  the  H"  production  rate.  This  observation 
is  interesting  in  light  of  the  study  by  Gauyacq,  et  al>,  which  indicates 
that  charge  transfer  and  collisional  detachment  processes  should 
reduce  H"  production  in  Na  seeded  plasmas.  We  have  analyzed  the 
electron  attachment  to  LiH  in  terms  of  calculated  potential  energy 
curves.  In  agreement  with  previous  studies,  wc  find  that  the  ground 
state  of  LIH"  is  thermodynamically  bound  relative  to  LiH  +  e,  with  a 
calculated  electron  affinity  of  ~  0.3  eV.  The  first  excited  state  of 
LiH">  which  asymptotically  correlates  to  Li"  +  H,  exhibits  repulsive 
behavior  in  the  region  3.0  £  R  £  6.0  A.  Based  on  these  preliminary 
studies,  DA  of  e  +  UH  will  yield  Li"  ions  for  low-energy  collisions 
H"  10ns  are  not  formed  by  DA  of  e  +  LiH  but  may  be  formed  by 
energetically  allowed  electron  capture  into  high  vibrational  states  of 
LiH  which  lie  above  the  U  +  H"  asymptotic  limit. 

•Supported  in  part  by  AFOSR  under  Contract  F49620'88*C*0019  and 
Grant  AFOSR-87-0342. 

*J.  P.  Gauyacq,  et  al,  Phys.  Rev.  A,  22,  2284  (1988) 


"High-field  time-dependent  positron  velocity  distribution  functions", 

Bull.  Am.  Pnys.  Soc.  M.  409  (1989) 

presented  at  the  1989  annual  meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics, 

Windsor,  Ontario,  May  17-19, 1989. 


FX  46  anhjagld  Tice -Dependent  Poamsit..  VelgGltz 
niit-rlhuHon  Functions.  P.J.  DRALLOS  ^n^  J,H.  VADEHRA, 
Uavne  Stare  has  been  reported  '  that  for  a  noon 

gaa  dona Icy  N  -  1  anagat,  cho  positron  avarn  paraaeter 
Z  „  varies  only  »t  vary  low  flolda  (up  Co  E  -  5  V/cn) 
«ndfchon  beconet  alnost  lndopondonc  of  E  (up  co  E  -  50 
V/cn)  with  a  value  of  abouc  6.3.  In  the  proaonc  calcula- 
clona,  a  nunorlcal  cachnlquo  for  evaluating  the  exact 
tine-dependent  behavior  of  elecc-on  veloclcy  distribution 
functions1  has  boon  adapted  to  the  case  of  poaltron 
svarns  In  gaaooua  neon.  In  these  calculations  an  F/H  of 
0.9  Td  (E  -  240  V/cn  and  H  -  1  anagat)  vas  used,  With 
this  technique,  the  tine  dependence  of  the  average 
energy,  annihilation  rate,  and  Z  „  were  obtained.  The 
equlllbrlun  value  for  Z  {{  under  ttieee  eondltlona  vat 
about  6.3.  This  result  extends  the  range  of  electric 
field,  for  vhlch  Zeff  is  a  constant  to  E  -  240  V/cn. 

*The  support  of  the  AF0SR  Is  gratefully  acknowledged. 

1.  p.s.  Grover,  J,  Fhya.  B  10,  2269  (1977) 

2.  K.V.  Slnha  and  P.S.  Grover.  Phys.  Rev,  A  35.  3309 
(1987). 

3.  P.J.  Dralloa  and  J.M.  Vadehra.  J.  Appl.  Phys.  63.  5601 
(1988) 


"Energetics  of  neutral  and  negative  ions  of  LimHn  clusters", 

Bull.  Am.  Phys.  Soc.  21,  1 177  (1990) 

presented  at  the  1990  annual  meeting  of  the  Division  of  Atomic,  Molecular  and  Optical  Physics, 
(  Monterey,  California,  May  21-23, 1990. 
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DB  9  Energetics  of  Neutral.Aa<l  IScgaiiyc  Iona  of  LimHnXlmisis*.  H  H 

Michels.  UTRC.andJ  M  Wadchra.  Wavne  State  U. —  The  formation  of  H~(Lr) 
by  dissociauve  attachment  (DA)  of  e  *  H2(Li2)  is  thought  to  be  the  dominant 
volume  process  in  discharge  type  negative  ton  sources  The  role  of 
L^iyCSnll,,)  molecules,  which  could  be  formed  from  seeding  an  allcali  into  a 
hydrogen  plasma,  is  presently  not  well  understood,  but  the  addition  of  an  alkali 
such  as  Cs  appears  to  enhance  the  H"  production  rate  This  observation  is 
interesting  m  light  of  the  study  by  Gauyacq.  ct  al1.  which  indicates  that  charge 
transfer  and  colltsional  detachment  processes  should  reduce  H"  production  in  Na 
seeded  plasmas.  In  order  to  examine  the  structure  and  stability  of  L^H,,  clusters 
and  their  anions,  ab  mtio  calculations  were  earned  out  for  several  species  at  the 
MP2  level  of  theory  The  basis  set  chosen  was  the  Gaussian  6*3 11G  triple  split* 
valence  set.  augmented  by  d-poianzauon  functions  for  la  and  p- polarization  for  H 
and  diffuse  functions  to  baler  describe  the  negative  ion  charge  distribuuons  Our 
studies  to  date  indicate  that  several  Lx^H^  (and  by  analogy.  clusters  axe 

thermodynamically  stable.  In  particular,  the  LijHj  species,  as  a  C2v  structure, 
may  be  an  important  component  of  alkali-hydrogen  mixtures  This  species  can 
dissociauvely  auach  an  electron  to  form  L»2H  ♦  H’  for  E^u  £  2  0  cV  The  Lhll 
(C2v)  species  should  also  exhibit  DA  to  form  Lij  -f  H“.  but  the  concentn.uon  ol 
this  molecule  will  be  lower  than  that  of  the  more  stable  Li2H2  cluster 


•Supported  in  part  by  AFOSR  under  Contract  F49620-89-C-O019  and  Grant 
AFOSR-87-0342 

OP  Gauyacq.etal  Phys  Rev  A.  2$.  2284  (1988) 
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